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ABSTRACT 


Theoretical  analyses  of  flow  with  rotating  passages 
or  turbo-machinery  have  become  necessary  for  the  proper 
design  of  turbo-pump  elements  in  liquid  fuel  boosters  and 
power  conversion  units.  One  such  theory  is  presented,  and 
a  numerical  solution  derived.  This  theory  develops  a  meth¬ 
od  for  the  analysis  of  steady,  inviscid,  adiabatic  flow 
through  arbitrary  rotors.  A  detailed  analysis  in  a  merid¬ 
ional  plane  is  given,  assuming  axial  symmetry.  A  simplified 
approach  to  the  blade  to  blade  solution  is  also  presented. 
The  merits  of  these  theories  are  compared  with  other  pro¬ 
posed  methods.  The  inverse,  or  design,  approach  is  con¬ 
sidered,  and  found  to  be  unnecessary. 

A  numerical  solution  for  incompressible  flow  is  derived 
and  applied  to  the  flow  solution  in  the  impeller  of  a  mixed 
flow  compressor  with  backwards-bent  blades  of  arbitrary 
shape.  Meridional  streamlines  and  relative  velocity  dis¬ 
tributions  are  progressively  calculated  on  a  CDC  1604  com¬ 
puter,  using  FORTRAN  program  language.  Data  are  measured 
from  a  detailed  presentation  of  the  blade  shape  in  a  merid¬ 
ional  plane.  Blade  to  blade  relative  velocity  distributions 
are  calculated  from  the  meridional  plane  analysis. 

It  is  concluded  that  the  results  completely  define  the 
flow  and  are  sufficiently  accurate  for  engineering  applica¬ 
tions.  Validation  is  based  upon  the  reliability  of  the  the¬ 
ory,  and  upon  comparisons  with  results  of  other  methods. 
Extensions  of  the  scope  of  this  approach  are  recommended, 
which  include  the  compressible  solution  and  the  solution  of 
flows  in  unbladed  passages. 

The  writer  wished  to  express  his  appreciation  for  the 
assistance  and  encouragement  given  him  by  Professor  M.  H. 
Vavra  of  the  U.  S.  Naval  Postgraduate  School.  His  guidance 
in  the  development  of  this  solution  was  invaluable. 
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A  METHOD  FOR  THREE  -  DIMENSIONAL 


FLOW  ANALYSIS  IN  A  ROTOR 
USING  A  HIGH  SPEED  DIGITAL  COMPUTER 

INTRODUCTION 

There  is  an  increasing  demand  for  compact,  high  speed, 
high  power  output  turbomachines  for  use  in  liquid  fuel  roc¬ 
ket  engines  and  power  conversion  units.  These  units  are 
being  designed  to  work  with  such  media  as  cryogenics  and 
liquid  metals.  The  proper  design  of  such  machinery  largely 
depends  upon  the  accuracy  of  the  theoretical  analysis,  used 
prior  to  fabrication.  The  method  used  should  reasonably 
predict  aerodynamic  forces  and  tendencies  toward  flow  sep¬ 
aration  and  cavitation. 

The  more  simplified  approaches,  using  conditions  ahead 
of  and  after  the  rotor,  ignore  conditions  within  the  rotor 
passage  itself,  and  therefore,  are  unsatisfactory.  Poten¬ 
tial  flow  analyses  are  useful  for  stationary  cascades,  but 
cannot  be  accurately  applied  to  flows  in  rotating  passages. 
The  need  for  more  precise  methods,  required  for  the  proper 
analysis  of  modern  rotor  designs,  was  anticipated.  Several 
three-dimensional  theories  have  been  postulated.  In  gen¬ 
eral,  these  theories  are  quite  complex  and  are  difficult 
to  apply. 

One  such  theory  was  developed  by  the  NACA's  Lewis 
Laboratory  in  the  early  1950's.  The  three  dimensional 
problem  was  reduced  to  an  iterative  process  between  two- 
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dimensional  solutions  on  hub  to  tip,  and  blade  to  blade 
stream  surfaces.  Ref.  1  presents  an  analysis  on  a  merid¬ 
ional,  or  hub  to  tip,  surface.  A  similar  approach  is  de¬ 
veloped  in  Ref.  2.  A  blade  to  blade  analysis  is  described 
in  Ref.  3.  These  treatments  were  combined  into  a  general 
theory  in  Ref.  4.  Subsequent  efforts  were  directed  towards 
applications  of  this  theory  to  specific  examples.  The  com¬ 
plex  nature  of  the  theory  required  the  introduction  of  cer¬ 
tain  simplifying  assumptions.  The  development  of  the  high 
speed  digital  computer  made  solutions  more  practical,  as 
demonstrated  in  Refs.  5  through  7. 

Prior  to  the  NACA’s  work,  a  method  of  solution  in  a 
meridional  plane,  similar  to  that  of  Ref.  1,  was  developed 
by  Meyer  in  Ref.  8.  This  method  uses  an  iterative,  graph¬ 
ical  process  to  solve  two  simultaneous,  linear  differential 
equations  along  the  characteristics  of  the  equations.  The 
equations  are  derived  from  the  Eulerian  equations  of  motion, 
and  the  characteristics  are  defined  by  the  geometry  of  the 
flow  channel.  In  Ref.  9,  Vavra  reorganized  Meyer’s  scheme 
into  a  more  general  theory,  reducing  the  equations  of  mo¬ 
tion  to  one  non-linear  equation  which  is  also  solved  along 
characteristics.  In  addition,  a  simplified  blade-to-blade 
analysis  was  developed. 

The  purpose  of  this  thesis  is  to  transform  the  theo¬ 
retical  development  of  Ref.  9  into  a  method  of  solution 
for  arbitrary  rotors,  using  the  CDC  1604  digital  computer. 

A  completely  contained  computer  solution  was  presumed  to 
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be  quite  long.  It  was  anticipated  that  the  amount  of  in¬ 
put  data,  required  for  acceptable  accuracy,  would  exceed 
computer  capacity  and  that  long  preparation  and  computer 
run  time  would  be  required.  It  was  decided  to  use  a  ser¬ 
ies  of  short  computations,  which  could  be  repeated  as  often, 
or  at  such  intervals,  as  required  to  obtain  the  desired 
accuracy. 

The  theoretical  development  may  be  applied  to  both 
the  design  and  inverse  problems.  The  major  portion  of  this 
thesis  is  devoted  to  the  latter  application j  the  analysis 
of  flow  in  a  given  rotor.  A  design  attempt  is  made  to  pro¬ 
vide  a  physical  model  for  the  evaluation  of  the  analytical 
methods.  This  attempt  is  not  completed  because  of  time 
considerations,  and  an  actual  rotor  is  used  for  the  model. 
This  impeller  is  a  part  of  a  compressor  test  rig  located 
at  the  U.  S.  Naval  Postgraduate  School.  Details  of  this 
machine  are  enumerated  in  Ref.  10.  The  rotor  is  of  the 
mixed  flow  type,  with  non-radial  blades.  The  blades  have 
a  constant  thickness,  and  are  of  the  deloaded  type,  with 
a  reversal  in  curvature  over  the  rear  section.  This  con¬ 
figuration  presents  an  arbitrary  design  which  does  not 
conform  to  standard  impeller  types. 

The  methods  developed  in  this  thesis  are  applied  to 
this  model,  using  simplified  initial  conditions,  which  do 
not  reflect  actual  operating  conditions,  but  are  within 
the  compressor " s  operating  range.  No  correlation  is  made 
between  theoretical  and  test  results,  however  theoretical 
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results  are  presented  for  validation.  Both  the  theory 
and  the  results  of  this  analysis  are  compared  with  those 
of  Refs.  1  through  7,  to  substantiate  the  applicability 
of  this  method. 
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II.  THEORY 


A  complete  derivation  of  the  theory  is  presented  in 
Ref.  9.  This  section  contains  those  equations  and  devel¬ 
opments  that  are  considered  necessary  for  clarity  and 
continuity,  in  following  the  subsequent  transformation 
from  theory  to  numerical  methods. 


A.  Meridional  Plane  Analysis 

1.  The  following  assumptions  are  made  to  establish  a 
model  for  the  hub  to  tip  analysis.  A  more  thorough  treat¬ 
ment  of  these  assumptions  is  conducted  in  Section  IV. 

a.  The  rotor  cascade  contains  an  infinite  number  of 
infinitely  thin  blades.  Thus,  stream  surfaces  and 
fluid  motion  are  axi symmetric. 

b.  Flow  is  inviscid,  steady,  and  isentropic. 

c.  Entropy  changes  due  to  discontinuities  at  the  lead¬ 
ing  and  trailing  edges  are  acknowledged  but  ignored. 

d.  Flow  is  incompressible. 

2.  The  solution  of  the  equation  of  motion  is  based  upon 
complete  definition  of  the  geometry  of  the  blade  surface. 

A  set  of  meridional  streamlines,  m,  are  assumed  and  the 
normals,  n,  are  drawn,  establishing  an  orthogonal,  axi sym¬ 
metric  coordinate  system.  This  system  is  corrected  by 
successive  approximations.  The  blade  surface  is  repre¬ 
sented  by  the  circular  projection  on  the  meridional  plane 
of  the  lines  of  intersection  of  the  blade  surface  with 
planes  ©=constant.  0  is  the  angle  measured  in  the 
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peripheral  direction.  The  ©  constant  planes  are  planes 
extending  in  a  radial  direction  from  the  axis  of  rotation, 
perpendicular  to  that  axis.  The  peripheral  angle,  ©  ,  is 
referenced  to  an  arbitrary  point  on  the  blade  surface, 
usually  on  the  leading  edge.  The  ©  =constant  curves  are 
shown  in  Fig.  1.  Assumed  coordinate  systems  are  illus¬ 
trated  in  Fig.  4. 

The  angles  defined  by  this  family  of  curves,  at  an 
arbitrary  point,  P,  are  shown  in  Figs.  2  and  3.  The  angle 
6  is  the  inclination  of  a  ©  =constant  curve  from  the  nor¬ 
mal.  The  angle  A  is  the  inclination  of  a  meridional 
streamline  from  the  axial  direction,  Z  .  The  angle  (3 
is  the  flow  angle  of  the  relative  velocity,  W,  on  the 
stream  surface. 

The  deviation  of  a  blade  section,  £  ,  from  the  radial 

direction,  R,  is: 

tar\  ©  =  tav\  (5  Sm  ( 2 ) 

COb  <5" 

Thus,  the  system  of  streamlines  and  0  =constant  curves 
are  sufficient  to  completely  define  the  blade  shape. 

The  relative  velocity  on  a  streamline  may  be  ex¬ 
pressed  vectorially  as: 

V)  «  U  VJ^ior,  (A  **  ''fc  ^  3  j 

where:  u  =  unit  vector  in  peripheral  direction 
m  ss  unit  vector  in  stream  direction 
The  flow  problem  is  reduced  to  a  solution  for  the  merid¬ 
ional  component  of  the  relative  velocity,  Wm# 
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3.  The  Euler ian  equation  of  motion  is: 


&5  +  VH  "  V  x  (^^V)  +  Tvs  +  1 

For  steady ,  isentropic,  relative  flows: 

VHr  =  Vjx(v*W  +  2<2>)  +  ? 


where. 

Ho  »  h  V  W2  _  coxR2 
«  1Z 


(4) 

(5) 

(6) 


For  inviscid  flow,  the  friction  force,  f,  is  zero.  How¬ 
ever,  the  effect  of  infinitesimal  pressure  changes  across 
an  infinite  number  of  blades  is  accounted  for  by  intro¬ 
ducing  the  blade  force,  F0,  which  is  normal  to  a  blade 
element.  Thus: 

VHr  =  wx  (  V*  w  +  2  to')-*-  Fb  (7 


Eq.  (7)  is  reduced  to  its  scalar  compenents,  using: 

fi!  - 0  (s.) 

Fb  -  Fy  (0- 4-f\-ta'ofb-ta>o5N)  (8b) 

where:  ft  =  unit  vector  in  normal  direction 
The  u'  component  yields: 

RFu  =  £>■  (RVVu  ^  00  (9a) 

Z)  rr\ 


The  m  component,  with  Eq.  (9a)  gives: 

=  o 

2>tY| 


(9b) 


This  relation  only  holds  within  the  rotor  for  the  assumed 
theoretical  model.  Flows  outside  the  rotor  are  unaffected 
by  such  a  model. 
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For  all  but  design  conditions,  there  is  a  flow  dis¬ 
continuity  at  the  leading  edge,  since  flows  ahead  of  the 
rotor  will  not  meet  the  rotor  at  blade  entry  angles.  It 
has  been  assumed  that  entropy  changes  caused  by  these  dis¬ 
continuities  are  neglected  by  ignoring  the  changes  in  total 
enthalpy  across  the  leading  edge.  Conditions  ahead  of  the 
rotor  are  derived  from  Eq.  (4)  for  steady,  isentropic  flow. 

c)H  _  Vy  —  o  (10) 

b™  R  ‘b'PcN 

A  relation  between  conditions  on  a  streamline  ahead  of  the 
rotor  and  within  the  rotor  is  established  as: 

=  ^  R  u  (  Hu  -i-  <^o  R  j  X1  j 

where  e  and  L  denote  entrance  and  leading  edge  stations. 

The  change  of  HR  along  a  normal  is: 

—  C>Ue  b  As.  —  CO  Hu  ^  ^  cl  L 

bA  bAe.  c>A  <dvi  (12) 


The  n  component  of  Eq.  (7)  becomes: 

w ^  -V  +  Wu  c)(R\NJu  -4-  c*Rz') 

£>A  ^  bA 

(9c) 

.  Ns/m  NRWu-v  doRx>)  -Lahfet-aA  S  =5  \<(Ab 

R  bvn  v 


where  k  is  the  curvature  of  the  meridional  streamline, 
m 

From  Fig.  3: 


(13) 
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Eq.  (9c)  is  reduced  to  a  relation  in  W  s 

in 

£ 

^  (j  +  -La -v-  (i_a>n2A  t/an15>) 

■&Y\  v  ' 

+  Wl (2 +  2tp ft  ^ +  far.  5  | tR  t**)  fffj  ^  (14) 

•*Wm  (M-oo+a«(i[  |&  +  tar,  S  ) 

=  aKt'n') 

The  partial  derivatives  in  brackets  are  reduced  to  a  total 
derivative  along  a  ©  =constant  curve  by: 


dlO  ■=■  Sm  6  COS.  S 

<=Hl  ^  ™ 

E)q.  (14)  becomes: 


(15) 


-*-sm*(2>tenS  4-  Wm  (2  cos* (3 

bn  v  1 

+  iinM  d(Rtan(^\^\u  ^/aWac^  dR\  (16) 
RcosS  dfi  i  ^  I  "cos  6  dl  ) 

=  2  Co^2p)  K(r^ 


Eq.  (16)  is  reduced  to  an  ordinary  differential  equation 
along  a  characteristic#  whose  tangents  at  any  point  satisfy 
the  relation: 

tanTT*  2>»m*  fc-Vam  5  (17) 


According  to  Ref.  8#  these  characteristics  are  regarded  as 
circular  projections  on  a  meridional  plane  of  unique  spa- 
cial  lines  lying  on  the  blade  surface.  Introducing  Eq. 
(17)  into  Eq.  (16): 

-+  \M^  V,  +  Ylm 60 V,  =  Y,  (18) 
d* 
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where  x  denotes  distance  along  a  characteristic,  andj 


V, 

y2 


Cos  H  [2  Cos2fe  +  Sin  (2^ 
L  RcosS 


d(R  "t3n 

d* 


(19a) 


2cos^  (19b) 

cos  S  du 


\s  -  2  cos 'y  VC (.tT) 


(19c) 


4.  The  value  of  K(n)  is  dependent  upon  conditions  ahead 
of  the  rotor,  as  shown  by  Eq.  (12).  The  normal  component 
of  the  equation  of  motion  for  steady,  isentropic  flows 
ahead  of  the  rotor  blading  is: 


V, 


vn 


-  v*  -  0 


(20) 


Eq.  (10)  shows  that  for  isentropic  conditions  at  the  en¬ 
trance,  the  total  enthalpy  and  the  product  RVu  are  constant 
along  a  meridional  streamline.  This  condition  is  proved 
from  the  u  and  m  components  of  Eq.  (4).  Thus: 

,2 


1^-"  +v;  x,  +  x2  =0 

6)t\ 


(21) 


X,  =  2km 


(22a) 


X 


_  Vu 


MRVuN  _  bW 
Tri 

Re  r^RVu^  _  ^VUTI  4H* 

R2 L dVii'  ci^e]  d n 


The  solution  of  Eq.  (21)  is: 


(22b) 


(23) 
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This  solution  is  solved  for  given  entrance  conditions 
and  the  volumetric  flow  rate,  Q,  by  iterating  the  hub 
velocity,  VmH,  until  the  flow  rate  is  satisfied  by: 


Q  =  2 tt  (  RVrn  tamp  dsn 

(24) 

M 

At  the  leading  edge: 

(25) 

Thus,  K(n)  is  solved  for  given  entrance  conditions. 

5.  The  remaining  terms  of  Eqs.  (19)  are  determined  from 
blade  geometry  and  the  orientation  of  the  meridional 
streamlines  within  the  rotor.  These  streamlines  must  be 
assumed  and  corrected  by  successive  approximations.  The 
location  of  the  streamlines  at  the  leading  edge  are  de¬ 
termined  from  entrance  conditions.  Eq.  (18)  is  solved 
along  a  characteristic  until  the  flow  rate  is  satisfied  by: 

Q •  2.-TC  r  R\NmC05-S-  dx  (26) 

*o 

B.  Blade  to  Blade  Analysis 

1.  It  is  assumed  that  the  relative  velocity  varies  line¬ 
arly  across  the  blade  channel.  This  assumption  is  taken 
from  thin  airfoil  theory,  where  blade  profiles  are  replaced 
by  bound  vortices.  The  relative  velocity  from  the  axisym- 
metric  solution  is  considered  to  be  the  mean  velocity 
along  the  periphery  at  any  point  in  the  meridional  plane. 
These  assumptions  lead  to: 

\Vp  *  Vs/-  AW 
W3  *  W  -v  AW 


(27a) 

(27b) 
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where  p  and  s  refer  to  the  pressure  and  suction  sides 
of  the  blade. 

2.  The  difference  between  relative  velocities  across  the 
blade  is  related  to  the  static  pressure  difference  byj 

^4-  f/2  Na/2  =  COm+emV  (28) 

Thus: 


A\a/  =  *£ 


(29) 


The  pressure  difference  across  a  blade  element,  which 
is  projected  onto  a  meridional  plane,  is  related  to  the 
moment  about  the  axis  of  rotation,  exerted  by  the  flow  on 
that  element,  by? 

=  -  g  (R  A-P  clndm)  (30) 


This  moment  is  also  derived  from  the  momentum  theorem  and 


Eq.  (9a): 

Atf  - 


--^Fr  Aedn  dm  P  +  (31) 

00  L  c)  V'o  J 


The  minus  sign  indicates  that  the  moment  oppose  the  rota¬ 
tion  of  the  rotor.  The  arc,  A©  /  is  expressed  in  terms 
of  an  equivalent  blade  thickness.  At/  ,  which  is  measured 
perpendicular  to  the  meridional  plane. 

-v/|  -+  C.0S  S 

Co*  (3  Y 


A©  - 


2JT  _  At! 

N  R 


(33) 


where:  N  =  number  of  blades 
At  =  blade  thickness 
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Equating  Eqs.  (30)  and  (31)  and  combining  with  Eqs.  (29) 
and  (33): 


Aw  CQ%<$>  (34) 

^  R  >  2  ^ - 


III.  NUMERICAL  SOLUTIONS 


A.  Design  or  Inverse  Solution 

1.  The  application  of  the  preceding  theory  to  impeller 
design  was  investigated  in  an  attempt  to  produce  a  physi¬ 
cal  model  for  the  direct  solution,  or  flow  analysis.  Over¬ 
all  dimensions  of  an  axial  entry,  mixed  flow  pump  impeller 
were  provided  by  Professor  Vavra,  which  included: 

a.  Hub  and  tip  profiles  of  meridional  contour 

b.  Leading  and  trailing  edge  contours  in  a  merid¬ 
ional  plane 

c.  Blade  entry  and  exit  angles 

The  problem  was  reduced  to  that  of  determining  blade  shapes, 
which  would  not  only  satisfy  these  boundary  conditions,  but 
would  conform  to  practical  structural  and  aerodynamic  limi¬ 
tations. 

2.  A  drawing  was  made  of  the  meridional  plane  and  a  series 
of  approximate  normals  were  constructed.  A  meridional 
streamline  was  constructed,  dividing  the  rotor  annulus  in 
half. 

A  'te  =  +  (  35a) 

Iterations  were  made  on  R_  along  a  normal  until: 

m 

A|/2»  S.  (R-r  +  Rh)  ^-r  (35b) 

The  slope  of  this  streamline,  \  ,  was  measured,  plotted, 

smoothed,  and  the  streamline  reconstructed.  Normals  were 
corrected. 

The  distribution  of  the  flow  angle,  (3  ,  along  the 

mean  streamline  was  assumed  in  the  form: 


14 


(2>  =  A  +  bW 


(36) 


A  and  B  were  evaluated  from  initial  conditions  at  the  lead¬ 
ing  and  trailing  edges.  Values  of  ©  were  calculated  by 
numerical  integrations  (Ref.  11).  Prom  Eq.  (1): 


(37) 


where : 


The  slope  of  the  ©  =constant  curves  at  any  point  on 
a  streamline  may  be  calculated  by  Eq.  (2),  if  £  is  known. 
The  choice  of  an  £  distribution  is  governed  byj 

a.  Maximum  £  is  limited  by  structural  considerations 

b.  £  distribution  must  be  compatible  with  fabrica¬ 
tion  tooling  practices.  An  example  is  given  in 
Ref.  9. 

The  design  problem  could  have  been  simplified  by  using 
radial  blades,  where  £  =  0.  However,  this  is  contrary 
to  the  intent  of  this  thesis  in  presenting  methods  applic¬ 
able  to  arbitrary  designs.  Therefore,  several  £  distri¬ 
butions  were  assumed,  and  the  angles,  S  ,  computed  on 
the  mean  streamline. 

3.  The  problem  of  analytically  generating  a  complete 
family  of  ©  =constant  curves  was  as  yet  unsolved.  The 
(3  and  £  distributions  on  streamlines  other  than  the  mean 
are  not  independent  of  the  mean  distributions.  Arbitrary 
choice  of  these  distributions  would,  in  all  probability. 
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generate  impractical  blade  shapes.  These  distributions 
should  be  assumed  as: 

£  -  6. 

It  appeared  that  the  design  problem  had  become  more  complex 
than  originally  anticipated.  It  was  decided  that  a  complete 
treatment  would  detract  from  the  main  objective  of  this 
thesis.  Therefore,  the  design  analysis  was  discontinued. 

3.  Meridional  Plane  Analysis 

1.  The  approach  used  to  solve  3q.  (18)  along  a  chosen  num¬ 
ber  of  characteristics  is,  briefly: 

a.  Assume  meridional  streamlines  in  the  vicinity  of 
the  leading  edge  which  divide  the  flow  into  ap¬ 
proximately  equal  increments. 

b.  Construct  normals,  thus  establishing  a  coordinate 
grid  system. 

c.  Measure  the  necessary  blade  physical  characteris¬ 
tics  at  each  grid  point. 

d.  Generate  a  characteristic  curve  within  the  grid 
system. 

e.  Compute  the  coefficients,  Y,  of  Eqs.  (19). 

f.  Solve  Eq.  (18)  for  at  the  intersections  of  the 

characteristic  with  the  streamlines,  iterating 

W  „  until  the  flow  rate  is  satisfied. 
mH 

g.  Correct  these  intersections  until  the  streamlines 
divide  the  flow  rate  into  the  prescribed  increments. 

h.  Recompute  W  at  the  new  intersections. 
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i.  Project  the  corrected  streamlines  farther  into 
the  rotor  channel ,  and  repeat  the  process  for  a 
new  characteristic. 

Thus,  the  meridional  streamlines  are  generated  from  lead¬ 
ing  to  trailing  edge,  and  the  distribution  of  the  relative 
velocity  on  a  meridional  plane  of  the  rotor  is  computed. 
Steps  a,  b,  c,  and  i  are  solved  by  the  computer.  The  re¬ 
maining  steps  are  carried  out  by  graphical  means. 

2.  There  are  two  categories  of  initial  data  required. 

The  first  is  a  complete  physical  description  of  the  rotor. 
Finally,  flow  conditions  ahead  of  the  rotor  must  be  pre¬ 
scribed.  Normally,  sufficient  data  are  available  from 
drawings  of  the  impeller  to  construct  a  meridional  plane. 
The  ©  =constant  curves  are  best  produced  by  orthogonal 
projections  from  a  three- view  drawing.  The  meridional 
plane  for  the  impeller  used  in  this  analysis  is  shown  in 
Fig.  1. 

The  flow  conditions  at  an  entrance  station  must  either 
be  prescribed  from  an  analysis  of  the  machine  installation, 
or  assumed.  The  specific  data  required  are: 

a.  Thermodynamic  data  for  the  fluid 

b.  Velocity  distribution 

c.  Flow  rate 

d.  Inlet  channel  contours 

e.  Impeller  RPM 

Conditions  for  this  analysis  are  simplified  by  assuming: 


17 


VL=  Vw,_=  co-nit. 

O 

e>^e 

This  flow  is  somewhat  impractical  since  it  does  not 
reflect  actual  velocity  and  energy  distributions  imposed 
by  intake  ducting  or  guide  vanes.  However ,  these  distri¬ 
butions  do  present  a  definite  off-design  condition  at  the 
leading  edge.  Design  RPM,  and  a  flow  rate  compatible  with 
actual  machine  operating  conditions ,  are  used.  Initial 
data  are  listed  in  Table  I.  The  meridional  velocity  at 
the  leading  edge  was  computed  byj 

Vnra.  =  (38) 

A^  was  computed  along  the  leading  edge  by  Eq.  (35b). 

3.  Streamlines  are  extended  from  initial  points  on  the 
leading  edge,  such  that  the  flow  rate  is  divided  into  equal 
increments.  Eight  divisions  were  used  for  this  analysis, 
as  shown  in  Fig.  4. 

A  grid  system  is  constructed  in  the  vicinity  of  the 
leading  edge.  A  point  on  the  hub  streamline  is  chosen  as 
the  starting  point  of  the  first  characteristic.  This  point 
is  determined  by  predicting  the  approximate  alignment  of 
the  characteristic.  Eq.  (17)  indicates  that  the  slope  of 
the  characteristic,  ,  is  proportional  to,  but  less  than, 
the  slope  of  the  ©  =constant  curve  which  originates  from 
the  starting  point.  H  also  has  the  same  sign  as  $  . 

These  guidelines  fix  the  characteristic  within  a  desired 
region.  These  considerations  also  assist  in  estimating 
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the  proper  grid  width.  A  "starting  normal"  is  constructed 
from  the  starting  point  on  the  hub  to  the  tip  streamline. 

A  second  normal  is  constructed  adjacent  to  the  first,  in 
the  predicted  direction  of  the  characteristic.  The  dis¬ 
tance  between  the  two  normals  on  each  streamline  is  used 
to  fix  the  remaining  grid  points.  These  points  do  not 
coincide  with  true  normals,  however  equal  grid  spacing  on 
streamlines  is  used  to  good  advantage  in  subsequent  calcu¬ 
lations. 

The  values  of  X  ,  R,  and  Z  are  measured  at  each  grid 
point.  The  values  of  ©  were  read  at  the  intersections  of 
the  ©  =constant  curves  and  the  streamlines,  plotted,  and 
©  at  the  grid  points  found  by  interpolation.  The  angles 
5  are  determined  in  a  similar  manner.  The  deviations  of 
the  ©  =constant  curves  from  the  radial  direction,  5-X  , 

are  measured.  \  is  added  to  the  interpolated  results  to 
obtain  8  . 

4.  Each  point  on  a  characteristic  must  satisfy  Eq.  (17), 
therefore,  the  angles  /3  and  5  must  be  calculated  for  any 
point  on  a  streamline.  ft  is  defined  by  Eq.  (1).  The 
derivative  of  ©  at  a  grid  point  on  the  j  streamline  is 
approximated  by  five-point  difference  formulas  from  Ref.  12. 
Forward  Differences 

=  -22  -4-  u-fi,  ©l-m  ~  3fo  ©L+a~*~ 1  k ©,-■>- 39^ ) 

dYni  12  ANj 
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Central  Differences 


-  (®L-2  -  ~  8  (e>L-> _ SLL±A.  (39b) 

d  12  A  Nj 

Backward  Differences 

-  3  <5>;  -4.  ~  \  (o  Q|-^.  -t-  3  4>  Q  L-  Z  —  4-8  ©£_->+ 2. 5iS>i(  3qc) 

cl™j  12.ANj 

j.*!* 

A  Nj  is  the  grid  spacing  on  the  j  streamline.  The  de¬ 
rivative  of  ©  at  points  other  than  grid  points,  are  cal¬ 
culated  by  interpolating  linearly  between  adjacent  normals. 

C\Q*  -  c\e»L  .  MQl-h  _  d£Lu\  tAnS 

6^  dro^  \  ciYn^  J  ANj  ^40' 

x  y  ^ 

dN  is  the  distance  from  the  i  grid  point  to  the  point  P  . 
A  consequence  of  Eqs.  (39)  and  (40)  is  that  the  minimum 
width  of  the  coordinate  lattice  is  six  grid  points. 

The  values  of  R,  \  ,  and  &  at  a  point  P  are  also 

calculated  by  linear  interpolation.  (3  and  'X  are  calcu¬ 
lated  by  Eqs.  (1)  and  (17). 

5.  A  characteristic  is  approximated  by  a  polygon  with  each 

side  terminated  by  adjacent  streamlines.  Each  intersection, 
* 

P  ,  is  located  from  information  calculated  at  the  preceed- 
ing  intersection.  In  this  derivation,  the  lines  between 
grid  points  are  assumed  to  be  straight.  The  derivation  is 
shown  graphically  in  Fig.  5. 

Assume  that  PQ  is  established,  and  the  angle  'X  , 
known.  A  straight  line  from  PQ  to  P'  on  the  next  stream¬ 
line  at  the  angle  .  The  distance,  dN'  iss 

dN'  =  T  (41) 
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where : 


ors 


dM  = 


dM  = 


dR 
CC£>  7s 


Sin  A 


(42a) 


(42b) 


Eq.  (42b)  is  used  for  A  greater  than  45  degrees.  The 
angle  '?f/  is  calculated  by  the  methods  of  paragraph  4. 
A  line  is  extended  from  PQ  to  P"  at  this  angle.  The 


point  P  is  established  by: 


d  M  -  d  Kl;  -V 


(43) 


P*  is  the  reference  point  for  the  extension  of  the  char¬ 
acteristic  to  the  next  streamline.  The  polygonal  approx¬ 
imation  of  the  completed  characteristic  is  smoothed  through 
the  points  P*.  (Fig.  4) 

There  are  a  number  of  considerations  that  must  be 
accounted  for  in  generalizing  this  method  for  computer  pro¬ 
gramming.  dN*  originates  at  the  starting  normal.  The 
sign  of  dN'  is  the  same  as  that  of  /"X'  .  The  sign  of  5N 
depends  on  the  difference  between  "S'  and  .  The  sign 

of  dN*  may  be  plus  or  minus,  and  its  magnitude  may  be  great¬ 
er  or  less  than  that  of  dN".  The  value  of  dN  in  Eq.  (40) 
for  a  positive  dN*  is: 

dN  -  dWK  .  . 

(44a) 

while,  for  negative  dN*: 

dN/=  (44b) 

6.  The  coefficients,  Y,  of  Eq.  (18)  contain  four  elements 
which  are,  as  yet,  undetermined. 

a. 
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b.  ^  (  R 

c.  c\R/d5 

d.  KCyTN 


The  curvature,  k,  of  a  given  curve  is  extremely  difficult 
to  calculate.  If  the  curve  can  be  expressed  analytically: 


=  - & 


(45) 


where  R'  and  R"  are  the  first  and  second  derivatives, 
respectively,  of  the  function  R  =  f^Z). 

An  attempt  was  made  to  determine  the  equations  of 
the  hub  and  tip  contours  in  cartesian  coordinates:  R  =  f(Z). 
Two  polynomial  approximations  of  these  curves  were  made 
using  CDC  Cooperative  Library  routines.  The  coordinates 
of  the  curves  at  each  quarter  inch  along  the  axis  were  in¬ 
troduced  as  data.  Only  high  order  polynomials  fit  these 
data  points  properly.  The  derivatives  of  these  equations 
were  difficult  to  obtain  without  frequent  interruptions  in 
computer  operation.  When  calculated,  the  second  deriva¬ 
tives  reflected  the  sinuous  nature  of  polynomial  approxi¬ 
mations,  resulting  in  inaccurate  curvatures  which  varied 
in  sign  along  a  given  curve.  These  methods  were  considered 
to  be  unacceptable  for  this  study. 

The  derivatives  of  the  hub  and  tip  contours  were  cal¬ 
culated  with  various  finite  difference  equations.  (Ref.  12) 
The  three  and  four  point  calculations  were  inaccurate.  The 
higher  order  solutions  resulted  in  erratic  second  deriva¬ 
tives.  It  became  apparent  that  all  methods  were  very  sen¬ 
sitive  to  the  accuracy  of  the  data.  It  was  found  that 
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inaccuracies  in  the  third  decimal  place  of  radius  data 
were  sufficient  to  introduce  gross  errors  in  the  second 
derivative. 

A  semi-graphical  approach  was  used  to  determine  the 
derivatives.  The  first  derivative  was  calculated  by 
finite  differences,  and  the  results  plotted  and  smoothed. 
This  process  was  repeated  to  calculate  the  second  deriva¬ 
tive,  using  the  plot  of  the  first  derivative  for  input. 
Data  taken  from  the  two  plots  were  used  to  calculate  the 
curvature  in  Eq.  (45).  Hub  and  tip  curvatures  are  plot¬ 
ted  in  Fig.  6. 

This  procedure  is  obviously  incompatible  with  con¬ 
tinuous  computer  operations.  The  computation  of  the 
streamline  curvatures,  k^,  at  points  P*,  would  be  quite 
tedious  using  this  method.  Therefore,  k^  is  approximated 
within  the  flow  passage  by  linearly  interpolating  between 
hub  and  tip,  using  the  data  of  Fig.  6.  k^  is  assumed  to 
vary  in  equal  increments  between  streamlines: 

kn>j  =  <46> 

Linear  interpolation,  similar  to  Eq.  (40),  is  used  to  cal¬ 
culate  k  between  grid  points  on  a  streamline, 
m 

The  following  derivation  defines  (R  tan{3  )  in 

terms  of  known  quantities.  From  Eq.  (15): 

dlStan@ )  „  3fnS  +  c<»  S 

4T  c>Y\ 

(47) 

-  R  (smS  $  4  S  4"^*  p 
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From  Fig.  2: 


dR 

dir 


Cos  —  o, 


(48^ 


From  Eq .  ( 1 ) : 


c)(Rt.sn^  __  R  ?4©  4®.  . 

5  YV>  c> To2  2>Yn  c)m 


a(R4yo  <y>  =  p  d^Q  de  c)R 

d  v\  r'  £  ml)Y\  dvo  ciN) 


d  (Rtan  o’)  iR2(irtl5  cP©  +  C0Sl5  i£a.  \  +  2&,t*n& 

1  dwdh  )  ^ 

Grouping  terms  of  Eq.  (19a) : 


_-!— 5  d  (R tan  _  q  -  fj/tanS  d®,  +  \ 

RCcr.5  dQ  ^2-  \  dYn2  dYftdn ) 


(49) 


■4  8  D>  ~bs>r> 

R  C-OS^ 

2  2 

The  term  d  ©/dm  at  a  point  P*  is  computed  by  as¬ 
suming  a  parabolic  distribution  of  do /dm  between  grid 
points. 

dS> 

dw 


=  a  4-  brvi  -v  CVTi 

(50) 

~  d&  -4  P*  —  "b$n  ^ 

'  aw  ot  -  r 

“  dT^x  at  K 

following  the  derivation  for 

(51) 

b. 

bYn(  +CVr\})  TO* 

(52a) 

4  brn2  -4-  cyo|^ 

(52b) 
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(53) 


Subtracting  3q.  (52b)  from  (52a): 

b=  W  [(l&V  (|™Q 

"try, TO 2  -  Yr^mf- 

m^  and  m2  are  defined  in  Fig.  8. 

j  2 

The  function,  ,  is  calculated  by  dif ferentiat- 


j  r\ 

ing  along  the  normal  nearest  the  point  P*.  A  para¬ 

de 


bolic  distribution  of  over  three  streamlines  along  the 
normal  is  assumed.  The  coefficient,  b,  is  derived  in  the 
same  manner  as  that  described  in  the  preceding  paragraph. 
In  this  derivation  the  two  succeeding  grid  points  from  the 
origin  are  used,  instead  of  adjacent  points.  Calculations 
for  points  on  the  last  two  streamlines  are  made  by  using 
the  two  preceding  grid  points. 

Forward  Differentiation 


(54a) 


b  = 


-  rC,  [(f&u- 

-  n.j’n.2; 


(54b) 


The  increments  n^  are  functions  of  the  grid  spacing,  dM. 

The  values  of  K(n)  are  required  for  the  calculation 
of  Y^.  Initial  conditions  have  simplified  Eq.  (12).  With 
Eq.  (25): 


k(ri> 


__  _  (  RuVm  dt_ 

“  at-  dm 


oz  =  3J  • 

2>0 

dL  _  dL  d  R  _  COS  X 
d  y>  dR  'dm  5m  c* 
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55) 

(56) 


(57) 


where  cX,  is  the  slope  of  the  leading  edge.  is  cal¬ 

culated,  using  the  forward  difference  formula  of  Eq.  (54a) 
in  conjunction  with  Eq.  (1).  The  terms  in  parentheses  in 
Eq.  (55)  are  grouped,  and  the  derivative  calculated  using 


the  parabolic  distribution  methods  of  Eqs.  (54). 
F(LX  *  RLVm  oo  F?* 

In  summary,  the  coefficients,  Y,  are  reduced  to: 

V,  =  cos2  (2,  -v  siv>  (2  (3^  •  D^j 


Yo  =  2  Cos 


s'\v>  (2 
Cos  5 


D, 


Y-i=E.cos'Y  Cos*  (2>  K(rT) 


(53) 

(59a) 

(59b) 

(  59c) 


7.  The  equation  of  motion,  Eq.  (18),  is  solved  by  succes¬ 
sive  approximations,  until  the  distribution  of  Wm  along 
the  characteristic  satisfies  the  flow  rate  of  Eq.  (26). 

The  steps  in  this  solution  are: 

a.  Assume  WmR  2 

Solve  Eq.  (18)  for  ^  ^  ^  h 


b. 

c.  Let: 


d  X 


d  X?  = 


<dx 


= 

^  ~  cos'Yh 


(60) 

(61) 


d.  Solve  Eq.  (18)  for 

e.  Let: 


cl  V/ 


m2 


cl  X 


=  \Nwh  -V-  -i  (  ■+  (62) 

f.  Repeat  for  Wm3,  using  Wm2  in  (a). 

The  variables  under  the  integral  sign  in  Eq.  (26)  are 
grouped. 
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(63) 


EM  -  r  co^n( 


Q  -  2TT  (  h  M  cl* 


Numerical  integrations  are  performed  over  two  adjacent 
streamlines  at  a  time,  and  the  results  are  successively 
summed.  A  parabolic  distribution  is  assumed  for  E(x). 

EM  =  a*  b*  -V  C  X*  (64) 

The  coefficients  are  derived  in  the  same  manner  as  that 
used  for  the  development  for  the  difference  equations. 

a-  E0 


where : 


X.X^  -  X-Z.xf 

c  =  ME.-E„VX,(Ei.-E^ 

x?xz-x^  x, 

Xj®  dX> 

“dX,4-  d  Xz. 

The  integral  of  E(x)  is: 

j  Ewdx  =  a(xte42-  |(x^+l-x^ 

+  3  (XX+2  ~  XV) 

where :  ^ 

Xt.  =  Z  dx-L 

K  L  =  2 

In  case  there  exists  an  odd  number  of  strearatubes,  the 
final  integral  is  calculated  by: 

E  M  «  a  -v  bx 
d  »  Lo 


(65a) 

(65b) 

(65c) 


(66) 


(67) 


(68) 

(69a) 
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(69b) 


b  = 


E»-H.o 

dX, 


J  ew  ax  -a(nh+r  +  k  (x^+l  -  xV) 


(70) 


8.  The  intersection  of  the  characteristic  with  each  stream¬ 
line,  excluding  hub  and  tip,  is  connected  until  adjacent 
streamlines  divide  the  flow  rate  into  M-l  equal  increments. 
Calculations  progress  from  streamline  M-l  to  streamline  2. 
Each  calculation  iterates  the  increment  Sx  until: 

VO  1 


Q (w) -  2TT 


ECA'VSx')  (dx+cSxV 


VC’ 

The  variables  in  E(x)  are  corrected  by  Sx. 

Rk  ■»  ■+  (ar^t 

-  Rl  "V  5  X  COS  (Yi.  — 

vw*  w*** 

/  aw^\ 

~  +  2>N*hA  dx  )l 


(71) 


(72) 


(73) 


Y 


A  »  C»w ,  -  TTt)  to 


(74) 


L  ) 


K  indicates  the  new  point  of  intersection,  and  i,  the 
original  point.  Eq.  (71)  is  solved  by  introducing  the 
proper  limits  into  Eqs.  (66)  or  (70). 

9.  Preliminary  calculations  for  the  blade  to  blade  solu¬ 
tion  are  included  in  the  computer  program  for  the  merid¬ 
ional  plane  analysis.  Variables  of  Eq.  (34)  are  grouped. 
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and  the  groupings  computed  at  the  corrected  intersections 
of  the  characteristic  and  streamlines. 

=  (  23T-  ^  (*  (75a) 

+  R  (75b) 

C.  Computer  Programs?  ROTOR  1  and  LEDGE  1 

1.  The  preceding  methods  are  translated  into  FORTRAN  com¬ 
puter  language  compatible  with  the  CDC  1604  digital  computer. 
Program  ROTOR  1  calculates  a  complete  solution  for  one 
characteristic,  from  the  location  of  the  characteristic 
to  the  correction  of  assumed  streamlines  and  the  relative 
velocity  profile.  The  solution  procedes  in  an  orderly 
fashion,  much  like  the  development  of  the  preceding  section. 
Subroutines  are  used  for  repetitive  calculations.  Control 
of  decisions,  iterations,  and  progressive  development  is 
maintained  in  the  main  program.  Subroutine  PSTAR  locates 
the  characteristic,  accounting  for  all  signs  of  dN*.  Sub- 
routine  DDM  calculates  the  first  derivative  of  ©  ,  , 

by  finite  differences.  Subroutine  ANGLE  calculates  (2>  and 
'X  .  Subroutine  COY  calculates  the  elements  of  the  coef¬ 
ficients,  Y,  except  K(n).  The  coefficients,  Y,  are  calcu¬ 
lated  in  the  main  program.  The  relative  velocities  are 
calculated  in  subroutine  RELVEL,  and  the  flow  rate  in  sub¬ 
routine  FLOW.  Control  of  the  iterations  is  maintained  in 
the  main  program.  Iterations  for  streamline  corrections 
are  made  in  the  main  program,  using  subroutine  FLOW. 
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2.  Entrance  conditions  are  calculated  in  a  separate  pro¬ 
gram,  since  these  computations  are  not  required  after  the 
locations  of  streamlines  on  the  leading  edge  are  deter¬ 
mined.  Program  LEDGE  1  calculates  K(n)  as  defined  by  Eqs. 
(55),  (56),  and  (57),  with  the  exception  that  cos  A  is 
omitted.  This  variable  is  a  function  of  position  on  a 
streamline,  and  therefore  is  included  in  the  calculations 
for  Y^.  The  results  of  LEDGE  1  are  used  as  inputs  to 
ROTOR  1. 


It  should  be  noted  that  the  first  term  of  K(n)  in  Eq. 
(12)  is  omitted,  since  0  in  this  analysis.  This  term 

is  introduced  as  input  (DHEDN)  into  ROTOR  1.  It  is  intend¬ 
ed  that  this  term  would  be  calculated  separately  from  given 
thermodynamic  data  at  the  entrance. 

3.  Definitions,  flow  diagrams,  and  program  listings  are 
included  in  the  Appendix.  Many  control,  indexing,  and 
grouping  names  are  undefined.  Their  usage  may  be  inter¬ 
preted  from  the  developments  of  the  preceding  section  or 
from  the  program  listings. 


D.  Results  of  ROTOR  1 

This  program  computed  properly,  but  the  results  were 
unacceptable.  Streamline  corrections  were  large  (up  to 
1/2  inch),  imposing  improbable  velocity  distributions. 

The  characteristics  appeared  to  be  properly  located,  how¬ 
ever,  the  flow  angles,  /3  ,  were  somewhat  arbitrarily 
distributed  along  normals  and  characteristics. 


30 


A  number  of  test  print-outs  were  made  in  order  to 
locate  the  inaccuracies.  (The  listing  of  Rotor  1  in  the 
Appendix  includes  test  print-out  instructions. )  It  was 
discovered  that  the  derivatives  of  ©  ,  calculated  by  Eqs. 
(39),  (40),  (53),  and  (54)  were  inaccurate.  The  second 
derivatives  were  particularly  incorrect.  The  methods  used 
are  considered  to  be  sound,  however,  the  degree  of  accuracy 
primarily  depends  on  the  accuracy  of  the  input  data.  In 
this  analysis,  the  0  =constant  lines  were  originally  plot¬ 
ted  from  approximations.  Errors  introduced  in  plotting 
were  compounded  by  graphical  interpolations  (and  extrapola¬ 
tions)  for  ©  values  at  grid  points.  The  erroneous  ©  de¬ 
rivatives  were  introduced  in  the  computations  for  the  co¬ 
efficients,  Y,  resulting  in  an  improbable  flow  solution. 

E.  Alternate  Meridional  Plane  Analysis 

1.  It  was  decided  to  eliminate  the  derivatives  of  ©  from 
the  calculations,  therefore,  it  was  necessary  to  introduce 
known  values  of  (2>  ,  It  is  reasonable  to  assume  that  these 
data  would  be  available,  or  could  be  computed,  from  detailed 
drawings  of  the  impeller.  For  this  analysis,  /2>  distribu¬ 
tions  along  the  hub  and  tip  contours  were  obtained  from  data 
used  in  the  original  design  work.  However,  cannot  be 
specified  on  the  internal  streamlines  until  the  streamline 
is  constructed.  A  distribution  of  from  hub  to  tip  must 
be  assumed. 

2.  (h  along  the  hub  and  tip  contours  is  shown  in  Fig.  9. 
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(76) 


It  is  assumed  that  (3>  varies  linearly  along  normals  be¬ 
tween  these  contours. 

The  function  D2  of  Eq.  (49)  becomes? 

^ —  i  d-i  F?  ±&r\  [£) 


a 


RCo^<5  ai 

—  g  ~crhay>  (^>  -^"crLan  (b  t>t  tav>  ft> 
dm  dvT  RcoS>3 


(77) 


The  distribution  of  tan  (3>  between  two  grid  points  on  a 
streamline  or  normal  is  assumed  to  be  linear. 


At  the  ith  grid  points 


^)~LaY\  p>  ~  ~h&Y> 

ANj 

^“Ldr*  -  ~tdr(  fij-vv 

b  ric  dM'i 


3 L  -  tav>  f^i-vv  —  fel-\ 

At  a  point,  P*s 


b  V) 


V  2)  L-n 


(78a) 


(78b) 


Z  ANJ^ 

(79) 

2)-tan  (h\  ^ N 

(30; 

2>  n l  ) 

The  angle  S  is  still  measured  on  the  0  =constant  curves. 

Leading  edge  calculations  are  also  modified.  A  linear 
distribution  of  along  the  leading  edge  is  assumed  in  the 
same  manner  as  Eq.  (76). 


F.  Computer  Programs  ?  ROTOR  2  and  LEDGE  2 

1.  Frogram  ROTOR  2  is  a  modification  of  ROTOR  1,  with  the 

new  calculations  for  fi  replacing  the  ©  derivatives.  The 
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order  of  some  solutions  has  been  changed  to  simplify  the 
mechanics  of  the  program.  Subroutines  DTD  and  ANGLE  are 
eliminated.  The  values  of  the  angles,  ^  ,  are  computed 

in  the  main  program  or  in  subroutine  PSTAR.  Subroutine 
PSTAR  is  modified  to  reflect  these  changes.  Subroutine 
DDL  replaces  COY,  and  calculates  the  derivatives  along  the 
©=constant  curves.  Curvatures  are  calculated  in  PSTAR. 
Subroutines  RELVEL  and  PLOW  are  unchanged.  Program  LEDGE 
2  is  a  simplification  of  LEDGE  1,  reflecting  the  introduc¬ 
tion  of  /3  data. 

Both  programs  are  diagrammed  and  listed  in  the  Appen¬ 
dix.  Only  new  or  modified  variable  names  are  defined. 

G.  Results  of  ROTOR  2 

1.  The  complete  meridional  plane  analysis  was  solved  by 
ROTOR  2.  Nine  characteristics  were  generated.  Fig.  4 
is  a  tracing  of  the  actual  construction  for  the  first  two 
characteristics.  The  grid  networks  for  both  solutions  are 
indicated.  Two  modifications  are  made  to  the  first  set  of 
estimated  streamlines.  The  first  modification  reflects  the 
changes  effected  by  the  calculations  for  characteristic  Cl. 
Each  change  required  recalculations  in  LEDGE  2.  The  re¬ 
sults  of  the  last  calculation  were  held  constant  throughout 
the  remaining  solutions. 

Data  for  each  characteristic  are  listed  in  Tables  II-l 
through  II-9.  The  "starting  normal"  originates  from  the 
"starting  point"  for  each  calculation.  The  "starting  point 
number"  indicates  the  position  of  the  normal  in  the  grid 
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system  for  each  solution.  The  intersection  of  a  character¬ 
istic  with  each  streamline  is  located  by  measuring  the  dis¬ 
tance  DNSTAR  (dN*)  from  the  intersection  of  the  starting 
normal  with  the  streamline.  Streamlines  are  corrected  by 
measuring  the  distance  DELTA  ( 5  x)  along  the  characteris¬ 
tic  from  the  original  intersection,  P*.  Relative  velocities 
are  listed  for  the  final  intersection  point.  Additional 
data  are  listed  as  checks  on  the  calculations. 

For  example,  the  calculations  for  characteristic  C8 
showed  a  discrepancy  between  the  correction  length,  &  x, 
and  the  computed  radius,  in  locating  the  new  streamlines. 
Test  print-outs  revealed  errors  in  the  calculated  grid 
lengths,  dM.  These  errors  had  accumulated,  resulting  in 
somewhat  radical  streamline  deviations  and  questionable 
velocities.  dM  had  been  calculated  by  Eq.  (42a),  which, 
at  the  time,  was  used  for  angles,  \  ,  of  60  degrees  or 
less.  This  limit  was  changed  to  45  degrees  and  Eq.  (42b) 
applied.  The  results  were  satisfactory. 

2.  The  completed  system  of  computer  characteristics  and 
streamlines  is  shown  in  Fig.  10.  True  normals  are  con¬ 
structed  perpendicular  to  these  streamlines.  Streamlines 
have  been  faired  and  extended  to  the  trailing  edge.  Data 
points  are  indicated. 

The  distributions  of  relative  velocity  along  odd  num¬ 
bered  streamlines  are  plotted  in  Fig.  11.  Values  at  the 
leading  edge  were  hand  calculated.  Velocities  at  the  trail¬ 
ing  edge  were  extrapolated  on  a  large  scale  plot.  Curves 
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are  faired  through  data  points  without  smoothing.  The  rela¬ 
tive  velocity  profiles  on  the  odd  numbered  normals  of  Fig. 

10  are  cross  plotted  in  Fig.  12. 

3.  Computer  run  time  for  one  ROTOR  2  solution  was  approxi¬ 
mately  two  minutes.  Construction  and  data  preparation  time 
was  about  four  hours.  Total  preparation  and  computer  time 
for  LEDGE  is  approximately  two  hours. 

H.  Blade  to  Blade  Analysis 

I.  Most  of  the  preliminary  calculations  have  been  accomp¬ 
lished  in  the  hub  to  tip  solution.  The  values  of  W,  Dwcoef/ 
and  DWfUnc  are  plotted  from  the  results  in  Tables  II ,  and 
the  curves  smoothed.  Data  are  read  from  these  curves  at 
equal  increments  along  each  streamline.  Extrapolated  data 
are  used  at  the  leading  and  trailing  edges.  An  increment 
of  1/2  inch  is  used  in  this  solution. 

The  derivatives  of  DWfunc  are  computed  by  the  five- 
point  difference  formulas  of  Eqs.  (39).  The  relative 
velocity  difference,  AW,  is  computed  by  a  form  of  Eq.  (34): 

AVJ  =  (81) 

2  dvr* 

The  relative  velocity  distributions  on  the  pressure  and 
suction  blade  surfaces  (driving  and  trailing  surfaces  in  a 
compressor)  are  calculated  by  Eqs.  (27). 

2.  Program  BLADE  performs  these  calculations  along  each 
streamline,  from  leading  edge  to  trailing  edge.  Computer 
run  time  was  approximately  30  seconds.  Preparation  time  was 
about  six  hours. 
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3.  Results  are  compiled  in  Tables  III-l  through  III-9. 
The  relative  velocity  distributions  along  the  hub,  tip, 
and  mean  streamlines  are  plotted  in  Figs.  13,  14,  and  15. 
The  linear  velocity  profiles  across  the  blade  channel,  at 
various  stations  along  the  mean  streamline,  are  shown  in 
Fig.  16. 
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IV.  DISCUSSION 


A.  Theoretical  Assumptions 

1.  The  assumption  of  axial  symmetry  is  in  keeping  with 
sound  engineering  practice.  Although  changes  in  the  cir¬ 
cumferential  direction  are  ignored,  incremental  pressure 
changes  across  the  blades  are  considered  by  introducing 
the  blade  force,  Pfi.  This  assumption  is  also  made  in 
Refs.  1  and  2.  The  general  theory,  postulated  by  Wu  in 
Ref.  4,  does  not  admit  to  axial  symmetry.  Wu  modifies  the 
development  of  Ref.  2  to  allow  for  deviations  from  axial 
symmetry,  accounting  for  these  deviations  with  a  thickness 
factor,  B,  This  factor  is  conveniently  used  as  an  inte¬ 
grating  parameter  in  defining  a  flow  function.  B  is 
interpreted  as  being  proportional  to  the  thickness  of  a 
stream  sheet  containing  an  arbitrary  stream  surface,  and 
to  blade  thickness  distribution. 

This  theory  is  sound,  and  would  produce  accurate  re¬ 
sults  if  explicitly  applied.  However,  the  factor  B  must 
be  assumed,  either  arbitrarily  or  from  available  data.  An 
example  of  an  assumed  distribution  of  B  is  found  in  Ref.  6. 
A  blade  to  blade  analysis  is  made,  using  the  distribution 
as  input.  This  report  states  that  such  data  are  obtained 
from  a  meridional  plane  analysis,  presumably  by  the  methods 
of  Ref.  4.  Yet,  Ref.  4  states  that  B  is  best  calculated 
from  the  results  of  a  blade  to  blade  analysis.  Wu,  in  Ref. 
13,  uses  a  value  of  B  equal  to  one,  which  imposes  axial 
symmetry.  He  states  that  methods  for  estimating  B  were 
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unavailable  at  the  time.  This  reference  also  shows  that 
the  changes  effected  by  this  assumption  are  small.  The 
advantage  of  using  approximate  values  of  B,  over  the 
axi symmetric  approach,  is  doubtful. 

2.  The  assumption  of  a  steady,  inviscid,  adiabatic  flow 
is  generally  accepted.  Exclusion  of  any  one  of  these  fac¬ 
tors  introduces  highly  complex  theoretical  considerations, 
which  are  difficult,  if  not  impossible  to  apply.  The 
assumption  of  isentropic  flow  within  the  rotor  is  reason¬ 
able,  since  known  exceptions  are  considered.  The  specific 
exceptions,  included  in  this  theory,  are  the  entropy  changes 
caused  by  discontinuities  of  the  leading  and  trailing  edges. 
The  analysis  of  this  thesis  is  based  upon  the  premise  that 
these  entropy  changes  are  small  and  can  be  ignored.  This 

is  also  implied  in  Ref.  1,  which  analyzes  flows  which  con¬ 
form  to  the  prescribed  blade  angles  at  the  leading  and  trail¬ 
ing  edges.  The  analysis  of  this  thesis  is  not  limited  by 
this  implication.  Entrance  conditions  are  not  restricted 
to  design  conditions,  and  Eq.  (25)  is  generally  applicable. 

3.  An  incompressible  solution  was  developed  in  this  thesis 
in  order  to  simplify  the  presentation.  Compressible  solu¬ 
tions  can  be  solved  by  this  method  if  the  volumetric  flow 
rate  is  replaced  by  the  mass  flow  rate.  Density  changes 
along  a  characteristic  can  be  computed  by  thermodynamic 
relations  which  equate  the  equilibrium  conditions  at  any 
point,  P*,  to  known  or  computed  entrance  conditions. 

4.  A  simplified  blade  to  blade  solution  is  used  in  this 


38 


thesis.  The  assumption  of  a  linear  velocity  distribution 
is  not  unique.  Ref.  1  assumes  a  linear  pressure  variation 
across  the  blade  channel.  Ref.  6  validates  this  assump¬ 
tion  for  design  conditions.  The  results  of  Ref.  3  show 
near-linear  velocity  distributions  for  both  the  compressible 
and  incompressible  solutions. 

B.  The  Design  Analysis 

1.  The  design  analysis  was  discontinued  because  of  time 
considerations.  However,  it  was  shown  that  the  theory 
could  be  applied  to  the  design  solution  as  well  as  the 
direct  solution.  Difficulty  arose  in  prescribing  practi¬ 
cal  blade  characteristics.  One  method  would  be  to  assume 
various  (3  and  £  distributions  as  functions  of  stream¬ 
line  and  normal  directions;  then  generate  a  blade  shape 
for  each  combination,  and  select  the  most  practical  solu¬ 
tion.  Selection  would  be  contingent  upon  the  results  of 
a  direct  analysis,  similar  to  the  method  developed  in  this 
thesis.  This  approach  appears  to  be  quite  long. 

The  question  is  whether  or  not  this  development  is 
warranted.  A  design  method,  based  on  the  theoretical  de¬ 
velopment  of  Ref.  1,  is  presented  in  Ref.  5.  The  hub  pro¬ 
file  is  calculated  from  successive  flow  solutions.  This 
approach  is  modified  by  Ref.  7  to  provide  for  arbitrary 
entrance  and  discharge  conditions.  In  effect,  the  initial 
conditions  establish  the  ©  =constant  curves.  The  solution 
provided  the  upper  limit  of  these  curves.  The  design 
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problem  is  practically  solved  by  the  initial  conditions. 


All  that  remains  is  a  system  of  repetitive  calculations , 
which  produces  an  acceptable  tip  profile.  In  comparison, 
the  method  proposed  in  this  thesis  would  generate  blade 
shapes,  rather  than  tip  contours.  It  seems  more  logical, 
and  perhaps  more  expeditious,  to  establish  feasible  physi¬ 
cal  contours  and  analyze  the  flow  under  prescribed  condi¬ 
tions.  The  contours  could  then  be  altered  to  eliminate 
undesirable  flow  phenomena,  and  to  produce  the  desired 
performance.  This  is  precisely  what  is  done  in  the  direct 
approach.  Therefore,  the  design,  or  indirect,  method  is 
considered  to  be  unnecessary. 

C.  Methods  of  Solution 

1.  The  developments  of  Refs.  8  and  9  present  a  relatively 
simple  approach  to  the  direct  problem  when  compared  to  Ref. 
4.  The  establishment  of  characteristic  solutions  provides 
a  method,  compatible  with  high  speed  machine  computations. 

A  limited  number  of  applications  of  Ref.  4  have  been  pre¬ 
sented.  None  have  attempted  a  solution  of  the  complete 
theory,  but  have  reduced  the  problem  to  solutions  for 
specific  applications.  Ref.  6  deals  with  the  blade  to 
blade  analysis.  Refs.  5  and  7  are  concerned  with  the  d*- 
sign  problem,  and  consider  the  methods  proposed  in  Ref.  4 
to  be  impractical  for  engineering  applications.  In  Ref. 

13  the  methods  of  Ref.  4  are  simplified  by  assuming  axial 
symmetry,  and  an  incompressible  solution  is  used  as  a  basis 
for  the  compressible  solution.  It  is  concluded,  therefore. 
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that  the  approach  presented  in  this  thesis  is  theoretically 
sound. 

2.  The  accuracy  of  the  solution  of  any  theoretical  analy¬ 
sis  depends  upon  the  precision  of  the  numerical  methods  used 
in  computations,  and  in  the  accuracy  of  the  input  data.  It 
is  difficult  to  assess  the  inaccuracies  introduced  by  the 
numerical  approximations  in  this  method.  It  would  be  im¬ 
possible  to  isolate  each  error,  even  if  test  results  were 
available  for  comparison.  Gross  inaccuracies,  such  as  those 
revealed  by  ROTOR  1,  can  be  detected.  Small  errors  may  can¬ 
cel,  or  large  errors  may  ensue  from  an  accumulation  of  small 
inaccuracies. 

For  example,  the  assumptions  of  linear  distributions 
of  curvature  and  flow  angles  are  questionable  in  the  vicinity 
of  the  axial  midpoint  of  the  test  impeller.  Hub  and  tip  cur¬ 
vatures  (Fig.  6)  follow  different  patterns  in  this  region. 

The  odd  distribution  of  hub  curvature  is  difficult  to  ex¬ 
plain.  The  difference  between  hub  and  tip  blade  angles  is 
a  maximum  in  this  region.  (Fig.  9)  These  two  quantities 
are  multiplied  in  the  calculations  of  the  coefficient,  Y. 
Difficulties  were  encountered  in  obtaining  accurate  results 
along  characteristics  in  this  region,  indicating  that  the 
errors  introduced  by  these  assumptions  may  have  multiplied. 
Such  reasoning  is  not  conclusive,  however,  since  the  results 
of  the  complete  solution  do  not  indicate  excessive  errors. 

Many  interpolations  were  used  in  this  solution.  The 
accuracy  of  each  depends  on  the  grid  spacing.  The  spacing 
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of  normals  was  limited  to  a  maximum  of  .3  inches.  Nine 
streamlines  are  considered  sufficient.  A  finer  network 
would  certainly  increase  the  accuracy  of  each  calculation , 
however,  the  process  of  measuring  data  would  become  tedious. 

The  reliability  of  all  calculations  also  depends  upon 
the  accuracy  of  the  measured  data.  The  meridional  plane 
was  drawn  to  double  scale.  Lengths  were  measured  to  an  ac¬ 
curacy  of  about  .05  inches.  The  angles,  \  ,  were  measured 
with  a  drafting  machine  to  an  accuracy  of  5  minutes.  (.1 
degree)  This  accuracy  is  reduced  to  about  .3  degrees,  since 
measurements  were  made  on  curved  lines.  Interpolations 
were  made  for  &  to  the  nearest  tenth  of  a  degree,  but  an 
accuracy  of  1/2  degree  is  more  reasonable. 

3.  The  general  approach  to  the  meridional  plane  solution 
can  be  applied  to  obtain  the  desired  degree  of  accuracy, 
within  the  limits  outlined  in  the  preceding  discussion.  A 
finer  grid  will  increase  accuracy  for  a  particular  charac¬ 
teristic  solution.  A  solution  can  be  repeated,  using  the 
previous  results  to  establish  new  initial  conditions.  The 
number  of  characteristic  solutions  used  is  arbitrary. 
Therefore,  a  detailed  investigation  can  be  conducted  in  dis¬ 
crete  regions  in  a  rotor  channel,  or  a  complete  solution  can 
be  rapidly  calculated,  using  a  few  characteristics. 

Theoretically,  the  method  can  be  used  to  calculate 
characteristic  solutions  outside  the  rotor  channel.  In  this 
case  the  defining  angles  of  the  blade  are  set  equal  to  zero, 
and  the  characteristics  coincide  with  normals.  In  stationary 
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cascades ,  the  rotational  effects  are  eliminated.  These 
solutions  are  excluded  in  the  computer  solutions  of  this 
thesis ,  because  of  the  schemes  adapted  in  performing  nu¬ 
merical  interpolations.  Solutions  can  be  calculated  in 
close  proximity  to  the  leading  and  trailing  edges  if  fic¬ 
titious  blade  characteristics  are  used  outside  these 
boundaries. 

D.  Results 

1.  The  results  of  the  first  method  of  solution  (ROTOR  1) 
were  unacceptable.  This  was  partially  due  to  the  numerical 
methods  used  in  calculating  the  ®  -derivatives.  However, 
inaccuracies  of  the  0 =constant  curves  are  considered  to 

be  the  primary  fault.  It  is  believed  that  accurate  ©  =con- 
stant  curves  can  be  constructed  from  three-view  drawings. 
Under  this  assumption,  ROTOR  1  can  be  utilized  with  accep¬ 
table  accuracy. 

2.  The  results  of  ROTOR  2  are  considered  valid.  Although 
test  results  are  not  available  for  comparison,  the  logic 
of  the  preceding  discussion  is  sufficient  proof,  in  that: 

a.  The  theory  is  sound. 

b.  The  assumptions  are  justified  for  engineering 
applications. 

c.  Numerical  methods  are  designed  to  minimize  inac¬ 
curacies  in  each  calculation,  so  that  errors  will 
not  accumulate. 

d.  Data  are  accurately  measured. 

3.  The  computed  system  of  streamlines  is  regular  and  follows 
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a  logical  development  through  the  rotor.  Fig.  10  is  a 
smooth  plot  of  this  system.  Deviations  from  data  points 
are  small.  Maximum  deviations  occur  in  the  region  of  rapid 
curvature  and  blade  angle  changes. 

4.  The  acceptable  system  of  streamlines  demonstrates  the 
consistency  of  the  solution.  The  accuracy  of  the  method  is 
measured  by  the  resultant  relative  velocity  distribution. 

This  distribution  completely  defines  the  flow.  Both  pres¬ 
sure  distribution  and  power  delivered,  or  required,  are 
calculated  directly  from  the  velocities.  Variations  in 
relative  velocity  indicate  areas  of  possible  flow  separation. 

The  accuracy  of  the  relative  velocities  at  the  lead¬ 
ing  edge  are  computed  directly  from  entrance  conditions, 
and  are  accurate.  Variations  of  velocities  along  selected 
streamlines  are  shown  in  Fig.  11.  Proper  trends  are  indi¬ 
cated  for  a  compressor  with  deleaded  type  blades.  Most  of 
the  work  input  is  accomplished  in  the  forward  section  of  the 
rotor,  diminishing  towards  the  trailing  edge.  It  should  be 
noted  that  the  energy  level  of  the  fluid  is  not  affected 
by  velocity  alone.  Eq.  (6),  together  with  Eq.  (9b),  shows 
that  the  local  enthalpy  is  also  a  function  of  peripheral 
speed,  go  R.  Therefore,  the  deloaded  blade  effectively 
reduces  the  deceleration  of  the  relative  velocity  at  high 
radii,  so  that  the  work  input,  thus  blade  loading,  will  not 
become  excessive.  The  irregular  nature  of  the  velocity 
distributions  in  Fig.  11  might  be  attributed  to  inaccura¬ 
cies  in  the  calculations;  however,  similar  presentations 
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in  Refs.  5,  6,  and  13  support  the  existence  of  these  irreg¬ 
ularities.  In  contrast,  the  velocity  profiles  along  nor¬ 
mals  in  Fig.  12  are  quite  regular. 

5.  The  applicability  of  these  results  can  be  illustrated 
by  an  investigation  of  flow  separation  along  a  contour. 
Relative  velocity  changes  along  the  hub,  near  the  leading 
edge,  appear  to  be  conducive  to  separation.  The  rapid  de¬ 
crease  in  velocity  implies  an  adverse  pressure  gradient. 

A  separation  parameter,  K,  is  derived  in  Ref.  14.  A  value 
of  K  greater  than  .045  indicates  probable  separation. 


.  1ST 


f  YJ* 

Reynold's  Number  is  defined  as: 

Re~  ^ 
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where:  1  =  length  along  blade  from  stagnation  point 

(3  =  average  blade  angle 
Along  a  streamline,  from  Eqs.  (6)  and  (9b): 
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For  incompressible  flow: 
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At  a  point  on  the  hub,  ra  =  2  inches: 
W  =  106  ft. /sec. 

H  =  -216  ft. /sec. -ft. 
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165 


=  sin  X  = 
con 


R  =  .5  ft. 

V  =  .00016  ft. 2/sec.  for  air 
( 3  =  49°  from  Fig.  9 


Therefore:  K  =  .0138 

Thus,  flow  at  this  particular  point,  which  appears  to  be 
in  one  of  the  more  critical  regions  for  separation,  should 
not  separate.  The  relatively  high  value  of  K  does  indicate 
possible  separation  farther  along  the  hub;  however,  the 
change  of  velocity  decreases  beyond  this  point.  The  ef¬ 
fects  of  peripheral  speed  and  curvature  on  separation  are 
illustrated  in  Eq.  (84).  The  value  of  I<  is  approximately 
,0092  at  the  point  of  maximum  curvature,  at  a  distance  of 
about  .4  inches  along  the  tip  streamline.  The  velocity 
gradient  at  this  point  is  practically  zero.  Although  the 
value  of  K  is  low,  it  can  be  surmised  that  a  nominal  nega¬ 
tive  velocity  gradient  might  induce  separation. 

6.  The  results  of  the  blade  to  blade  analysis  are  exten¬ 
sions  of  the  more  exact  meridional  plane  analysis.  The 
relative  velocity  distributions  along  streamlines  (Figs. 

13,  14,  and  15)  can  also  be  analyzed  for  local  phenomena, 
which  may  be  aggrevated  by  the  correction,  AW.  This  veloc¬ 
ity  difference  is  used  to  calculate  the  fluid  forces  on  the 
blade  with  the  relations  of  Section  II.  A  presentation  of 
the  type  shown  in  Fig.  16  is  adequate  for  these  calculations. 

An  attempt  was  made  to  indicate  the  correlation  between 
the  trailing  edge  velocities  in  Fig.  14,  in  order  to 
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establish  some  initial  conditions  for  a  wake  analysis. 

The  dashed  line  indicates  a  possible  velocity  distribu¬ 
tion,  but  this  is  only  a  guess.  Figs.  13  and  15  do  not 
demonstrate  this  convergence.  In  addition,  the  velocities 
at  the  trailing  edge  are  calculated  from  extrapolated  data. 
Therefore,  this  type  of  correlation  is  inconclusive. 

S.  Extensions  of  the  Method 

1.  The  meridional  plane  analysis  sufficiently  describes 
the  flow  for  an  incompressible  analysis.  The  results  might 
be  extended  to  include  the  calculations  of: 

a.  Pressure  distribution 

b.  Power  input  or  power  required 

c.  Separation  parameter 

d.  Cavitation  parameter 

e.  Blade  forces 

Most  of  these  calculations  depend  upon  the  results  of  sev¬ 
eral  characteristic  solutions,  and  would  be  included  in 
a  separate  computer  program.  Many  preliminary  calcula¬ 
tions  of  point  functions  (pressure,  enthalpy,  etc. )  can  be 
included  in  the  main  solution. 

2.  The  methods  should  be  modified  to  include  character¬ 
istic  calculations  outside  the  rotor  passage.  This  would 
permit  continuous  solutions  through  rnulti-stage  machines 
and  in  unbladed  passages.  Major  modifications  would  be 
required,  since  a  number  of  decisions  would  be  needed  to 
insure  the  use  of  proper  blade  properties  in  the  vicinity 
of  blade  boundaries.  Input  requirements  would  necessarily 
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become  more  complex. 

3.  The  theory  and  methods  could  be  modified  to  use  non- 
dimensional  variables.  This  would  be  particularly  applic¬ 
able  to  design  studies,  where  non-dimensional  results  of 
preliminary  rotor  configurations  could  be  directly  compared 

4.  This  method  should  be  extended  to  include  the  compres¬ 
sible  solution.  A  method  is  outlined  in  Ref.  9. 

5.  The  blade  to  blade  method  used  in  this  thesis  is  prac¬ 
tical  in  its  simplicity,  and  is  sufficiently  accurate  when 
used  in  conjunction  with  the  meridional  plane  analysis. 

It  would  be  interesting  to  compare  the  results  of  this  solu 
tion  with  a  more  exact  approach,  such  as  the  method  of  Ref. 

6.  The  extensions  to  the  meridional  plane  analysis  should 
also  be  included  in  the  blade  to  blade  solution. 
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V.  CONCLUSIONS  AND  RECOMMENDATIONS 


A.  It  is  concluded  that: 

1.  Theoretical  assumptions  are  based  on  practical 
consideration,  and  do  not  impose  excessive  limitations  on 
the  solution. 

2.  The  meridional  plane  analysis  is  developed  from 
sound  theoretical  derivations,  which  provide  a  relatively 
rapid  and  accurate  method  of  solution. 

3.  The  methods  which  are  applied  in  solving  this 
problem  are  good  approximations.  The  degree  of  accuracy 
of  each  characteristic  solution  depends  on  the  size  of  the 
coordinate  grid  and  on  the  accuracy  of  data  measurements. 

4.  The  progressive  method  employed  in  the  meridional 
plane  solution  can  be  applied  at  arbitrary  intervals  to 
obtain  the  desired  accuracy. 

5.  The  accuracy  of  the  results,  based  primarily  upon 
theoretical  considerations,  is  sufficient  for  engineering 
purposes. 

6.  Results  are  not  sufficiently  accurate  to  compute 
detailed  analyses  at  the  trailing  edge. 

7.  The  results  of  this  method  are  sufficient  to  com¬ 
pletely  describe  the  flow  within  the  rotor,  and  predict  the 
performance  of  the  machine. 

8.  This  solution  is  restricted  to  flows  within  the 
rotor  passage. 

9.  The  computer  program,  ROTOR  1,  can  be  used  if 
Q=constant  curves  are  accurately  specified. 
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10.  The  use  of  a  direct  flow  analysis  is  prefer- 
rable  to  the  inverse  method  in  design  studies. 

11.  The  blade  to  blade  analysis  is  limited  by  the 
assumption  of  a  linear  velocity  distribution,  but  the 
results  are  acceptable  for  engineering  purposes. 

3.  It  is  recommended  that: 

1.  The  method  be  non-dimensional ized. 

2.  Computations  be  added  to  describe  pressure  dis¬ 
tribution,  blade  force,  power,  and  critical  flow  phenomena. 

3.  The  solution  be  modified  to  include  compressible 
flows,  and  calculations  beyond  rotor  boundaries. 

4.  A  more  accurate  blade  to  blade  analysis  be  used 
in  conjunction  with  this  meridional  plane  solution. 
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TABLE  I 

INITIAL  DATA  FOR  TEST  IMPELLER 

No.  of  Blades . 23 

Blade  Thickness,  in . 3 

RPM . 1800 

Volumetric  Flow  Rate,  cu.  ft. /sec . 128.212 

Annulus  Area  at  L.E.,  sq„  in . 186.22 

Meridional  Velocity,  VmL,  ft. /sec . 94.09 

Air  at  Standard  Temperature  and  Pressure 
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TABLE  U-  1 


CATA  FDR  CHARACTERISTIC  Cl 

LOCATION  OF  CHARACTERISTIC  CURVE 
STARTING  NORMAL  NO.  =  4 


M  -  1  2  3  4  5 


DNSTAR  = 

.OOCO 

-.0342 

-. 0778 

-. 1268 

-.1842 

GAM-LAM  = 

-11. 3400 

-1C. 9150 

-10.4652 

-10.4767 

-  1C. 1 7  C  1 

GAMMA 

-2 . 59C0 

-3.  7924 

-5.2378 

-5.7088 

-7.3777 

BETA (  P )  = 

-49.30C0 

-49.5236 

-49.7943 

-50.0899 

-50.4328 

RADIUS  = 

6. 0100 

6.6160 

7.1639 

7.6701 

8.1727 

LOCATION  OF  NEW 

STREAMLINES 

M  = 

1 

2 

3  4 

5 

DELTA  X  = 

.OOCO 

.0443 

.0805  .0220 

-.01 66 

BET  A ( X  )  = 

-49.30CO 

-49.5414  - 

49.8367  -50.1027 

-50.4215 

RADIUS  = 

6. 01  CO 

6.6595 

7.2431  7.6917 

8. 1563 

VELOCITY 

PROFILE 

M  = 

1 

2 

3 

4 

5 

WM  ( P  ) 

86.4590 

88. 7083 

90.0108 

90.6751 

90.81  16 

WM  (  X ) 

86.4590 

88.  8395 

90.1501 

90.6938 

90.8169 

REL  VEL  = 

132.5857 

136. 9081 

139.7742 

141.3968 

142.5352 

DW  COEF  = 

.1281 

.  1320 

.1345 

.  1359 

.  1368 

DW  FUNC  = 

-44  1  . 

35. 

605. 

1  139. 

1787. 

5  Ar 


TABL  E  TI-  \ 


DATA  FOR  CHARACTERISTIC  Cl 


LOCATION  OF  CHARACTERISTIC  CURVE 
STARTING  NORMAL  NO.  =  4 


M  = 

6 

7 

8 

9 

DNSTAR  = 

-.2471 

-.3166 

-.3881 

-.4679 

GAM-LAM  = 

-9.8294 

-9.7783 

-9.5601 

-9.5577 

GAMMA 

-8.5532 

-9.00  37 

-  10.8228 

-12.3382 

BETA { P  )  = 

-50.8637 

-51.3793 

-51.8894 

-52.4967 

RADIUS  = 

8.6200 

9.0700 

9.4878 

9.9023 

LOCATION  OF  NEW 

STREANL 

INES 

M  = 

6  7 

8 

9 

DELTA  X  = 

-.0664  -.0935 

-.1196 

.0000 

BETA ( X )  = 

-50.7988  -51.2735  - 

51.7424 

-52.4967 

RADIUS  = 

8.5526  8.9779 

9.3699 

9.9023 

VELOCITY 

PROFILE 

M  = 

6 

7 

8 

9 

WM(P) 

90.3799 

89.4412 

88.1390 

86.3790 

WM(X) 

90.4889 

89.6787 

88.5868 

86.3790 

REL  VEL  = 

143.  1683 

143.3472 

143.0669 

141.8821 

DW  COEF  = 

.  1370 

.  1368 

.  1363 

.1349 

DW  FUNC  = 

2401  . 

3145. 

3917. 

5  108. 
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TABLE  IT-2 


DATA  FCR  CHARACTERISTIC  C2 


LOCATION  OF  CHARACTERISTIC  CURVE 
STARTING  NORMAL  NO.  =  3 


M  = 

1 

2 

3 

4 

5 

DNSTAR  = 

.  COCO 

.0221 

.0311 

.0320 

.0326 

GAM-LAM  = 

-9.4547 

-8.8207 

-7. 1870 

-7.  1506 

-5.6883 

GAMMA 

2.2953 

1  .6850 

.1174 

.1087 

.0426 

BETA ( P )  = 

-43.70C0 

-43.6597 

-43.6947 

-43.7616 

-43.8372 

RADIUS  = 

6. 21  CO 

6.8346 

7.3839 

7.8343 

8.2730 

LOCATION  OF  NEV* 

STREAMLINES 

M  = 

1 

2 

3  4 

5 

DELTA  X  = 

.0000 

.0691 

.1202  .1068 

.1069 

BETA ( X )  = 

-43.70C0 

-43.6548 

-43.7039  -43.7778 

-43.8600 

RADIUS  = 

6. 2100 

6.9029 

7.5031  7.9402 

8.3813 

VELOCITY 

PROF  ILE 

M  = 

1 

2 

3 

4 

5 

WM(P) 

77.3222 

81.5205 

84.3795 

86.1041 

87.3981 

WM  ( X ) 

77.3222 

81.9309 

84.8900 

86.4604 

87.6684 

REL  VEL  = 

106.951 1 

113. 2405 

1 17.4265 

1 1 9.7464 

121.5882 

DW  COEF  = 

.  1491 

.1542 

.1575 

.1595 

.1612 

DW  FUNC  = 

1763. 

2506. 

3307. 

3990. 

4768. 
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table  ir-2 


DATA  FCR  CHARACTERISTIC  C2 


LOCATION  OF  CHARACTERISTIC  CURVE 
STARTING  NORMAL  NO.  =  3 


M  = 

6 

7 

8 

9 

DNSTAR  = 

.0329 

.  0303 

.0217 

.0048 

GAM-LAM  = 

-4.7525 

-5.31 10 

-6.0772 

-5.9846 

GAMMA 

.0575 

-.7670 

-1.7137 

-2.5034 

BETA { P  )  = 

-43.9161 

-44.0091 

-44 . 1269 

-44.2784 

RADIUS  = 

8.6330 

9.0528 

9.4521 

9.9105 

LOCATION  OF  NEW 

STREAML 

INES 

M  = 

6 

7 

8 

9 

DELTA  X  = 

.1141 

.1013 

.0315 

.0000 

BET  A ( X )  = 

-43.94CC 

-44.0321 

-44.1350 

-44.2784 

RADIUS  = 

8.7467 

9. 1537 

9.4834 

9.9105 

VELOC  ITV 

PROFILE 

M  = 

6 

7 

8 

9 

WM  (  P ) 

88. 1297 

98.5673 

88.5827 

88.0851 

WM  { X  ) 

88.3047 

88.6192 

88.5680 

88.0851 

REL  VEL  = 

122.6340 

123.2617 

123.4050 

123.0314 

DW  COEF  = 

.  1624 

.  1636 

.1644 

.  1653 

DW  FUNC  = 

5489. 

6383. 

7173. 

8299. 
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TABLE  U>'3 


CATA  FOR  CHARACTERISTIC  C3 


LOCATION  OF  CHARACTERISTIC  CURVE 
STARTING  NORMAL  NO.  =  2 


M  = 

1 

2 

3 

4 

5 

DNSTAR  = 

.ooco 

.0728 

.  1  355 

.  1760 

.  22  C6 

GAM-LAM  = 

-11.0921 

-10.3759 

-9.5373 

-9.0372 

-7.8418 

GAMMA 

5. 6079 

5.9430 

5.9C27 

5.5052 

5.98C3 

BETA ( P )  = 

-38 . 20  CC 

-38.2507 

-38.3586 

-38.5683 

-38.7540 

RADIUS  = 

6 . 47  CO 

7. 1808 

7.7735 

8.  1744 

8.6075 

LOCATION  OF  NEW  STREAMLINES 


M  = 

1 

2 

3 

4 

5 

DELTA  X  = 

.oocc 

.  0072 

.0242 

.0569 

.  C480 

BET  A ( X )  = 

-38.20CC 

-38.2513 

-38.3650 

-38.5950 

-38.7793 

RADIUS  = 

6.47CC 

7.1879 

7.7973 

8.2306 

8.6550 

VE LOCI  TV 

PROFILE 

M  = 

1 

2 

3 

4 

5 

WM  ( P  ) 

70.90  1C 

76.6468 

80.8738 

83.2348 

85.4015 

WM  {  X ) 

70.90  1C 

76.  7004 

8 1.0299 

83.5263 

35.62C7 

REL  VEL  = 

90.2212 

97. 6698 

103.3449 

106.8691 

109.83 14 

DW  COEF  = 

.  1677 

.  1722 

.1752 

.  1766 

.  1778 

DW  FUNC  = 

3559. 

4523. 

5458. 

6  185. 

6976  . 

TABLE  H-  3 


CATA  FOR  CHARACTERISTIC  C3 


LOCATION  OF  CHARACTERISTIC  CURVE 
STARTING  NORMAL  NO.  =  2 


M  = 

6 

7 

8 

9 

DNSTAR  = 

.2578 

.2997 

.3324 

.3751 

GAM-LAM  = 

-7.0684 

-5.6366 

-5.6294 

-5. 6032 

GAMMA 

6.  1473 

6. 1015 

6.0736 

5.8575 

BETA ( P )  = 

-38.9747 

-39.2207 

-39.5066 

-39.6997 

RADIUS  = 

8.9586 

9.3394 

9.6484 

10.0450 

LOCATION  OF  NEW 

STREAMLINES 

M  = 

6  7 

8 

9 

DELTA  X  = 

.0492  .0318 

.0331 

.0000 

BETA ( X )  = 

-39.0023  -39.2461 

-39.5297 

-39.6997 

RADIUS  = 

9.0074  9.3710 

9.6814 

10.0450 

VELOCITY 

PROFILE 

M  = 

6 

7 

8 

9 

WM  (  P ) 

86.8675 

88.2065 

69.0256 

89.7898 

WM  (  X ) 

87.05C9 

88.3039 

89.1115 

89.7898 

REL  VEL  = 

1 12.0172 

1  14.0234 

1  1  5.5348 

1  16.7006 

DW  COEF  = 

.  1785 

.  1791 

.  1793 

.  1799 

DW  FUNC  = 

7673. 

8441. 

9124. 

1 0034. 
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TABLE  3C-4- 


CATA  FOR  CHARACTERISTIC  C4 


LOCATICN  OF  CHARACTERISTIC  CURVE 
STARTING  NORMAL  NO.  =  2 


M  = 

1 

2 

3 

4 

5 

ONSTAR  = 

.  ooco 

.1156 

.2  132 

.2938 

.  37C9 

GAM- LAM  = 

-15. 1727 

-14. 49  58 

-12.7895 

-11 . 6315 

-  1C. 39C7 

GAMMA 

9.5773 

8. 7551 

1C. 1 163 

10.9301 

10. 85C6 

BETA(P)  = 

-32.00CC 

-32.2773 

-32.6731 

-33. 1846 

-33.7C38 

RADIUS  = 

6 . 99  CO 

7.6851 

8.2632 

8.6847 

9.C74  1 

LOCATION  OF  NEW  STREAMLINES 


M  = 

1 

2 

3 

4 

5 

DELTA  X  = 

.OOCC 

-.0004 

.0228 

.0156 

.  CC86 

BETA ( X )  = 

-32. COCO 

-32.2771 

-32.6935 

-33.2041 

-33.7162 

RADIUS  = 

6.99C0 

7.6847 

8.2854 

8.6999 

9.C825 

VELOCITY 

PROFILE 

M  * 

1 

2 

3 

4 

5 

WM  ( P  ) 

66.5165 

72.9600 

77.6404 

80.8896 

83.6173 

WM  (  X ) 

66.5165 

72.9568 

77.8154 

80.9992 

83.6713 

REL  VEL  = 

78.4348 

86.2903 

92.4644 

96.8051 

100.59C9 

DW  COE  F  = 

.  1866 

.  1903 

.1917 

.  1920 

.  1923 

DW  FUNC  = 

5723. 

6882. 

7974. 

8732. 

9464  . 
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TABLE  31-4- 


DATA  FOR  CHARACTERISTIC  C4 


LOCATION  OF  CHARACTERISTIC  CURVE 
STARTING  NORMAL  NO.  =  2 


M  = 

6 

7 

8 

9 

DNSTAR  = 

.  4394 

.  5076 

.5698 

.  6407 

GAM-LAM  = 

-9.6181 

-9. 4840 

-9.0488 

-8.3707 

GAMMA 

1  1 . 0479 

1 1 . 5039 

1  1 .9500 

12.9363 

BETA ( P  )  = 

-34 . 38  C  2 

-35.1 071 

-35.9127 

-36.7628 

RADIUS  = 

9.4277 

9.7770 

10.0699 

10.3860 

LOCATION  OF  NEW  STREAMLINES 


M  = 

6 

7 

8 

9 

DELTA  X  = 

-.0056 

-.0078 

-.0114 

.0000 

BET  A ( X  )  = 

-34.3697 

-35.0909 

-35.8826 

-36. 7628 

RADIUS  = 

9.4222 

9.7694 

10.0586 

10. 3860 

VELOCITY 

PROFILE 

M  = 

6 

7 

8 

9 

WM  {  P ) 

85.7169 

87.5236 

89.0085 

90.484 1 

WM  (  X ) 

85.6873 

87.4843 

88.9556 

90.4841 

REL  VEL  = 

103.81 14 

106.9175 

109.7919 

1  12.9469 

DW  COEF  = 

.  1918 

.1910 

.  1896 

.1879 

DW  FUNC  = 

10108. 

10784. 

11304. 

11908. 

\ 


6)1 


TABLE  H-5 


CATA  FOR  CHARACTERISTIC  C5 


LOCATION  OF  CHARACTERISTIC  CURVE 
STARTING  NORMAL  NO.  =  2 


M  = 

1 

2 

3 

4 

5 

DNSTAR  = 

.COCO 

.lull 

.2609 

.3466 

.4271 

GAM-LAM  = 

-19.6662 

-19.0712 

-  1 8.391 0 

-16.4404 

-15.7630 

GAMMA 

1  1 . 53  3  8 

1  1.5338 

1  1.3148 

12.0517 

12.8821 

BETA ( P  )  = 

-27.80CC 

-28. 53L3 

-29.3978 

-30.4472 

-31.6027 

RADIUS  = 

7.59C0 

8.2556 

8.8105 

9.2064 

9.5618 

LOCATION  OF  N  Eb 

STREAMLINES 

M  = 

1 

2 

3  4 

c 

DELTA  X  = 

.OOCO 

.  0046 

.0032  .0092 

.  C  4  1  9 

BET  A ( X )  = 

-27.80C0 

-28. 5399 

-29.4043  -30.4730 

-31.7408 

RADIUS  = 

7.59CC 

8.2599 

8.8135  9.2152 

9.6C21 

VELOCITY 

PROFILE 

M  = 

1 

2 

3 

4 

5 

WM(P) 

63.0930 

70.0131 

75.2679 

78.5352 

81  .  1981 

WM  ( X  > 

63.0930 

70.0553 

75.2940 

78.6023 

81  .49  14 

REL  VEL  = 

71  .3253 

79.7456 

86.4278 

9 1 . 1999 

95.8230 

DW  COEF  = 

.  1985 

.2005 

.2012 

.2001 

.1982 

DW  FUNC  = 

7829. 

9084. 

10154. 

10892. 

11571. 

(oZ 


TABLE  XT-  5 


CATA  FOR  CHARACTERISTIC  C5 


LOCATION  OF  CHARACTERISTIC  CURVE 
STARTING  NORMAL  NO.  =  2 


M  = 

6 

7 

8 

9 

DNSTAR  = 

.5103 

.6037 

.6952 

.8083 

GAM-LAM,  = 

-15.2153 

-14.5683 

-14.7957 

-14.1805 

GAMMA 

1 4.4943 

1 5. 7970 

18.0082 

19.5109 

BET  A ( P )  = 

-32.8655 

-34.2872 

-35.8854 

-37.72C8 

RADIUS  = 

9.8964 

10.2519 

10.5742 

10.8850 

LOCATION  CF  NEW  STREAMLINES 


M  = 

6 

7 

8 

9 

DELTA  X  = 

.  03  £9 

-.00  40 

-.0084 

.0000 

BET  A ( X  )  = 

-33.0031 

-34.2713 

-35.8425 

-37.7208 

RADIUS  = 

9.9339 

10.2480 

10.5660 

10.8850 

VELOCITY 

PROFILE 

M  = 

6 

7 

8 

9 

WM(P) 

83.5503 

85.8377 

87.81 1 5 

89.7423 

WM(X) 

83.7991 

85.8 120 

87.  7596 

89.7423 

REL  VEL  = 

99.9224 

103. 8407 

108.2609 

1 13.4539 

DW  COEF  = 

.  1956 

.  1930 

.1890 

.  1  842 

DW  FUNC  = 

12113. 

12605. 

1  3006. 

1 3267. 

TABLE  lT-<o 


CAT  A  FCR  CHARACTERISTIC  C6 


LOCATION  OF  CHARACTERISTIC  CURVE 
STARTING  NORMAL  NO.  =  2 


M  = 

1 

2 

3 

4 

5 

DNSTAR  = 

.ooco 

.  1  4  62 

.27  16 

.3765 

.4888 

GAM-LAM  = 

-28.8995 

-27.4739 

-25.9103 

-24.  1102 

-24.5668 

GAMMA 

1  1 .05C5 

12.8221 

14.0721 

15.8992 

16. 1C83 

BETA(P)  = 

-23. 70CC 

-25.5592 

-27.5695 

-29.7262 

-31 .5938 

RADIUS  = 

8.70C0 

9. 3214 

9.8C34 

10.1755 

10.5488 

LOCATION  OF  NEW  STREAMLINES 


M  = 

1 

2 

3 

4 

5 

DELTA  X  = 

.0000 

-.0620 

-.0409 

-.0278 

-.0712 

BETA { X )  = 

-23. 70C0 

-25. 3956 

-27.4 165 

-29.5774 

-31 . 5970 

RADIUS  = 

8.70CC 

9 . 2  664 

9.7665 

10.1501 

10.4842 

VELOCITY 

PROFILE 

M  = 

1 

2 

3 

4 

5 

WM(P) 

61.6410 

68.8717 

74.0780 

77.6503 

80.5685 

WM  ( X  ) 

61.6410 

68.2392 

73.7027 

77.4126 

80.41 59 

REL  VEL  = 

67.3185 

75.  5386 

83.0281 

89.0113 

94.4121 

DW  COEF  = 

.21  18 

.2102 

.2076 

.2035 

.2000 

DW  FUNC  = 

11442. 

12583. 

13499. 

14068. 

1  4496. 

6>4 


TABLE  IT-O, 


DATA  FCR  CHARACTERISTIC  C6 


LOCATION  OF  CH AR AC T ER I S T I C  CURVE 
STARTING  NORMAL  NO.  =  2 


M  = 

6 

7 

8 

9 

DNSTAR  = 

.5852 

.  6687 

.7762 

.  9022 

GAM-LAM  = 

-23.5355 

-22.0885 

-21 .0109 

-  19.2679 

GAMMA 

17. 5715 

19.0663 

2 1 .3290 

22.7489 

BETA(P)  = 

-34.3762 

-36.8106 

-39.4102 

-42.3090 

RADIUS  = 

10.8624 

11.1154 

11.4175 

11.7516 

L  CCA  T ION  OF  NEW 

STREAMLINES 

M  = 

6 

7 

8 

9 

DELTA  X  = 

-.0661 

.  0530 

.0501 

.0000 

BETA ( X )  = 

-33. 90C9 

-37.2273  - 

39.7999 

-42.3090 

RADIUS  = 

10.8018 

11.1645 

1  1 . 4642 

1  1.7516 

VELOCITY 

PROFILE 

M  = 

6 

7 

8 

9 

WM(P) 

83. 4560 

85.2295 

87.0504 

88.5593 

WM  (  X  I 

82.9763 

35. 5546 

87.3327 

88.5593 

REL  VEL  = 

99.97 1 1 

107.4487 

1  1  3.6722 

119.7514 

DW  COEF  = 

.  1949 

.  1865 

.  1791 

.  1719 

DW  FUNC  = 

1 4766  . 

14786. 

14764. 

14664 . 

(o5 


TABLE  H-l 


CATA  FGR  CHARACTERISTIC  C7 


LOCATION  OF  CHARACTERISTIC  CURVE 


STARTING 

NORMAL  NO. 

=  2 

M  = 

6 

7 

8 

9 

DNSTAR  = 

.  6225 

.7263 

.8381 

.9475 

GAM-LAM  = 

-29.5988 

-27.9648 

-28.2480 

-28.9740 

GAMMA 

18.5752 

20.  1474 

21.3246 

2 1 .9428 

BETA ( P  )  = 

-37.6563 

-40.4622 

-43.3651 

-46.2584 

RADIUS  = 

1  1  .8069 

12.  1087 

12.4042 

12.6870 

LOCATION  OF  NEW 

STREAMLINES 

M  = 

6 

7 

8 

9 

DELTA  X  = 

-.0185 

.0040 

-.0122 

.0000 

BETA { X  )  = 

-37.5028 

-40.4982  - 

43.2526 

-46.2584 

RADIUS  = 

1  1  .  79 C 8 

12.  1  122 

12.3935 

12.6870 

VELOCITY 

PROFILE 

M  = 

6 

7 

8 

9 

WM  (  P  ) 

79.6489 

81 . 5500 

33.2322 

84.4272 

WM(X) 

79. 53CO 

81 .5733 

83.  1  723 

84.4272 

REL  VEL  = 

100.2490 

107.2731 

1  14.1944 

122. 1092 

DW  COE F  = 

.1877 

.  1793 

.17  11 

.1618 

DW  FUNC  = 

1 75 7C. 

17  52  8. 

17315. 

169C9. 

TABLE  TT-7 


CATA  FCR  CHARACTERISTIC  C7 


LOCATION  OF  CHARACTERISTIC  CURVE 
STARTING  NORMAL  NO.  =  2 


M  = 

1 

2 

3 

4 

5 

DNSTAR  = 

.ooco 

.  1578 

.  3039 

.4  101 

.5172 

GAM-LAM  = 

-32.2487 

-33.2812 

-31.5112 

-30.3829 

-29.9231 

GAMMA 

15.0513 

14.L762 

1 5.776C 

16.4842 

16.9608 

BETA  I P )  = 

—2  4 . OOCO 

-26.6517 

-29.3615 

-32 . 0475 

-34.8028 

RADIUS  = 

9. 81  CO 

10. 3442 

10.8314 

11.1799 

1  1  .  5C4  9 

LOCATION  OF  NEW  STREAMLINES 


M  = 

1 

2 

3 

4 

5 

DELTA  X  = 

.OOCC 

-.0096 

-  .04C0 

-.0250 

-.0108 

BETA! X  )  = 

-24. COCO 

-26.6106 

-29. 1 671 

-31.8716 

-54.7224 

RADIUS  = 

9. 81  CO 

10. 3361 

10.7972 

11.1583 

1  1.4955 

VELOCITY 

PROFILE 

M  = 

1 

2 

3 

4 

5 

WM  ( P  ) 

58. 037C 

65. 4298 

71  .  1 C  53 

74.5236 

77.3497 

WM  (  X  ) 

58.0370 

65.3268 

70.7192 

74.3173 

77.2744 

REL  VEL  = 

63.5294 

73.0666 

80.9884 

87.5111 

94.0167 

OW  COE  F  = 

.2129 

.2  107 

.2065 

.2012 

.  1949 

DW  FUNC  = 

1 5098. 

16C79. 

16861. 

17282. 

17522. 

(ol 


TABLE  H  -& 


DATA  FOR  CHARACTERISTIC  C8 

LOCATION  OF  CHARACTERISTIC  CURVE 
STARTING  NORMAL  NO.  =  2 


M  =  1  2  3  4  5 


DNSTAR  = 

.0000 

.1863 

.3537 

.4777 

.5875 

GAM-LAM  = 

-38.3072 

-36. 3681 

-36.9234 

-36.9877 

-37.0683 

GAMMA 

17.4928 

18.9027 

18.1034 

18. 1761 

17.5196 

BETA(P)  = 

-26.4000 

-29. 3276 

-32.2429 

-35.0748 

-37.9632 

RADIUS  = 

10.9200 

1 1.3957 

1  1.8267 

12.1581 

12.4577 

LOCATION  OF  NEW 

STREAMLINES 

M  = 

1 

2 

3 

4 

5 

DELTA  X  = 

.0000 

.0100 

.0013 

-.0190 

-  .0090 

BETA(X)  = 

-26.4000 

-29.3831  - 

■32.2524 

-34.9396 

-37.8908 

RADIUS  = 

10.9200 

1 1.4037 

1  1.8277 

12.1429 

12.4505 

VELOCITY 

PROFILE 

M  = 

1 

2 

3 

4 

5 

WM  I  P  1 

52.9880 

61.1133 

67.2976 

71.1384 

73.8746 

WM  (  X ) 

52.9880 

61.2430 

67.31 12 

70.9712 

73.81  68 

REL  VEL  = 

59.1574 

70. 290? 

79.5913 

86.5760 

93.5358 

DW  COEF  = 

.21  10 

.2050 

.2013 

.  1955 

.1885 

DW  FUNC  = 

19031 . 

19793. 

20341. 

20569. 

20637. 

68 


TABLE  H-e 


DATA  FOR  CHARACTERISTIC  C8 


LOCATION  OF  CHARACTERISTIC  CURVE 
STARTING  NORMAL  NO.  =  2 


M  = 

6 

7 

8 

9 

DNSTAR  = 

.6883 

.7985 

.8831 

.9746 

GAM-LAM  = 

-37.0218 

-37. 5553 

-38.2758 

-40.9817 

GAMMA 

18.8401 

19.8310 

19.3823 

17.7669 

BETA(P)  = 

-40.7806 

-43.7543 

-46.4951 

-49.3487 

RADIUS  = 

12.7323 

12. 9993 

13.2607 

13.5496 

LOCATION  OF  'JEW  STREAMLINES 


M  = 

6 

7 

3 

9 

DELTA  X  = 

.0038 

.  0275 

.0503 

.0000 

BETA ( X )  = 

-40.8126 

-44. 0789 

-46.9727 

-49.3487 

RADIUS  = 

12.7354 

13.0215 

1 3. 3002 

13.5496 

VELOCITY 

PROFILE 

M  = 

6 

7 

8 

9 

WM(P) 

75.7367 

77. 3583 

78.3727 

79. 1673 

WM( X )  = 

75.7561 

77.485+ 

78.5444 

79. 1673 

REL  VEL  = 

1 00.0937 

107. 8603 

1  15. 1 C  9 1 

121.5238 

DW  COEF  = 

.  1802 

.  1703 

.1614 

.  1539 

DW  FUNC  = 

20574. 

20237. 

19914. 

1961 5. 

TABLE  H-9 


DATA  FOR  CHARACTERISTIC  C9 


LOCATION  OF  CHARACTERISTIC  CURVE 
STARTING  NORMAL  NO.  =  2 


M  = 

1 

2 

3 

4 

5 

DNSTAR  = 

.ooco 

.  1820 

.3319 

.4518 

.5673 

GAM-LAM  = 

-45 .4285 

-45.  3287 

-44  .7492 

-45.0907 

-46.5900 

GAMMA 

18.37  15 

17.8140 

18.3051 

18.6397 

17.2144 

BETA ( P )  = 

-31 .80C0 

-34.7985 

-37.6903 

-40. 3692 

-43.0088 

RAOIUS  = 

12. 42CO 

12.8306 

13.1770 

1 3.4616 

13.7247 

LOCATION  OF  NEW  STREAMLINES 


M  = 

1 

2 

3 

4 

5 

DELTA  X  = 

.OOCO 

-. 0607 

-.0757 

-.0746 

-.1057 

BETA ( X )  = 

-31.8000 

-34.4885 

-37.2367 

-39.8407 

-42.2680 

RADIUS  = 

12.42C0 

12.7880 

13.1233 

13.4089 

13.6521 

VELOCI  TY 

PROFILE 

M  = 

1 

2 

3 

4 

5 

WM  (  P ) 

55 .  1540 

61.8969 

66. 3989 

69.3353 

7 1  . 64  4  1 

WM  (  X  ) 

55. 1540 

61.2727 

65.7637 

68.8022 

71 .1038 

REL  VEL  = 

64. 8952 

74.3385 

82.6028 

89.6061 

96.0853 

DW  COEF  = 

.2036 

.  1982 

.1915 

.  1  845 

.  1779 

DW  FUNC  = 

2398C . 

24365. 

24591  . 

24654. 

24544. 
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TABLE  H- 9 


DATA  FOR  CHARACTERISTIC  C9 


LOCATION  OF  CHARACTERISTIC  CURVE 
STARTING  NORMAL  NO.  =  2 


M  = 

6 

7 

8 

9 

ONSTAR  = 

.6681 

.  7596 

.8385 

.  9061 

GAM-LAM  = 

-47.61  LI 

-47.3437 

-45.6471 

-44.5153 

GAMMA 

16.3507 

1 6.3615 

16.4971 

16.7489 

BETA(P)  = 

-45 .5308 

-47.9953 

-50.2033 

-52. 3203 

RADIUS  = 

13.9669 

1  4.  2071 

14.3925 

14.5853 

LOCATION  OF  NEW 

STREAMLINES 

M  = 

6 

7 

8 

9 

DELTA  X  = 

-.  1353 

-.  1049 

-.0744 

.  0000 

BETA  I  X )  = 

-44.5559 

-47. 1992 

-49. 61 45 

-52.3203 

RADIUS  = 

13.8757 

14.  1360 

14.3405 

14.5853 

VELOCITY 

PROFILE 

M  = 

6 

7 

8 

9 

WM  (  P ) 

73.2304 

74.4063 

75.2378 

75.7349 

WM  ( X ) 

72.6960 

74.06  11 

75.0393 

75.7349 

REL  VEL  = 

102.0200 

109. 0013 

115.8142 

123.9022 

DW  COEF  = 

.17  11 

.  1628 

.  1  546 

.1450 

DW  FUNC  * 

24374. 

24100. 

23583. 

22936. 

T  A  GL  E  HL-  | 
VELOCITY  PROFILE 


MERIDIONAL  STREAMLINE  NO.  1 


M 


.00 
.50 
1 .00 

1.50 

2.00 

2.50 
3.00 

3. 50 
4.00 

4.50 

5 .00 

5.50 
6.00 
6.  50 
7.00 

7.50 

8.00 

8.50 
9.00 

9.50 

10.00 

10.50 
11.00 
1  1 .50 
12.00 

12.50 
13.00 

13.50 
14.00 


W 


160. 2000 
1 44. 70  CO 
130.20CC 
1  1  7. COCO 
106.80CC 

97.2000 
91 .3000 
86.40CC 
82.4000 

78.8000 
75. 60C0 
72. 9CCC 
70.9000 
69.6000 

68. 5000 
67.50CC 

66.8000 
65.8000 

64.5000 
63.20CC 
61 . 70CC 
60.0000 

59.2000 
60. COCO 

62.2000 

65. 5000 
69. COCO 
73. COCO 
77. COCO 


D  W 


10.  1043 
1 1. 3940 
1 1 .4028 
1 1.8866 
11. 5423 
10.  9433 
11.1128 
10. 5962 
10. 4334 
1  1.  6713 
12.6603 
13. 3047 
14. 2076 
15.  1219 
16.5328 
1  8.  302  1 
19. 1616 
19. 6617 
21. 3834 

23.  592  7 
2u. 9452 

24.  954  1 
24. 9375 
24.938  1 
24. 7177 
23.  9489 
23. 0312 
23.  660  1 
22. 2040 


W  SUCTICN 


W  PRES  SIRE 


170. 
156. 
141  . 
128. 
1  18. 
108. 
102. 
96. 
92. 
90. 
88. 
86. 
85. 
04. 
85. 
85. 
85. 
85. 
P  p 


86. 

86. 

84. 

84. 

84. 

86. 

89. 

92. 

96. 

99. 


3043 
0940 
6028 
8866 
2423 
1433 
4128 
9962 
8334 
47  1  3 
2603 
2047 
1076 
72  19 
0328 
802  1 
9616 
46  17 
8834 
7927 
6452 
954  1 
1375 
938  1 
9177 
4489 
03  13 
660  1 
2040 


150. 
123. 
1  18. 

1  C  5  . 
95. 
86. 
80. 
75. 
71. 
67. 
62 . 
59. 
56. 
54. 
51. 
49. 
47. 
46. 
43. 
39. 
36. 
35. 
34. 


3  C 

3? : 

4  1. 

45. 

49. 

54. 


0957 

2C6C 

7972 

1134 

2577 

2567 

1872 

8037 

9666 

1287 

93  9  7 

5953 

6924 

4781 

9672 

1979 

6384 

1383 

1166 

6073 

7548 

0459 

2625 

0619 

4823 

5511 

9688 

3399 

7960 
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TABLE  HT -2 


VELOCITY  PROFILE 

MERIDIONAL  STREAMLINE  NO.  2 


M 

W 

D  to 

W  SUCTICN 

V»  PRESSURE 

.00 

1 56 . 70 CO 

1  1  .  2554 

167.9554 

145.4446 

.50 

1 43.90C0 

1  1. 8943 

155.7948 

132.C052 

1  .00 

133.00CC 

12.7219 

145.7219 

120.2781 

1.50 

121 .COCO 

13.  24  14 

134.2414 

1 C  7 . 7586 

2.00 

111. 40  CO 

13.  3102 

124.7102 

98. C898 

2.50 

103. 10C0 

1  3.  1  943 

1 16.2943 

89.9057 

3.00 

97. 1CCC 

12. 5589 

1 09.6589 

84.54  1  1 

3.50 

92.40C0 

12. 3531 

104.7531 

8C.C469 

4.00 

88.60CC 

13.  C22 5 

1C  1 .6225 

75.5775 

4.50 

85.5000 

13. 7728 

99.2723 

7  1 .7272 

5.00 

82.7000 

14. 2125 

96.9125 

68.1)875 

5.50 

80 . 30CC 

14.7261 

95.026 1 

65.5739 

6.00 

78. 50C0 

15. 5826 

94.0826 

62.9174 

6.50 

77.20CC 

1  6.  64  9  8 

93.8498 

60. 55C2 

7.00 

76.2CCC 

18. 3017 

94.5017 

57.8983 

7.50 

75.40CC 

19. 5005 

94.9005 

55.8995 

8.00 

74. 60CC 

19. 6812 

94.2813 

54.9187 

8.50 

73.80CC 

20. 22C2 

94.0202 

53.5798 

9.00 

72 . 9CC0 

2  1.  5  4  7  C 

94.4470 

5  1  .  353C 

9.50 

72.0000 

23.  C 5 5 6 

95.0556 

48.9444 

10.00 

71  .OOCO 

23.  4652 

94.4652 

47. 5348 

10.50 

70.20CC 

2  3.  5125 

93.7125 

46. 6875 

1  1.00 

70.40CC 

23.  6 Li  1  C 

94.0410 

46 .759C 

11.50 

71  .LOCO 

23. 3800 

94.7800 

48.C2CC 

12.00 

74. 30C0 

23.  1866 

97.4866 

5 1 . 1 1 34 

12.50 

78.20CC 

22. 2397 

ICC. 4397 

55.96C3 

13.00 

83. COCO 

22. 5898 

105.5898 

6C. 4102 
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TABL  E  HD"  3 


VELOCITY  PROFILE 

MERIDIONAL  STREAMLINE  NO.  3 


M 

W 

D  Vs 

W  SUCTICN 

W  PRESS LRE 

.00 

1 55.20CC 

1  3. 2  4  L  4 

1 68.4444 

141 . 9556 

.50 

143. 80CC 

12. 9910 

156.7910 

1 3  C . EC9C 

1  .00 

1 33.30C0 

1 4. C6CC 

147.3600 

1 19 .24CC 

1  .50 

123. 30CC 

14. 6491 

137.9491 

1  C  8  .  6  5  C  9 

2.00 

1 14.50C0 

1  4.  6679 

129.  1679 

99.8321 

2.50 

107. 50CC 

14. 3500 

121.8500 

93.  1  5CC 

3.00 

102.00CC 

1  3.  9765 

115.9765 

88.0235 

3.50 

97. 70C0 

1 4. C631 

1  1  1.7631 

83.6369 

4.  00 

94.  10C0 

14. 3938 

1  C  8 . 49  38 

79.7062 

4.50 

90. 70C0 

14. 4722 

1C5.  1722 

76.2278 

5.00 

87. 90CC 

1  4.  7451 

1C2 .645  1 

73.  1549 

5.50 

85.80CC 

1 5. 571 3 

1 C 1 .3713 

7C.2287 

6.00 

84.60C0 

16. 3908 

ICC. 9903 

68.2C92 

6.50 

83. 70C0 

17. 1630 

ICO. 8630 

66. 537C 

7.00 

83. 10CC 

18.  9291 

102.0291 

64.  1  7C9 

7.50 

82. 30CC 

19. 35C4 

1  Cl .6504 

62. 9496 

8.00 

81 .60C0 

18. 7687 

ICO. 3687 

62.8313 

8.50 

81 .OOCO 

20. 4359 

1C  1  .  4359 

6C.5641 

9.00 

80.50C0 

22. 4865 

1C2.9865 

58.C135 

9.  50 

80. OOCO 

23.  2  1  2 e 

103.2128 

56.7872 

10.00 

79. 50CC 

22. 5787 

1C2.0787 

56.9213 

10.50 

79.70C0 

22. C  30  1 

1C  1 .7301 

57.6699 

11.00 

80. 60C0 

22.  0323 

1C2.6323 

58.5677 

1  1.50 

02.5CCC 

21. 7964 

1 C4 .2964 

6C . 7C36 

12.00 

85.70C0 

21. 6038 

1C  7. 30  38 

64 .C962 

12.50 

90. OOCO 

20. 9741 

1  1C. 9741 

6  9. C  25  9 
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TABLE  HL-4- 


VELCCITY  PROFILE 


MERIblCNAL  STREAKLINE  NC.  4 


M 


W 


D  V\ 


W  SLCTICN  W  PRESSURE 


.00 

.50 

1.00 

1.50 
2.00 

2.50 
3.00 

3.50 
4.00 

4.50 
5.00 

5.50 
6.00 

6.50 
7.00 
7.5C 
8.00 

8.50 
9.00 

9.50 
10.00 

10.50 
11.00 

11.50 
12.00 


1 54. 60CC 
143.40CC 
1  32.80CC 
123. 20 CO 
1 1 6. OOCO 
1 10.C0CC 
1 05. OOCO 
101 .  10CC 
97.60C0 
94 .60CC 
92. 30  C  C 
90.80C0 
89. 90C0 
89. 20C0 
88. 70C0 
88. 30C0 
87. 80CC 
87.50CC 
87. COCO 
86.70CC 
86. 60C0 
87 .30CC 
88. 80CC 
91 . 80C0 
96. 10C0 


1 2. 6731 
1 3. 3948 
14. 5162 
15.  1 5  C  5 

1  5.  2259 
14. 8893 
14. 2543 
14. 4224 
15.  0151 
15. 3777 
15. 5216 
15.7555 
16.2322 
1 6. 8063 
17. 6029 
18. 6923 
19.1391 
20. 2892 

2  1.  72L7 
22. C298 
21. 6809 
21.2212 
21.1551 
20.  8495 
20. 8327 


167.275  1 
156.7948 
147.3162 
138.3505 
13  1 .2259 
124.8893 
1  19.2543 
1  1  5 .5224 
112.6151 
1C9.9777 
1  07.82  16 
1C6.5555 
1C6. 1322 
1  06 . 0063 
1  C6 . 3029 
1  06 . 9923 
106.939  1 
107.7892 
1C8.7247 
1 C  8 . 7  298 
1C8 .2809 
108.5212 
1C9.9551 
1  12.6495 
1 16.9327 


141 

130 

118 

ice 

ICC 
95 
9  C 
86 
82 


92  6  9 
CC5  2 
2837 
0495 
774  1 
1107 
7457 
6776 
5  E4  9 
79.2223 
76.7784 
0445 
6678 
3937 
C  97  1 
6C77 
6609 
21C8 
2753 
6702 
919  1 
0788 
6449 
9505 
,2673 


75  . 
73. 
72. 
71  . 

69. 
68. 
67  . 
65. 
64  . 
64. 
66  , 
67, 

70, 
75 


75 


TABL  E  UL-^ 


VELOCITY  PROFILE 


MERIDIONAL  STREAMLINE  NO.  5 


M 

W 

DW 

W  SECTION 

W  PRESSURE 

.00 

152. 70CC 

1  3.  C122 

165.7122 

139.6878 

.50 

143. COCO 

14. 6320 

157.6320 

128.368C 

1 .00 

1 33. 70CC 

15. 3426 

149.0426 

118.3574 

1.50 

125.00CC 

15.  7981 

140.7981 

1C5.2C15 

2.00 

1 1 8.2CCC 

1  5.  8833 

134.C833 

1C2.3167 

2.50 

112. 20CC 

15.  5870 

127.7870 

56.6130 

3.00 

107. 40C0 

14. 7262 

122.  1262 

52.6738 

3.50 

103. 70CC 

14.  3C23 

1  18.0023 

85.3577 

4.00 

1 0C.80CC 

14. 7846 

1  1  5.5846 

86.C154 

4.50 

98. 30C0 

15. 0312 

113.3312 

83.2687 

5.00 

96.  40  CO 

1  5.  2025 

1 1 1.6025 

8  1  .  1975 

5.50 

95.40C0 

1  5.  4751 

1 10.8751 

79.5245 

6.00 

95. COCC 

15. 4559 

1  10.4599 

75.54C 1 

6.50 

94. 70CC 

15.  8551 

1 1 C. 555  1 

78.8449 

7.00 

94. 40CC 

15.  8689 

1  1C.  2689 

78.5311 

7.50 

94. 1 OCC 

1  6.  9366 

1 1 1 .0366 

77  .  1  634 

8.00 

93 .9CCC 

1  9.  5536 

1 13.4936 

74 . 3C64 

8.50 

93.80C0 

20. 7099 

114. 5099 

73. C9C 1 

9.00 

93.  60  C  C 

20. 6857 

1 14.2857 

72.5143 

9.50 

93.80C0 

20.  52C8 

1  14.3208 

73.2752 

10.00 

94 . 30  C  C 

20. 43CC 

1 14. 73C0 

73. 87CC 

10.50 

95. 20C0 

21. 4977 

1 16.6977 

73.7C23 

11.00 

96. 60C0 

22.  9335 

1  19.5335 

73.6665 

1  1.50 

99. 10CC 

15.  7829 

114.8829 

83.3171 

T  A  PL  E  m-  (o 


VELHC I  TV  PROFILE 

MERIDIONAL  STREAMLINE  NO.  6 


M 

w 

D  W 

W  SUCTICN 

W  PRESSURE 

.00 

151  . 30CC 

1  4.  9754 

166.2754 

136. 3246 

.  50 

142.  10C0 

15. 0663 

1  57  .  1663 

127.C337 

1.00 

133. 30C0 

15.  7379 

149.0379 

117.5621 

1.50 

1 26 . 00  CO 

16. 2337 

142.2337 

1C9. 7663 

2.00 

1 18.90C0 

1  6.  3429 

135.2429 

1C2. 5571 

2.50 

1 13.80C0 

15.  7706 

129.5706 

96 . C294 

3.00 

109. 60CC 

1  3.  8396 

123.4396 

95. 76C4 

3.50 

106.50CC 

1  3.  1  38C 

1  19 .638C 

93 . 362  C 

4.00 

103.70CC 

13. 9661 

117.6661 

89.7339 

4.50 

101 . 5  0  C  0 

1  3.  9572 

1  15.4572 

e?  .5423 

5.00 

100. 40C0 

14.  1321 

114.5321 

86.2679 

5.50 

99. 80C0 

14. 4612 

114.2612 

85.5388 

6.00 

99. 50CC 

14. 5626 

1  14.C626 

84.9374 

6.  50 

99.50C0 

14. 8687 

114.3688 

84.6312 

7.00 

99. 50C0 

1  5.  6544 

115.  1544 

83.8456 

7.50 

99. 50CC 

1  6.  4922 

1 15.9922 

83.CC7S 

8.00 

99. 50CC 

17. 5373 

1 1 7.0373 

8  1  .  9627 

8.50 

99 . 70  CO 

1  7.  8737 

1 17.5737 

81.8263 

9.00 

100. OOCO 

18. C790 

1 18.0790 

81.921C 

9.5C 

100.40CC 

19. 4753 

119. 8753 

6C. 9247 

10.00 

101 .CCCC 

2  1.  1  58C 

122  .  1580 

79 .842C 

10.  50 

1C1 .9CC0 

24. 7289 

126.6289 

77.1711 

11.00 

103.00CC 

39.  8627 

142.8627 

63  .  1  373 

17 


TABLE  431  ~~f 


VELOCITY  PPCFILE 

MERICICNAL  STREAMLINE  NC.  7 


M 

w 

nw 

W  SUCTICN 

Vi  PRESSURE 

.00 

149.80CC 

7. 2847 

157.0347 

142.5652 

.50 

1 40. 8CCC 

20. 6986 

161.4986 

120.1014 

1.00 

132.40C0 

17. 4U89 

148.8469 

114.9511 

1 .50 

125.00C0 

1  5.  4627 

140.4627 

109.5373 

2.00 

1 19.50CC 

1 4. 91C4 

134.41C4 

104. 5886 

2.50 

1  15.30CC 

1 3. 2504 

128. 5504 

1 02 . 0496 

3.00 

1  1  1 .5000 

13. 9944 

125.4944 

87.5056 

3.50 

1 08.90CC 

15. 0027 

122.9037 

82.8963 

4.0C 

1 06 .8000 

13. 7  1  19 

120.5119 

92.0861 

4.50 

1 05.20CC 

12. 6254 

1 17.8254 

82 . 5746 

5.00 

104.0000 

12.  0089 

1  16.0089 

81.8911 

5.50 

104. 60CC 

1  1.  5322 

116.1 333 

92.0667 

6.00 

106.20C0 

1  2.  CCC  1 

1 18.2001 

84.  1899 

6.50 

1 07.20C0 

13.  1771 

120.3771 

84.C229 

7.00 

1 07.50C0 

14. 2603 

1 2 1 .7603 

82.2287 

7.50 

107. 50  CC 

16.  1911 

123.6911 

91 . 2088 

8.00 

107 .5000 

1 7. 8708 

125.4708 

89.5282 

8.50 

1 07.70CC 

17.  8189 

125.5199 

88.6801 

9.00 

108. COCO 

1  7.  7862 

125.7962 

80 . 2C38 

9.50 

1 08 . 30C0 

19. 2902 

127.6903 

68.8097 

10.00 

108.70C0 

22.  121 4 

120.8214 

86.5736 

10.50 

1 09.20CC 

30. 7063 

139.9068 

78.4932 

78 


TABLE  HI  ~  B 


VELOCITY  PROFILE 

MERIDIONAL  STREAMLINE  NO.  8 


M 

W 

D  h 

1*  SUCTICN 

V,  PRESSURE 

.00 

147.80CP 

1 5. 59C6 

163.3906 

122.2CS4 

.50 

139.20CC 

1  4.  9572 

154. 1572 

124.2428 

1.00 

131 . 30CC 

1 5. 7562 

1 47.0562 

1  1  5  .  5438 

1.50 

124.20CC 

15. 2227 

1 39.4227 

108.9773 

2.00 

1  19. 90CC 

1 3. 8629 

132.7629 

1C6.C371 

2.50 

1  16. 4CC0 

1 1. 8229 

128.2239 

1C4.5761 

3.00 

1 13. 30C0 

12. 9087 

126.2087 

1 CC . 29  1  2 

3.50 

111. 80C  0 

1  4.  5964 

126.3964 

97.2036 

4.00 

109.40C0 

1  3.  0625 

122.4625 

96. 2375 

4.50 

108. 60C0 

1  1.  7227 

12C.3237 

96.8762 

5.00 

105.80C0 

10. 2CC4 

1  16.0004 

95.5996 

5. 50 

109.80CC 

S.  4688 

1 18.2688 

1C  1  .2312 

6.00 

1  12  .COCO 

9. 2796 

121.2796 

1C2.72C4 

6.50 

1 13.30CC 

1  1. 9992 

125.2992 

1C  1  . 2CC8 

7.00 

1 1 3. 70C0 

14.2146 

127.9146 

99.4854 

7.50 

1  1  4  .  1  0  C  0 

1  5.  86  5C 

129.9650 

98. 235C 

8.00 

1 14 .50C0 

16. C56  2 

1 20.5562 

98.4438 

8.50 

1 14.90C0 

16.  1184 

131.0184 

98.7816 

9.00 

1 15. IOC C 

17. C09  6 

1  32.  1096 

98. C9C4 

9.50 

1 1 5. 50C0 

18.  5C3 7 

134.0038 

96.9962 

10.00 

1  15.70CC 

22.  1589 

127.8589 

92.5411 

10.50 

116. 10C0 

32. 4139 

148.5139 

82. 6861 

T9 


TAPLE  HE -3 


VELOCITY  PROFILE 

MERIDIONAL  STREAMLINE  NO.  9 


M 

W 

DW 

W  SLCTICN 

W  PRESSURE 

.00 

144. 10CC 

1  5.  1797 

159.2797 

128.9203 

.50 

136.20C0 

1  5.  1771 

15  1  .  377  1 

121 . 0229 

1.00 

1 29. 10C0 

14.7211 

143.8211 

1  14. 2789 

1 .50 

123. 10C0 

1  3.  4406 

126.5406 

1C9. 6594 

2.00 

1 1 9. 9CCC 

12. C  86  C 

131 .9860 

1  C  7  .  8  1  4  C 

2.50 

1 17.20C0 

1C.  1354 

127.3354 

1C7.C646 

3.00 

1 15. COCO 

1  1.8577 

126.8577 

1  C3  .  1422 

3.50 

1 13. 50C0 

12. 9722 

126.4722 

ICC. 5278 

4.00 

112. 50C0 

10. 2646 

122.8646 

1  02 .  1354 

4.50 

112. 60C0 

8.  9383 

121 .5383 

103.6617 

5.00 

1 14.0000 

6.  9374 

12C.9374 

1C7. 0626 

5.50 

1 16. 50 CO 

6.  5427 

123.0427 

109.9573 

6.00 

1 1 8. CCCC 

8. 5  C94 

126. 5094 

1 C9 .4906 

6.50 

120. 2000 

1  0.  8457 

131. 0457 

1C8.2543 

7.00 

121.40CC 

13.  284  5 

1 34.6845 

1C8.  1  155 

7.50 

121  . 90  C  C 

14. 6543 

126.5543 

1C7. 2457 

8.00 

122. CCCC 

1 5. 532C 

127.5320 

106.4680 

8.50 

122.  10CC 

1  7.  4360 

1 39.5360 

104.6640 

9.00 

122.60CC 

17. 8984 

1 40. 4984 

1C4.7016 

9. 50 

123.20C0 

16.  2431 

139.4431 

1C6. 9569 

10.00 

124 .20CC 

20.  4772 

144.6772 

1 03 . 7227 
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Fig.  a 

De  5  in  Hi  on  oi  Armies 
on  a 

Meridtonal  Plane 


8a 


Fig.  3 

Velocity  Triangle 
•  on  +he 

TangenT  Planeo-Sa  Stream  Surface. 


Blade  Surface 


8  3 


Wm  =  Vm 


Location  StreamUnes 


Characteristics  Nos.  1^2. 


.Starting  Normals 


5 

Local' ion  of  a  Point. 
P* 

on  a  Characteristic 

-<Vv*N 
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Curvature  -  /inch 
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Fig. 8 

De+evmindi-iori  of  Iincremer>+s 
•for 

Pa  r<a  bo\  ic  xnterpo  la*h’on 


K1S'  Starlit^  formal 


*  **7*  .  ’■  y,  j «?  '  T*- ~* -  '*'*  -  -*  4*  :  •.  m*?r  v  • 


°u 


i/. 


r; 


AT  Tv 


92. 


o 


rrr 


*  > 


1  Jr* 


* 


Distance  along  Stream  1  me -inches 


>1 1. 
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APPSNDIX 

COMPUTER  PROGRAMS 
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Name 

ALFA 

AALFA 

DELN 

DFL 

DTDM 

F 

MO 

OMEGA 

RL 

THETA 

ATHSTA 

VM 

Additions 

BETA 

ABETA 


PROGRAMS  LEDGE  1  AND  LEDGE  2 

VARIABLE  NAMES 

Equivalent  to  Units 


c* 

deg 

rad 

AN 

in 

dF(L)/dL 

d  &  /dm 

rad. /in 

F(L) 

mt 

GO 

rad. /sec 

rl 

in 

© 

deg 

© 

rad 

Vm 

ft. /sec 

for  LEDGE  2 

@  deg. 

rad. 


08 


P roa ram  LEDGE  / 


PROGRAM  LEDGE  1 

DIMENSION  RL(9).ALFA(9),0ELN(9)  ,F <9 ) , EL ( 9 ) , VM  (  9  )  , 
1THETA(9,3),  AALF.AI  91.ATHETAI9) 

READ  1C.M0,RPM 

10  FORMAT  I  I  10®  F 10*0) 

READ  11, RL 

11  FORMAT I 6F 1 0 . 0  ) 

READ  1  1 , ALF A 
READ  ll.LELN 
Rt  AC  11, VM 

READ  1  1  ,  (  ( THE  T  A(M ,K ) , K* 1 , 3 ) , M* 1 , MO ) 

CO  12  M=1,M0 

AALFA(M)  =  ALFAIM 1*3. 14159/180. 

DO  12  K=1 ,3 

ATHET  A ( M»  K )  =  THE T A ( M, K ) *3 . 14 159/1  80. 

12  CONTINUE 

OMEGA  =  RPM*3. 141 59/30. 

DO  13  Ms 1 »  MO 

LT  DM  =  <-3.*ATHET AIM, 1 ) +4 . *ATHET A  I M,2 ) -ATHETA I M, 3)  1 
1 / ( 2. *CELN( M  )  ) 

F( M)  = ( VM ( M  )  • 1 2 •* CTCM+OMEGA )»RL(M)*RL(M) 

13  CONTINUE 
ELI  1 )  =  0. 
no  14  M=2,M0 

EL(M)=EL(M-1 )*(RL (M1-RLIM-1 ) ) »S INF ( AALF A ( H- 1 ) ) 

14  CONTINUE 
PRINT  21. 

21  I URMAT I / / / 1 CX,2 1HLEACING  EDGE  FUNCT ION/ / / 11 X 1 HM 1 IX 3HDFL ) 
DO  22  M*1,V0 

IFIM-1)  15,15,16 

15  MA  =  1 
MB  =  2 
MC  =  3 
GO  TO  19 

16  IF  I M- MO)  17,18,18 

17  MA  =  M 
MB  *  PHI 
MC  =  M-l 
GO  TO  19 

18  MA  =  MO 
MB  =  MC- 1 
MC  =  MC— 2 

19  XA  »  CL (MB)-ELIMA ) 

XB  =  ELIMO-ELIMA  ) 

DFL  =  (XB»XB*(F(M8)-F(MA) )-XA«XA*(F(MC)-F(MA) )) 

1 / l ( XA*XB«XB-XA»XA*X8)*12.*SINF( AALF A ( M) )) 

PRINT  20, M, DFL 

20  FORMAT!/ 10X, I  2 , 5X ,F15.5) 

22  CONTINUE 
END 

END 


I0O 


Program  LEDGE  2 


© 


lov 
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PROGRAM  LECGE  2 

DIMENSION  RL(9),ALFA(9J,F(9),EL(9),VH(9), 

1  BE T A ( 9) , A BETA (9),  AALFA(9) 

READ  10,MC,RPM 
10  FORMAT ( I  1  C  ,  F  1  C  .  0) 

READ  11, RL 
U  FORMAT! 6r 10. C) 

READ  11, ALFA 
READ  11, VM 

READ  22, BETA! 1 ) , B  ET  A ( 9 ) 

32  FORMAT (2F 10.0 ) 

ABETA(l)  =  BETA!  1  )  *3.14159/180. 

ABE  T  A ( 9 )  »  RETA(9)*3.  1 4 1 59/ 1 8C . 

DO  12  M= 1 , MC 

AM  =  FL0ATF(M-1 ) / FLG ATF ( MO- 1 ) 

AALFA(M)  =  ALFA(M) *3.  14159/18C. 

12  APETA(M)  =  ABETA!  1  )  +  { APETA ( 9 ) -ABET A ( 1  ) )*AM 
OMEGA  =  RPM *3. 141 59/30. 

DO  91  M= 1 , MC 

G £ T A ( M )  =  APETA(M  )<  18C./3. 14159 
91  CONTINUE 

00  13  M= 1 , MO 

F(M)  =  VM  (  M  )  *  P.  L  ( M  )  *  T  A  N  F  (  A  B  E  T  A  (  M )  )  *  1  2  .  +C  M  EG  A  *  R  L  !  M )  *•  R  L  (  M ) 

13  CONTINUE 

EL!  1  )  =  0.  * 

DO  14  M»2,M0 

EL! M) =EL( M-l ) + ( RL ( M)-RL(M-1 ) ) *S INF! A ALFA (M-l ) ) 

14  CONTINUE 
PRINT  51 

31  FORMAT! 1H1) 

PRINT  21 

21  FORMAT  (///1CX,  2  1H  LEADING  EDGE  FUNCTION’///  11  X  1  HM  11X  3HCFL  ) 
CO  22  M= 1 ,MQ 

I F ( M- 1 )  15,15,16 

15  HA  =  1 
MB  =  2 
MC  =  3 
GO  TO  19 

16  IF(M-MO)  17,18,18 

17  MA  =  M 
MB  =  M+ 1 
MC  =  M-  1 
CO  TO  19 

18  MA  =  MO 
MB  =  MC- 1 
MC  =  MO-2 

19  XA  =  EL ( MB ) -EL ( MA ) 

XB  =  EL  (  MO-EL  (MA  ) 

CFL  =  (XB*X8*(F(HB)-F(MA)  )-XA *X A* ( F ( MC )- F ( MA )  ) ) 

1/( (XA*XB«XB-XA*XA*X3 1*12. *S INF! AALFA(M) ) ) 

PRINT  2 C , M , C F L 
2C  FORMAT!  /  10X,  12.,  5X  ,  FI  5.5) 

22  CONTINUE 
END 

END 
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PROGRAMS  ROTOR  1  AND  ROTOR  2 


Name 


VARIABLE  NAMES 

Equivalent  to  Units 


BETA 

PBETA 

T3ETA 

XBETA 

2BETA 

BLNO 

CURV 

PCURV 

D1 

D2 

D 2 TDM 2 

D2TDNM 

DELN 

DELTA 

ADELTA 

DDELTA 

PDSLTA 

YDELTA 

DHEDN 

DM 

DNS TAR 

DTDM,  DTI,  DT2 
DWCO 


(bat  P* 

£  at  P* 
tan  (^> 
at  (P*+  Sx) 
(b  at  (P*+  £x) 
N 


k  at  P* 
m 

Di 

°2 

d2©  /dm2 
d2©/dndm 
AN 

5 

5 


S  at  P* 

Sat  P* 

£  at  (P*+  5x) 


dM 

dN* 

d  ©  /dm 
DWcoef. 


deg. 

rad. 

rad. 

deg. 


in. 

deg. 

rad, 

deg. 

rad. 

rad. 

in. 

in. 


104* 


Name 

Equivalent  to 

Units 

DWFUNC 

DW_ 

func 

DX 

dx 

in. 

SN 

5N/2 

in. 

FX 

3(x) 

GAMMA 

'Kat  P* 

deg. 

GAMX 

If  at  (P*+  5x) 

deg. 

PGAMMA 

Yat  P* 

rad. 

TGAMMA 

tan  ^ 

XGAMMA 

^  at  (P*+  Jx) 

rad. 

AL  AM 

A 

rad. 

DLAM 

A 

deg. 

PLAM 

A  at  P* 

rad. 

XLAM 

A  at  (P*+  5 x) 

rad. 

OMEGA 

GO 

per  sec. 

QA 

Q  (actual) 

cu.  ft. /sec. 

q: 

Q  (calculated) 

cu.  ft. /sec. 

QDEL 

AQ 

cu.  ft. /sec. 

DQ 

^E(x)dx 

cu.  ft. /sec. 

R 

R 

in. 

AR 

R  at  P* 

in. 

RX 

R  at  (P*+  5x) 

in. 

TKETA 

© 

deg. 

ATHETA 

e 

rad. 

THICK 

At 

in. 

TKEQV 

At' 

in. 
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Name 


Equivalent  to 


Units 


W 

WDEL 

m 

wx 

WSQ 
DWSQ 
X,  AX 
XDEL 

Y1 ,  Y2,  Y3 

Z 


N,  NN,  etc. 

M,  MM,  etc. 

MO 

NO 

NS 


W 

dW 

at  P* 
m 

Wm  at  (P*+  5x) 

W2  at  P* 
m 

dW2/dx 

m 

x 

5x 


ft. /sec. 
ft. /sec. 
ft. /sec. 
ft. /sec. 
ft. 2/sec. 2 
ft. 2/sec. 2-in. 

in. 

in. 

in. 


INDEX  NAMES 

Normal  number 
Streamline  number 
Total  number  of  streamlines 
Total  number  of  normals 
Starting  normal 


Changes  and  additions  for  ROTOR  2 


ABETA 

APBETA 

DTBN 

DTBM 


(3 

(b  at  P* 
d(  tan  )/dn 
d(  tan/^  )/dm 


rad. 

deg. 
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Program  ROTOR  / 


Main  Program 


|£WP| 


Ygs 

i 

© 


S.R.  COY>  ^ 


No 


*E lements  of  Y9  s 

Calculat 

e  Y9s 

Yes 


Wn*  -  0 


dW  =•  iOO 


W  =  W  +  oV/ 
m 


S.tf.  RELVEL-* - *- 


Velocity 

Distribution 


r-(-) 


i  Aq|  -  ^ 
0) 


(+)- 


I nc remen t 
cW 


& 
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S.R.  FLO W  ~K-- 


i  10 


PROGRAM  RCTCR  1 

DIMENSION  CLAM(9,  10)  ,R(  9,  1C)  ,Z(9. 10)  ,  THETA(9, 10)  , 

1ALAMC9,  1C) »  ATHETA ( 9  » 1 C ) , DEL T A ( 9 , 1 C ) , ADEL TA ( 9 , 1 C ) , 

2 AR ( 9  )  ,  PL  AM  (  9  )  ,  GAMMA  (  9  )  , PGAMMA  (  9  )  ,  PBET  A  (  9  )  ,  CELM  9  )  , 
3PDELTA(9),D.N$TAR( 9) ,  DM.  (  9  ) »  D 1  ( 9) »D2(9) ,CURV(9, 1C) , 
4PCURV(9),Y1(9),Y2(9),Y3(9) , DHEDN ( 9 ) ,CFL(9),kM(9), 

5DWSQ(  9)  ,I,X(  9)  ,  AX(  9)  ,FX(9)  ,PETA(9)  ,00ELTA(9)  ,XLAM(9  )  , 

6RX ( 9 )  ,WX(  9)  ,  X  GAMMA  (9)  , XBETA( 9 )  ,  GAMX  ( 9  )  ,W(9),XCEL(9), 

7ZBETA(9)  ,GML(9)  ,YCELTA(9)  , CWCC(9)  ,D’aFUNC(9)  ,  TKECV( 9) 

COMMON  M, DELN, ATHETA, NO, D TDM , PP E T A , PD  EL TA , PGAMM A , P l AM , 

1 AR, R, CM,DNSTAR, ADELTAtNOCO, ALAMfNS, 01 , D2 , CURV , PCUR V , MO , 

2Y1 , Y2  ,Y3,  t«M,QOEL* DWSC , OMEGA, CX, AX ,QA,  QC, 

3WX,RX,XGAMMA,XBET A, YCELTA 
READ  1C, MC, NO, NS 

11  FORMA  TC6F10.0) 

1C  FORMAT (31 10) 

READ  11, ( DELN ( M ) , M= 1 , MC ) 

RE  AC  11  ,  (  (  R  (  M ,  N  ) ,  N  =  1  ,  NO  )  ,  M=  1 ,  M C  ) 

READ  11 , (  (Z ( M, N),  N=1 , NO ) , M  =  1 , MC) 

READ  11 , ( (DLAMtM, N) , N= 1 , NO ) , M= 1 , MG ) 

READ  1  1  ,  (  (THETA  (M  ,  N  )  ,  N=  1  ,  NO  )  ,  M=  1  ,  M.O  ) 

READ  11  , ( (OELTA(M, N ) ,N  =  1 ,N0 ) , M* 1 ,M0) 

READ  11,  (CURV( 1,N) ,N=1,NC) 

READ  11,  (  CL'RV  (  9,  N  )  ,  N=  1 ,  NC  ) 

READ  90  ,  OMEGA ,  L'A  s 

90  FORMAT ( 2F 15.0) 

READ  11,  (L’HEDN(M)  ,M=1,M0) 

READ  11, (OFL(M) ,M=1,MC) 

RE AC  11, THICK, BLNC 
DO  12  M  =  1 , M 0 
CO  12  N  =  1 , N 0 
C  =  3.14159/180. 

AL  AM ( M » N )  =  D  L  A  M ( M , N ) #C 

A  THET  A ( M, N )  =  T HE  T  A ( M , N ) *  C 
ADELTA  (MfN)  «  DE  LT  A ( M, N ) *C 

12  CONTINUE 
M  =  1 

CALL  CTD(NS) 

AR  (  M  )  =  R  (  M  ,  N  S  ) 

POULT  A ( K )  =  AOELT  A ( M, NS) 

CALL  ANGLE 

GAMMA(M)  =  PGAMMA  (  N\ )  /  C 
DNS  TAR ( 1 )  =0. 

DM ( 1 )  =  C. 

PLAM(M)  =  A  L  A  f  ‘  ( M , N  S ) 

CETA(M)  =  P  B  E  T  A ( M ) / C 
Di)ELTA(M)  =  PDE  LT  A  (  M  )  /  C 
X L A M  (  M  )  =  P  L  A M  (  M )  /C 
PRINT  929 
929  FORMAT ( 1 H 1 ) 

PRINT  13 

13  FORM ATI 1GX39HLCCATI0N  OF  POINTS  ON  CHARACTERISTIC  ///) 

PRINT  1l< 

14  FORMAT ( 2H  MOX 5HGA MM A  1 CX6H0NS T AR9X4HBET A  1 1 X5FDELT Al 0X5HLAMDA 

1  1  2X1  HR 14X2HDM//) 

PRINT  15, M,GAKMA(  M)  , U NST A R ( M ) , RET A ( M ) , CC EL T A ( M ) , XL  AM ( M ) , A R (M ) , 
1  DM  (Mi) 

15  FORMA T(  1 2, 7F  15.4//) 

DO  16  M  =  2 , M C 

NOGO  =  C 

NH  >»  f.NSl  AR  (  M-U/OELN  (  M- 1  ) 

N't  =  NS+NP 

I  F ( PL  AM ( M-  )  ) - .  785  4 )  17,  17,  18 

17  ADM  =  (R(M,NE)-R( M-1,NE) )/CCSF(PLAM(M-1) ) 

DM( M)  =  ABSF(ADM) 

GO  TO  19 

18  ADM  =  ( Z ( M  »  NE ) -Z( M- 1 , NE )  ) /SI N F ( PL AM ( M- 1 ) ) 

OM(M)  ~  APSF(APV) 

19  l)NSTAP(M)  =  DNSTAR  ( M-  1  )  +  DM  (  M  )  *T  ANF  (  PG  AMM  A  (  M-  1  )  ) 

IF(DNSTA'MM)  )  20,21,22 

20  CALL  FSTAR(NS,-1,-1. ) 

GO  TO  23 

22  CALL  P  S  T  A  R  (  N  C  -  N  S »  1,1.) 


Ill 


21 


24 

25 
23 

26 
27 

2G 


16 

951 

952 


970 


971 

31 

97? 


32 


33 

34 

35 


GO  TO  23 
CALL  CTO(NS) 

PLAM ( M)  *  ALAM(M.NS)  t 
POELT A ( Ml  *  ADELT  A(M,kS ) 

AR(M)  =  R(M,NS) 

CALL  ANGLE 

I F ( PGAMM A ( H ) )  24,23,25 
CALL  PSTAR(NS,-1,-1. ) 

GO  TO  23 

CALL  PSTAR(NC-NS  ,1,1.) 

CONTINUE 

IF(I-NCGO)  26,26,  26 
PRINT  27, M 

FORMAT ( 27HINSUFFI Cl  ENT  GRID  WIDTH  M«,I2) 

GO  TO  55 

GAMMA  ( M  )  =*  PGAMMA  !M)/C 
BETA(M)  =  PBETAIM ) /C 
COELTA(M)  =  PDELT  A ( M ) /C 
XL AM( Ml  =  PLAM<M)/C 

PRINT  1 5, M, GAMMA!  M) , DNSTAR(M) , BET AIM) ,OOELTA( H ) , XL AM ( M ) , AR ( M) , 

1  DM ( M  ) 

CONTINUE 
PRINT  951' 

FORMAT! 1H1, 20X1 2H TEST  FOR  COY//) 

PRINT  952 

FORMAT (2H  M3X2HNX8X6H02TCM29X6BD2TDNM 11X2 HO  1 1 3X2H02 12X4HCUPV// ) 
CALL  COY 
CONTINUE 
PRINT  97C 

FORMAT! 1H1 ,5X,6HY  TEST//) 

DO  31  M*1 ,M0 

Y1IM)  =  COSF(PGAMMA(M))*(2.»PCURV(M)*(COSF(PB£TA(M)))««2 
1+SINFI2.* PBETAIM) )*C2(M))«12. 

Y2 ( M I  =  2.*C0SF IPGAMMAIM ) )*SINF(2.*P8ETA(M))«CUM) 

1/COSF (PDELTAIM) ) 

Y3!M)  =  I OHECNIM ) -OMEGA* OF LI M) )*COSF IPLAMIM) )*2.» 

1C0SF! PCAMMA(M) )•(  CC SF ( PR ETA ( M ) ) ) «*2 
PRINT  971, M. Y1!M) .Y2(M),Y3(M) 

FORMAT! 110, 3F  15.5 ) 

CONTINUE 
PRINT  972 

FORMAT! 1H1,5X,8HVEL  TEST//) 

EE  = .  .001 
WMl 1 )  =0. 

1  =  1 

WCEL  =  100.»l  —  .  1«  «C  I  —  1) ) 

MM! 1 )  =  WMl  1  UWOEL 
CALL  RELVEL 

IF!  AH  SF  ( (.DEL  )-EE  )  38,38,33 
IFICDEL)  32,38,34 
I  =  I  •»  1 

WDEL  =  ICO. *!“•!•• I  I- 1 ) ) 

WM  (  1  )  =  NM(  1  )+WDEL 
CALL  RELVEL 

IF(ABSF(COEL)-EE)  38,38,36 
IF ( QUEL )  37,38,35 


13HTEST  FOR  XOEL//) 


36 

37  I  =  !♦! 

GO  TO  32 

38  CONTINUE 
PRINT  9ei 

981  FORMAT ( 1 H 1 
PRINT  982 
98?  FORMA  f ( 3H  M,3H  J,15H 
115H  DELTA  X  //) 
983  FORMAT ( 36 X, F 1 5. 4) 

WX(  1  )  =  WM( 1) 

WX(MO)  =  WMIMO) 

RX! 1 )  =  AR!  1) 

RXtMOJ  =  AR (MO ) 

XPETA!  1  )  =  PR  ETA!  1  ) 
XBETA (MO)  =  PBETA (MO) 
XGAMMA(I)  =  PGAMMA! 1) 
XGAMMA(MC)  =  PGAMMA(MC) 


FLOW  FRACTION  ,15H  CAL  FLOW  RATE 


PDELTA!  T  ) 

*  PDELTA(MO) 


YDELTA(l)  = 

YDELT  A ( MO  )  > 

MO  1  =MO  1 
00  950  K=2,M01 
KK  =  NO+T-K 

F  =  FLCATF(M0-K)/FL0ATF(M0-1) 

OAF  =  CA*  F 
EE  =  .001 
CALL  FLOW ( KK »  0. ) 

IF(A8SF(QC-QAF)-EE)  50,50,40 
XOEL ( KK  )  =  .2 
PRINT  985  » XCEL ( KK ) 

ODD  -  XOEL ( KK ) 

CALL  FLOW  IKK, ODD) 

IF(GC-CAF)  142,50,43 
XDEL(KK)  =  XCEL (K  K ) ♦. 1 
ODD  =  XCEL(KK) 

PRINT  98  3, XOEL ( KK  ) 

C...  DRIVES  CC  GREATER  THAN  QA  TO  START 
GO  TO  41 
J  =  1 

XOEL ( KK )  =  XDEL(KK)-. 1*  ( -  • 1  ** ( J- 1 ) ) 
000  =  XOEL(KK) 

PRINT  985,XCEL(KK) 

CALL  FLCW ( KK , 000) 

I F ( ABSF ( QC-CAF )-E  E  )  50,50,45 
IF(QC-CAF)  46,50,44 
J  =  J+1 

XDEL(KK)  =  XDEL(KK)-. ) ) 
ODD=XCEL ( KK ) 

PRINT  98J,XCEL(KK) 

CALL  F LCW ( KK , 000 ) 

IF (ABSF(GC-QAF)-EE)  50,50,48 
1 F ( GC-CAF  )  47,50,  49 
J  =  J+1 
GO  TO  44 
CONTINUE 


40 


41 

42 


43 

44 


45 

46 

47 


48 

49 

50 

980 

950 

984 

985 


PRINT-  98Q.KK,  J,F,QC,XDEL(KK) 
FORMAT! 21 3,3F15.4  ) 


CONTINUE 
PRINT  984 

FORMAT!  1H1,  16HVEL0CITY  PROFILE//) 

PRINT  935 

FORMAT ( 2H  M8X6HRA C I  US 4X 1 CHRL ACE  ANGLE8X5HGAMMA9X2HWM 
1 15H  REL  VELOCITY  W  // ) 

00  51  M- 1 , MC 
ZRETA(M)  =  XPETA(M)/C 
GAMX(M)  =  PGA KM A{ M ) /C 

WtM)  *  V*X(M)*SCRTF( 1 . 4 ( T ANF ( XPE T A ( M ) )**2) ) 

PRINT  986,M,RX(M) ,EETA(M) , G A M X { M ) ,WX(M)  ,W(M) 

986  FORMAT! I2,5F15.4) 

TKEGV(M)  = ! THICK/COSF ! XBETA(M) ) JMSQRTF! 1. -MS  INF 
1 (XPETA(M) )*TANF!YCELTA(M) ) )**2) ) 

DWCO(M)  *  (6.2832/0LNC-TKECV!M) /RX(M) )*C0SF (XBETA! M) ) 
OWFUNCIM)  =  RX(M) »UX(K)»TANF(XBETA(M) )*12.+GNEGA*RX(M) 
1  *  R  X  (  M  ) 

CONTINUE 
PRINT  880 

FORMAT! 1H1/////// /25X9H  TABLE  ///17X 
127H  DATA  FOR  CHARACTERISTIC  Cl///) 

CHANGE  CHAR.  NO.  FOR  EACH  COMPUTATION 


51 

880 

382 

883 

884 

885 
897 

886 


FORMAT ( 6H  M  =  8X2H  WlH2  9X1H3  9X1H4  9X1H5//) 
PRINT  883 

FORMAT! 13X33H  LOCATION  OF  CHARACTERISTIC  CURVE//) 
PRINT  834,  NS 

FORMA  T ( 23H  STARTING  NORMAL  NO.  =  ,12//) 

PRINT  88? 

PRINT  885, (GAMMA! M ) , M  =  1 .5) 

FORMAT! 10H  GAMMA  =,5F1C.4) 

PRINT  897, (GML(M) ,M«1 ,5) 

FORMAT!  1CH  GAM-LAM  =,5FTC.4) 

PRINT  886, !CNSTAR!M),M=1 ,5) 

FORMAT!  1CH  ON ST AR  =,5F1C.4) 

M3 


AX 


PRINT  887,(BETA(M)fM=1,5) 

887  FORMAT! 1CH  BETA(P)  =,5F10.4) 

PRINT  P80,(AR(M), M= 1,5) 

88P  FORMAT (  1 CH  RACIUS  =,5F1C.4///) 

PRINT  809- 

889  FORMAT! 16X28H  LOCATION  OF  NEW  STREAMLINES//) 
PRINT  882 

XDEL!  1  )=C. 

XDEL(MC)=0. 

PRINT  890, ! XDEL(M  ),M=1 ,5) 

890  FORMAT  (  10H  DELTA  X  =,5F10.4) 

PRINT  891  .(ZBETA(M),M*1't5) 

891  FORMAT (  1CH  BETA(X)  =,5Fl0.4) 

PRINT  88Q, (RX!M),M=1,5) 

PRINT  892 

892  FORMA T ( 2 1 X 1 7H  VELOCITY  PROFILE//) 

PRINT  882 

PRINT  8  9  3  ,  (  W  M  (  M  ) »  M  =  1 , 5 ) 

893  FORMAT! 1CH  WM(P)  =,5F1C.4) 

PRINT  894  , ( MX ( M ) ,  M=  1 , 5) 

894  FORMAT l 1 CH  WM ( X )  =,5F1C.4) 

PRINT  895,(W()A),M  =  1,5) 

895  FORMAT ( 1 CH  REL  VE L  =,5F1C.4) 

PRINT  875,  ICWCOIM  )  ,V=  1 ,5  ) 

875  FORMAT  (  1CH  OW.COEF  =,5F10.4) 

PRINT  874, (CWFUNC !M) ,M=1 ,5) 

874  FORMAT ( 1 CH  CW  FUNC  =,5F1C.4) 

896  FORMAT ( 6F  M  •=  flX 2H  69X1H7  9X1H8  9X1H9//) 
PRINT  88C 

PRINT  883 

PRINT  884, NS  " 

PRINT.  e96 

PRINT  885, (GAMMA! M),M=6.9) 

PRINT  897, (GML(M) ,M=6,9| 

PRINT  886, (DNSTAR(M) ,M=6.9) 

PRINT  887,(BETA(M ) ,  M  =  6 ,  9  ) 

PRINT  876,(AR(M), M=6 *  9 ) 

878  FORMAT ( 10H  RADIUS  =,4F10.4///) 

’ PRINT  889 
PRINT  896 

PRINT  89C,!XDEL(M),M=fc,9) 

PRINT  091, (/BETA! M),M=6,9) 

PRINT  873 , (RX(M), M=6,9) 

PRINT  892 
PRINT  896 

PRINT  893, (WM(M)»  M*6  »  9 ) 

PRINT  894, (WX(M),M=6,9) 

PRINT  895,  (  W  (  M  )-,M  ~6 , 9  ) 

PRINT  875, (CWCOIM ),M=6,9) 

PRINT  874, (DWFUNC(M) ,M=6,9) 

PRINT  860 

860  FORMAT! 1H  1 , 4H  END) 

55  CONTINUE 
END 
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Program  ROTOR  I 


Subroutine  PSTAP 


.  SUBROUTINE  PS  TAR<  NN ,  JO  ,  R  ) 

01  MENS  ION  DLAMI9,  10)  .R(9,  10)  ,2(9,  101  ,  THETAI9,  10)  , 

1 ALAM(5, 1C) . ATHETA (9, 10 ), CELIA (9. 10) ,A0ELTA<9, 10)  . 

2AR(9) ,PLAM(9) ,  GAM  M A < 9 ) , PG AMM A ( 9 ) , PBE T A ( 9 ) ,DELN(9), 
3PDELTAI9) ,DNSTAR( 9).CM(9).D1(9),D2(9),CURV(9,10», 
4PCURV(9),Yl(9).Y2 (9) ,  Y 3  (  9 ) , OH  EON < 9) ,DFL(9),WM(9), 

50WSG ( 9  ) »  D  X ( 9 )  ,AX(9),FX(9),BETA(9),DDELTA(9),XLAM(9), 

6RX ( 9  )  ,  WX  (  9  ) , X  G  A  MM  A ( 9 ) ,  XBE T A ( 9 ) , GAMX ( 9 ) , W ( 9 > , X C EL ( 9 ) . 
7ZBETA  (9)  ,GML(9»  ,YCELTA(9)  ,0WC0(9)  .DWFUNC  (  9  )  .T  KECV  (  9  ) 
COMMON  M,CELN, ATH ET A , NO , 0 TOM , PBETA, PDEL TA, PGAMMA ,PLAM , 

1 AR,R,CM,DNSTAR, AD  EL T A ,NOGO , AL AM , NS , 0 1 ,  D2  ,  CURV  , PCLR V , MC , 
2Y1 , Y2,Y3,UM,QOEL, DWSC , OMEGA, 0 X, AX ,QA , UC , 

3WX, RX .XGAMMA, XBET  A , YOELTA 
CO  30C  1=1,2 
CEL  =  OELN(M) 

00  30  1  K=  1  , NN 
NA=NS ♦  JO* ( K- 1 ) 

NB  =  NS* JG*K 

X  =  (CNSTAR(M)-B*CEL)/OELN(M) 

I F ( ABSFICNS TARIM)  ) -CEL ) 302, 30  3,304 

304  DEL  =  DELN(M)«FL0ATF(K*1 ) 

GO  TO  3C1 

302  CALL  CTO(NR) 

OT 1 =D  TOM 
CALL  CTD(NA) 

0T2  =  CTCM 

CTDM  =  CT1-BMDT2-CT1 )*X 

AR ( M )  =  R(M,NB)-B«(R(M,NA) -R (M,NB))*X 

POELTA(M)  =  ADELT  A( M , NB ) (A0ELTA(M,NA)-A0ELTA(M, NB) )*X 
GO  TO  305 

303  CALL  CTO ( NB ) 

AR(M)  =  R ( M , NB ) 

POELTA(M)  =  ADELT  A ( M , NB ) 

305  CALL  ANGLE 

PLAM(M)  =  ALAM(M,NR)-B*( ALAM(M,  NA ) -A L AM ( M, NB ) ) *X 
I F ( I-  1  )  306,306,300 

306  TG AMD  =  T ANF ( PG AM M A ( M ) ) - T ANF ( PG AMM A ( M- 1 ) ) 

IF(TGAMD)  3C7.300.3C7 

307  EN  =  (DM(M)»TGAMU )/2. 

DNSTAR(M)  =CNS  T  AR ( M  >  +  EN 
GO  TO  300 

301  CONTINUE 
NOGO  =  I 
300  CONTINUE 
END 


1  i  (o 


SUBROUTINE  DTD(N) 

DIMENSION  DLAM(9, 10  )  .  R ( 9 , 10 ) , Z ( 9. 1 0  I , THET A ( 9, 10)  , 

1 ALAM(9, 10  ) , ATHETA (9, 10) , DELTA (9, 10),A0ELTA(9,10>. 

2AR ( 9 ) » PL AM ( 9 ) , GAM K A ( 9  )  , PGAMMA ( 9 ) , PBET  A ( 9 ) ,CELN(9), 
3PDELT  A  {  9 ) , ON STAR (  9 ) . UM ( 9 ) .D 1 ( 9 ) , U2 ( 9 )  ,CURV(9,  1C), 

4 PC UR V (  9  )  ,  Y  1  ( 9  ) , Y2 (9) ,Y3( 9) ,DHEDN(9) , OFL ( 9 ) ,WM(9), 
5DWSCM9) ,DX( 9) , AX( 9) ,FX(9) , RET  At  9) ,  DDE  L  T  A  (  9  ) ,  X  L  AM (9 ) , 

6RX ( 9) ,KX( 9) ,XGAMM  A ( 9 ) , XB E T A ( 9 ) ,GAMX( 9 ) ,U (9) ,XDEL(9  ). 
7ZBET A (9 ) ,GML(9) ,YDELIA(9) , DWCC ( 9 ) , DWFUNC ( 9 ) , T KEC V(  9) 
COMMON  M.fJELN,  ATM  ETA,  NO.DTDM,  PBETA.PDLLTA,  PGAMMA.PLAM, 

1 AR , R , CM , CNST  AR, AD  ELT  A ,NOGO, AL AM, NS, D1 ,D2 ,CURV , PCLR V , MC , 
2 Y 1 , Y2,Y2,WM,OOEL, OWSC  ,  CM E GA  ,  D X  ,  AX ,  CA  ,  C.C  , 

3WX, RX  ,XGAMMA, XBET  A, YU EL  TA 
DELTA  =  1  2  .  «DELN(  M  ) 

I  F  (  N—  2  I  ICO.ICO.ICI 

100  OTDM  =  (-25. *  ATHE  TA(M,N)+4fl.*ATHETA(M,N+1 ) -36 . • 

1 ATHETA (M.N+21+ 16. "ATHETA ( M, N+ 3 ) -3 .  * A THE T A ( M, N  +  4 ) ) / DE LT A 
GO  TO  104 

101  NREM  =  NC-N 
IF(NREM-I)  102,102,103 

102  CTOM  =  (3.*ATHETA ( M, N-4 ) -  1 6 . * ATHET A ( M , N-3 > ♦ 36 . * 

1 ATHETA( M, N-2J-48. * ATH£T A ( M, N- 1 ) +2 5 . * ATHE T A ( M , N ) )  / OELT  A 

GO  TO  1C4 

103  OTOM  =  (ATHETA( M, N-2)-ATHETA( M,N*2)  -8 . * ( ATHE T A ( M, N- 1 1 
1-ATHETA(M,N+1 ) ) )/ DELTA 

104  CONTINUE 
END 


1 1  7 


SUBROUTINE  ANGLE 

DIMENSION  DLAM(9,  1C) ,R(9,  1C) ,Z(9. 10) , THETAI9, 10)  , 
1ALAMI 9, 1C) , ATHETA (9, 1C) ,DELTA(9, lO) ,ADELTA(9,  10), 

2 AR ( 9 ) ,PLAM(9),GAMMA(9),PGAMMA(9) , PBET A ( 9 ) , CELN ( 9 ) , 
3P0ELTAI9)  ,ONSTAR(  9  )  .  DM  (  9  )  .  0  1  ( 9  )  , 02  (  9  )  , CURV  ( 9 , 10  )  , 
4PCURV  I  9  )  ,  Y  1  (  9  )  ,  Y2 I9),Y3(9),DHEDN(9),DFL(9) ,WM(9  > , 

5DWSC  (  9  )  ,  UX.(  9  )  ,  AX(  9  )  ,  F  X  (  9  )  ,  BE  T  A  (  9  )  ,  DDE  LT  A  (  9  )  ,  XL  AM  (9  )  , 
6RX ( 9 ) ,WX<  9),XGAMMA(9)  ,XBETA(9)  ,GAMX(9)  ,M9) ,XCEL(9 ) . 
7ZRETAI9)  ,GML(9) .YDELTAI 9) ,DWC0(9) ,DUFUNC(9) , TKEQVI  9) 
COMMON  M,DELN, ATH ETA, NO, DTDM.PPET A, POELTA, PGAMMA.PL AM, 
1 AR,R,CM,DNSTAR, AO ELTA ,NOGC, AL AM, NS, 01 , 02 , CUR V , PCUR V , MO 
2 Y 1 ,Y2,Y3, WM.GDEL, OWSC , OMEGA, CX , AX ,CA, QC , 

3UX.RX  .XGAMMA, XPET  A, YDELTA 
TBETA  =  AR ( M ) «OTD  M 
PBETA { M )  =  AT  ANF ( TPETA) 

TGAMMA=(SINF( PRET AIM) )»*2)»  T  ANF ( PDELT  A ( M  )  ) 

PGAMMA(M)  =  ATANF (TGAMMA) 

END 


Subroutine  COY 


S.R..  DTD 


o 


END 


Yes 


US 


401 


SUBROUTINE  COY 

DIMENSION  DLAMI 9,  10) tR(9, 10).Z(9, 10)  .THETAI9,  10)  , 

1  ALAM(  S»1C  )  »ATHETA  (9,  1C)»C£LTA(9.10),A0ELTA(9,  10)  . 

2 AR ( 0 )  .PLAMI9) ,GAMMA(9).PGAMMA(9),PBETA<9) ,CELN(9), 
3PDELTAI9) , DNSTARI 9 ) . DM( 9 ) . D 1 ( 9 ) ,02 < 9 ) , CURV ( 9 , 10) , 
4PCURVI9) ,  Y 1 (9) fY2 (9) , Y3 1 9  )  ,DHEDN 1 9 ) ,DFL I  9 ) .WMI9  ) . 

5DWSQI 9) ,DXI9).AXI  9),FX(9) ,BE T A I  9 ) , DOELTA ( 9 ) , X  LAM  19  ) , 
6RXI9)  ,WX(9) ,XGAMMAI9)  ,  XBE  TA  ( 9 ) ,  GAMX  ( 9  )  ,.W  19  )  ,XCEL  19  ) . 
7ZBETA(9),GML(9) »YDELTA19) , CWCO ( 9) ,  DfcFUNC I  9 ) , TKEQ VI  9 ) 
COMMON  M »  OELN  , ATH  E  T  A, NC, C TDM , PBE T A , PC EL T A , PGAMMA , PL AM, 
1AR,R,CM,CNSTAR,ADELTA,NOGO,ALAN,NStDl,D2»CL'RV»PCLiRV,MO, 
2Y1,Y2»Y3»WM,QDEL»  CVvSC , OMEGA,  DX,  AX,  QA,  CC  , 

3UX,RX ,XGAMMA,XBETA,YCELTA 
DO  4CC  M  =  1 , MG 
NX  =  CNSTAR(M)/DELN(M) 

IFICNSTAR (M) )  401,402,403 


402  NA 
NB 


NA 
NE 
XA 
XB 
JO 
R  = 

GO  TO  404 
NS+1 


1 


NS  +  NX 
NA-1 

-DNSTARI M)+FLOAT FI  NX )*DELN (M ) 
XA-CELN(M) 

-  1 


NS-1 

CELN(M) 

-XA 

-1 


XA 
XB 
JO 

•  R  =  1 

GO  TO  4Q4 

403  NB  =  NS+NX 
NA  =  A  B  + 1 

XB  =  -DN$TARIM)+FLOATFlNX)*DELNIM) 

XA  =  CELN ( M )— XB 
JO  =.  1 
n  =  -i. 

404  CALL  CTC(NA) 

CA  =  CTCM 
CALL  CTD(NB) 

OB  =  CTDM 

DC  =  TANFIPBETAIM  M/ARIM) 

C2TCM2  =  IXA+XA*!  C8-0C )-XB»XB» I OA-DC) )/ 

1 (XA*XA*XB-XB»xe»X A) 

NN  =  NS+NX 

IFIM0-2-M)  405,406,406 

405  ZA  =  CMIM) 

MZ  =  M 

M  =  M-1 

ZB  =  ZA-CMIM) 

CALL  CTDINN) 

ON A  =  CTCM 
M  **  M  *  1 

CALL  CTC(NN)  .  . 

DNB  =  DTDM 
M  =  M2 
GO  TO  407 

406  MZ  =  M 
M  =  M+  1 

ZA  =  CM ( M  ) 

CALL  CTCI NN) 

DNA  =  CTCM 

M  =  M+1 

ZB  -  ZA+OMIM) 

CALL  CTCINN) 

DNB  =  OTDM 
M  =  MZ 

407  CALL  DTDINN ) 

CNO  **  C'T D K 

C2TCNM  =  (  ZB  +  ZBM CNA- CNO ) -ZA*ZA« ( DNB- CNO  )  ) / 
1CZR«ZB»ZA-ZA‘*ZA*ZB) 

ALFA  =  PCELTAIM)-PLAM(M) 

01 (M)  =  CCSFI ALFA ) 

02 ( M )  =  (2.«TANF(PBETA(M))«D1(M))/(AR(M)*CCSF(PDELTA(M))) 


120 


1  +  A  R  (  y  )*(  r  ANF(  PCEITMM  )  )*C2TDM2  ♦  D2TDNM  } 

IF(M-I)  408,408,409 

408  PCURV ( 1 )  =  CURVM ,NS) 

954  GO  TO  950 

C — REPLACE  M  TH  (GOTO  4C0) 

409  AM  =  FLCATF(M-l)/FLCATF{MO-1) 

CURV  ( M,NN )  =CURV(  1  ,  NN )  +  { CUR V ( MO , NN ) -CUR V (  1,NN))*AM 
CURV(  M ,  N  N  +  J  0  )  =  CURV(  1 ,  NN  + JO  )  +  (  CUR  V  (  *’0,  NN  + JO  ) -CURV  {  1 
1NN+JQ ) ) *  A  M 

X  =  DNSTAR(M)/DELN(M)-FLCATF(NX*JO) 

PCURV ( M  )  =  CURV(M,NN)-B*(CURV(M,NN+JC)-CURV(M,NN)) *X 
950  CONTINUE 

PRINT  953,K,NX,D2TCN2,D2TDNM,D1(M),D2{N)  ,  PCURV (M) 

953  FORMAT (  12 ,  15,  5F15.4// ) 

40C  CONTINUE 
END 
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Subroutine  RELVEL 


FLOW 


END 


SUBROUTINE  RELVEL 

DIMENSION  DLAM(9,  10) ,R(9, 1C) ,Z(9, 1C), THETA(9,  1C)  , 

1 ALAM! 9, 1C) , AT  FETA (9,  10),CELTA(9,l0) ,ACELTA(9, 1C) , 

2AR(  9)  t  PLAN.  (9)  ,  GAMMA  (  9  )  ,  PGAMMA  (9  ) ,  PP  ET  A  (  9  )  ,  OE  LN  (  9  ) , 
3PDELT  A ( 9 ) , D  N  S  T  A  R ( 9),DM(9),D1!9),D2(9),CURV(9,10), 

4PCURV (9) , Yl  (9  )  ,Y2 (9) , Y3( 9 ) , CHECN ( 9 ) , DFL( 9 ) ,WM (9 ) , 

5DWSG ( 9),LX(9),AX(9),FX(9),BETA(9),DCELTA(9),XLAM(9), 
6RX! 9)  ,WX( 9)  , XGAMM A ( 9 )  ,  XBE  TA ( 9 ) VGAKX(9 ) ,W(9 ) , XCEL(9 ) , 
7ZBETA (9 )  ,  G  M  L  (  9),YCELTA(9),  CWCC(9)  ,DWFUNC  (9),TKEGV(9) 
COMMON  M, CELN,ATHETA,NO,CTCM,PBETA,PDELTA,PGAMMA,PlAM, 

1 AR,R, CM, UNSTAR, AO E L T A , NOGC , AL AM , NS , 0 1 , C2 , CLR V , PCUR V , MC , 
2Y1 »  Y2  »  Y3  »  WM  ,  QDE  L*  CW SC , OME GA , DX , AX , GA , CC , 
3WX,RX,XGAMMA,XRETA,YDELTA 
PRINT  975 

975  FORMAT! 1CX, 1HM, 10X, 2HWM) 

AX (  1  )  =0. 

PRINT  973,M,WM( 1) 

DO  500  M=2,M0 

CX ( M )  =  DK ( M ) / CCSF ( PGAMMA (  M-  1  )  ) 

A X  (  M  )  =  AX  (  M-  1  )  +DX  (  M  ) 

CWSQ(M-I)  =  Y  3  (  M-  1  )  -  (  WM(  M- 1 ) *  «2 ) *Y  1 !  M—  1  ) 

1-WM( M-1 ) *GMEGA»Y2 ( M- 1 ) 

WSQ  =  WM (M-1 ) **  2  +CWSC(M-1 )*CX(M)/12. 

WM ( M )  =  SORT  F ( WSQ ) 

DWSQ(M)  =  Y  3 ( M ) -WSG  *  Y 1 (M )-WM( M) »0MCGA*Y2 (M)  * 

WSQ  =  W  M ( M - 1 ) *  *  2  +( (CWSQ(M-1)+0WSQ(M) )/2.)«CX(M)/12.. 

WM(M)  =  SQRTF ( WSQ ) 

PRINT  9  7  3  » M »  W  M ( M ) 

973  FORMAT!  I  1C,F15.U) 

50C  CONTINUE 

CALL  FLCW ( MO , C . ) 

ODE L  -  GC-CA 
PRINT  97R  , OCE L 

97U  FORMAT ( /8HQC-GA  =  ,F10.U) 

END 
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Si ' 


*  • 


Subrout j ne  F L CW 


SUBROUTINE  FLOW(MM.OELTAX) 

DIMENSION  0LAM(9, 10) ,R<9, 10) ,Z(9. 10) ,THETA(9, 10)  , 
1ALAM(9, 10) .ATHETA  (9 . 10  I  ..DELTA  (  9, 10  )  »  ADELT  A(  9,  10)  . 

2AR ( 9 )  , PL AMI  9) .GAMMA {9 ) , PG AMMA ( 9 ) , PBET A ( 9 ) , CELN ( 9 ) , 
3PDELTAI 9) . DNS  T  AR( 9).CM(9),D1(9),D2(9),CURV(9,10), 
4PCURV(9) ,Yl (9), Y2 (9) , Y3(9),DHECN( 9) ,DFL ( 9 ) . toM ( 9 ) . 
5DWSQ(S)  ,C.X(9) ,  A  X  ( 9 ) »  F  X ( 9 ) , BE T A ( 9 ) , DDE L T A ( 9 ) ,XLAM(9), 
6RX ( 9 )  ,WX ( 9) ,  XGAMM  A ( 9 ) , XB E TA ( 9 ) , G AMX ( 9 ) , W ( 9 ) , XCEL (9 ) . 
7ZBETA (9) ,GML( 9) ,Y CELT A ( 9 ) , OtoCC ( 9 ) , DWFUNC ( 9 ) , TKECVI 9) 
COMMON  M.DELN, ATHETA. NO, DTDM , P BE T A , PD EL TA , PGA  MM A  »PLAM 
1 AR. R. CM.DNSJAR.ADELTA.NOGO.ALAM.NS.Dl .D2.CLRV .PCLRV.M 
2YV,  Y2,Y3,WM,QDEL,  DtoSC, OMEGA,  DX,  AX,  CA,  CC  , 

3toX, RX  .XGAMMA.XBET A.YDELTA 
IF(DELTAX)  620,620,630 
620  JD=  0 

GO  TO  6l<0 
630  J D= 1 
640  CONTINUE 

RX(MM)  =  AR(MM)+DELTAX*COSF(PGAMMA(MM)+PLAM(MM) ) 
to  X  (  M  M  )  =  WM(MM)+(  CELTAX*DtoSQ( MM) ) / ( 24 . # toM ( MM)  ) 

XGAMM A ( MM )  =  PG  AM  MA ( MM) +DELT AX* ( PG AMM A ( MM )- PG AMMA 
.  1 (MM-1 ) )/DX(MM+JD) 

XBETA(MM)  =  PP,ETA(MM)+DELTAX*(PBETA(MM)-PBETA 
1  (MM-1  )  )  /OX ( MM+JD) 

YDELTA(MM)  =  POEL T A ( MM ) + C EL T AX* ( PDE LT A ( MM) - PD ELT A 
1 (MM-1 ) )/DX(MM+JD) 

FX(MM)  =  RX(MM)*WX(MM)«COSF(XGAMMA(MM)  ) 

CO  60C  M=  1 , MM- 1 

F  X ( M )  =  AR(M)*WM( M ) « CCSF ( PGAMMA ( M  )  ) 

600  CONTINUE 
GC  =  C. 

1  =  1 

606  IG  =  1+2 

IF (  IG-MM)  601 , 607 ,6C2 

607  AA  =  F*( I ) 

XI  =  C  X (  1  +  1  ) 

X2  =  X 1 +DX ( 1+2 ) 

PB  =  (X2*X2*FX( 1+ 1  )-  Xl«Xl  *FX( I+2)-AA*(X2**2-Xl*«2) ) 
1 / ( XI  * X2*X2-X2*X1* XI ) 

CC  = ( X2*(FX(I  +  1 )—  F  X ( I )  )-X  1 « ( FX ( I +2 )-F  X ( I ) ))/ 

1  ( X  1  *X  1*X2-X2*X2*X 1 ) 

XAX  =  AX(  1  +  21+CEL TAX  „  „  ' 

DO  =  (AA* (XAX-AXt  I  )  )  +  BB* ( XAX**2-AX(  I )  «*2)/2. 

1+CC*( XAX**3-AX( I) «*3)/3. ) / 1 44 • 

GO  TO  602 

60 1  AA  =  F  X (  1  ) 

XI  =  CX(  1  +  1  ) 

X2  =  X 1 +DX (  1  +  2  ) 

RB  =  ( X2*  X2*FX ( 1+ 1 ) -  X 1  *  X  1  *FX( I+2)-AA*( X2«*2-X1*«2) ) 
1/(X1*X2*X2-X2«X1*X1 ) 

CC  =( X2*(FX( 1  +  1  )  —  F  X ( I  ) )-Xl*(FX( I +2 ) -FX ( I)  ) )/ 

1 ( X 1  *  X 1*X2-X2*X2*X 1 ) 

dq  =( aa*( ax( i  +  2 ) -  a  x  c  i )  )+ee*  ( ax( 1+2 )*«2-ax( i ) **2)/2. 
1+CC*( AX( I+2)**3~AX( I  )**3)/3. ) / 1 44  • 

GO  TO  603 
A  =  FX (  I  ) 

XY  =  CX ( 1 +1  ) 

XX  =  CX ( I  +  1  )  +DELT  A  X 

CciF( aIxx  +  B«C 1 AX(  1+1 J+DELTAX )**2-(AX(  I))  **  2  )V  2  .  )  / 1  44  , 
CC  =  CC+LC+6. 2B31 6 
IF (  IG-MM  )  604,605 , 6C5 
1*1+2 

GO  TO  6C6 
CONTINUE 
END 
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PROGRAM  ROTCR  2 

DIMENSION  DLAM19,  10)  ,R19, 10) ,Z(9, 10) ,ABETA(9, 10) , 
1ALAM(9,10),BETA19,10).0ELTA19,10),ADELTA(9,10), 

2AR(9)  , PL AM  I  9) ,  GAM M A  1 9 ) , PG AMMA ( 9 ) , PBET  A ( 9 ) ,CELN(9), 
3PDELTA1 9) .DNSTARI 9),DM(9),D119),D2(9) ,CURV(9, 10)  , 

4 PCURV ( 9  )  , Y 1 (9)  ,  Y2  ( 9  ) , Y 3 ( 9 ) , DHEDNI 9 ) ,DFL ( 9 ) . WM ( 9 ) . 
SDWSQ(9),DX(9).,AX(9),FX19),APBETA(9)  , DOE  LT  A l 9 )  ,XLAM(9) , 
6RX l 9) , WXl 9 ) , XGAMM  A ( 9 )  , XBETA ( 9 ) , GAMX 1 9 ) , W 1 9 )  ,XDEL (9 ) . 
7ZPETA ( 9) , GML1 9) , YDFLTAl 9)  ,DWCD( 9) , DWFUNC (9)  ,TKECV( 9) 
COMMON  M, DCLN, ARE  T A, NC.PRET A, POELTA, PGAMMA.PL AM, 

1 AR, R.DM.ONSTAR, ADELTA.NOGO, ALAM.NS, 01 , C2.CURV , PCURV, MO 
.  2Y1 ,Y2  ,Y2,WM,C0EL, DW SC , OMEGA , CX , AX , QA , CC , 
3WX,RX,XGAMMA,XPET  A, YDELTA 
READ  1C, MG, NC. NS 

11  FORMA.Tl  6F  1C.0) 

10  FORMAT ( 3110) 

READ  11  ,  (OELN(M),  M= 1 , MO) 

READ  1 1 , 1 1 R ( M, N ) , N= 1 , NO ) , M= 1 , MO) 

READ  1  1 , (  (Z(M, N),N=1 .NO) ,  M  =  1 ,M0) 

READ  11,  (  ( DL AM ( M, N ) , N  » 1 , NO ) , M= 1  , MO ) 

READ  1  1  ,  ( BETA ( 1 ,N ) , N=1 ,N0) 

READ  1  1  ,  l  OETAl 9,N ) ,N=1 ,N0 ) 

READ  1  1 ,  (  ( DELTA (M , N), N=1 ,NC) ,M= 1  ,  MO) 

READ  11,  ( CUR V ( 1 ,  N ) , N  =  1 , NO ) 

READ  11,  (CURVI9, N),N=1,N0) 

READ  90, OMEGA, GA 
90  FORMAT ( 2F 1 5.0 ) 

READ  11,1 DHEDNI  M) ,M=1 ,M0 ) 

READ  11 , ( DFL( M) ,M=1 ,M0) 

READ  90, THICK, RLNO 
DO  12  M= 1 , MO  ' 

CO  12  N= 1 , NO 
C  =  3.14159/18C. 

ABETA(M.N)  =  BETA ( M , N  )  *C 
AL AMI M, N )  =  DLAM(  M , N ) «  C 
ADELTA  1 M , N  )  =  OELTA(M,N)*C 

12  CONTINUE  '  ' 

DO  29  M  =  2, MO-1 

AM  =  FL0ATF1M-1  )/FLCATF( MO-11 


DO  29  N= 1 »  NO 

A BE  TAIN, N)  =  ABET Al 1 , N) + 1 A8ETA 1 9, N )- ABET A ( 1 ,N ) ) « AM 
29  CONTINUE 

CO  997  M  = 1 , MO 
CO  997  N  = 1 , NO 
BETA!  M.N )  =  ABETA ( M.N  )  /C 
997  CONTINUE 
M  =  1 

AR(M)  =  R(M.NS) 

PCURV ( M l  =*  CURV1M.NS) 

PDELTA(M)  =  ADELT  A { M , NS ) 

PBETA(M)  =  ABETAi M,NS  > 

TGAMMA=(SINF(PBETA1H) )•*?)•  TANFl POEL  TA ( M ) ) 
PGAM.MA1MJ  =  ATANF  (  T  GAMMA  1 
GAMMA (M)  *  PG  AMMA ( M ) /C 
DNSTARI 1)  =C. 

DM ( 1 }  =  0. 

PLANT  Ml  =  A  L  A  M.  {  M » N  S  ) 

APRFTA(M)  =  PBETA (M)  /C 
CDCLTA(M)  =  PDELT  A(M)/C 
XLAM1N)  =  PLAMIMJ/C 
G M L 1 M  )  =  GAVMA(W) -XLAM(M) 

CO  16  f1  =  2  ,  MQ 
NtlOO  *  C 

’  \P  =  CNS1  AR  (M- U/ CcLM  M- 1  ) 

NE  -  NS+NP  „  , 

lFfPLAM(M-l)-.7854)  17,17,18 

17  Af)M  •=  (  R{M,NE  )-R(  M-  1  ,  NE)  )  /COSFl  PLAM1M-1  1  ) 

C  V  (  M  )  -  A  it  $  F  (  A  ■  v ) 

1C  ADM7-  (Z(M. M  )-?( M- 1 , ■ 1 ) )/SINF(PLAM(M-1 ) ) 
i  ..*  (  A  )  =  1  H  SF  (  At  *'  I 

19  D  N  S  T  • ,  R  t  M  )  --  OI;STTR*.H  1  >  +  I  Ml  M  )  »  T  ANF  1  P«J  AM  M  A  (  M-  I  )  ) 

I  H  ~ N  1  f  A  IN))  .  •  - 


2C 

22 

21 


24 

25 
23 

26 

27 

28 


16 


31 


32 


33 

34 

35 


36 

37 


38 


4C 


CALL  FSTAR(NS,-1,-1.) 

GO  TO  23 

CALL  PSTAR ( NC-NS,  1,1.) 

GO  TO  23 

PBETA(M)  =  ARETA( M,NS ) 

PLAM(M)  =  ALAM(M,  NS) 

PCURV(M)  =  CURV (M , NS ) 

PDELTA(H)  =  A CELT  A  ( M  ,  N  S  ) 

AR(M)  =  R(M,NS) 

T G A  N N  A=  (  S  I N F  (  P 8  ET  A  (  M )  )**2  )*  T ANF ( PDELTA ( M ) ) 
PGANMA(M)  =  AT ANF ( TGAMMA) 

I  F  f  PGAMKA  (M  )  )  24,  23,25 
CALL  PSTAR( NSt-1, -1  .  ) 

GO  TO  23 

CALL  FST  AR ( NC-NS  ,1,1.) 

CONT  INUE 

IF(I-NOGG)  26,26,28 
PRINT  27, M 

FORNATf 27HINSUFFI CIENT  GRID  WIDTH  M=,I2) 

GO  TO  55 

GAMMA (N)  =  PGAMMA ( M)/C 
APBETA(M)  =  PRETA ( M ) /C 
DDELTA  (  M)  =  PDELT  A  ( M. )  /C 
XL  AN { N  )  =  PLAN ( M )  /C 


GML(M)  =  GANNA (N)-XLAM(M)  .  * 

CONTINUE 

CALL  DDL 

CONTINUE 

DC  31  M= 1 ,NC 

V 1 ( M )  =  CCSF ( PG AMN A ( M ) )*(2.*PCURV(M)*(COSF(PBETA(M) ) ) 
1  +  SINF(2.«PEETA(M))*C2(N>))«12. 

Y2(M)  =  2.*COSF(PGANNA(K) )  *S  INF { 2  .*PBETA (M) ) »D1 (N) 

1 /COSF ( PDELT  A ( N  )  ) 

Y3(M)  =  (DHEr>N(N)-CNEGA*OFL(M)  )*COSF(PLAN(M)  )*2.* 
1COSFI PGANMA(M) )«•(  CCSF ( PPETA { M ) ) ) **2 
CONTINUE 
EE  =  .COI 
WM(  1  )  =  o. 


I  =  1 

WDEL  ' 
WM  (  1  ) 
CALL 


:  ICO. *(-.!**( 1-1  ) ) 

=  W  K  (  1  ) +WDEL 
RELVEL 

I F ( ABSFICDFL  )-EE  )  30,38,3: 
IF(CCEL)  32,38,34 
I  =  1  +  1 

WDEL  =  1CC.*(-.1»*( 1-1 ) ) 

W M  (  1  )  =  WM  (  1  J+WDEL 

CALL  RELVEL  „  ,  , 

I  F  (  ABSF  (QDEL  )  —  EE)  38,38,  36 

IF(QDEL)  37,38,35 

I  =  1  +  1 

GO  TO  32 

CONTINUE 

WX(  1 )  =  WM (  1  ) 

WX(NO)  =  WM(NO) 

RX ( 1 )  =  ARID 
RX(NO)  =  AR(NO) 

XRETA (  1  )  =  PPLTA(  1  ) 

=  PRETA (MO) 

=  PGAMMA  (  1  ) 
a  PGAMMA(MG) 

=  PDELT A (  1  ) 

=  PDELTA (NO) 


XBETA (NO) 
XGAMMA  (  1  ) 
XGANMA ( NO  ) 
YDEL1 A (  1  ) 
YDELTA(MO) 
NO  1 =  M0- 1 
DO  50  K  =  2 , 


N0 1 


KK  =  NC+ 1 -K 
F  =  F  LC A  I F ( NO-K ) / FLO  A  TF ( NO- 1 ) 
OAF  =  0  A  *  F 
EE  =  .OCI 

CALL  FLOW (KK,C. )  c  ^  „ 

I F ( ABSF(QC-CAF)-EE)  5C,5C,40 
XDEL(KK)  =  .2 
ODD  =  XOi-L  (  KK  ) 


*  *2 
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41  CALL  FLCMKK.CDD) 

IF(QC-QAF)  42,50,43 

42  XOEL(KK)  =  XOEL(KK)+.l 
ODD  =  XDEL(KK) 

C...  ORIVES  OC  GREATER  THAN  QA  TO  START 
GO  TO  41 

43  J  =  1  < 

44  XDEL(KK)  =  XOEL <K K ) 1 • ( - . 1 •• ( J-l ) ) 

ODD  =  XDEL(KK) 

CALL  FLOW ( KK , ODDI 

IF ( ABSF (GC-GAF  )-E  E)  50,50,45 

45  IF(QC-GAF)  46,50,44  t  * 

46  J  =  J4l 

47  XDEL(KK)  =  XOEL! K K ) - . 1  * ( - . 1 •• ( J- 1 )) 

CDD=XCEL ( KK ) 

CALL  FLOW ( KK , ODD) 

IFCABSFf QC-QAF )-E  E )  50,50,48 

48  IF(CC-CAF)  47,50,49 

49  J  =  J  4 1 
GO  TO  44 

50  CONTINUE 

CO  51  M=  1  »  MC 
ZBETA(M)  =  XBETA! Ml/C 
GAMX(M)  = -PGAMMA!  M)/C 

W ( M 1  =  WX(M)*SQRTF(1.4(TANF(XBETA(M))««2)) 

TKEGVIM)  =(THICK/CCSF(XBETA(N) ) ) • ( SQRTF ( 1.4 (S INF 
1 (XBETA! M) 1«TANF(Y DELTA! Ml  1  ) • *2 ) 1 
CWCC(M)  =  (6.2832/8LN0-TKEQV(M)/RX(M) )«C0SF(XBETA(M1 1 
DWFUNC(M)  =  RX(M1 «WX ( M ) »T ANF ( XB ETA ( M ) 1* 1 2 • 40M EGA»R X ( M) 
1 *RX ( M  > 

51  CONTINUE  . 

PRINT  R8C 

880  FORMAT ( 1H1 /////// /25X9H  TABLE  ///17X 
127H  DATA  FOR  CHARACTERISTIC  C9/// 1 
C  CHANGE  CHAR.  NO.  F CR  EACH  COMPUTATION 

882  FORMAT ( 6H  M  =  8X 2H  19X1H2  9X1H3  9X1H4  9X1H5//1 
PRIM  883 

883  FORMAT! 12X33H  LOCATION  OF  CHARACTERISTIC  CURVE//) 

'  PRINT  884, NS 

884  FORMAT ( 23H  STARTING  NORMAL  NO.  =  ,12//) 

PRINT  882 

PRINT  886, (CNSTAR(M),M=1,5) 

886  FORMAT (  1 CH  CNSTAR  =,5F1C.4) 

PRINT  P97,!GML!M) ,M=1,5) 

897  FORMAT ( 1 CH  GAM-LAM  -.5F10.4) 

PRINT  885, (GAMMA! Ml ,M=1 , 5) 

885  FORMA  T ! 1 CH  GAMMA  =,5F1C.4) 

PRINT  Pe7,  (APBETA(M»),M=1 .5) 

887  FORMAT ( 1  OH  BETA(P)  =,5F10.4) 

PRINT  P8E , ( AR ( M) , M= 1 , 5 1 

888  FORMAT ( 1CH  RACIUS  =,5F1C.4///1 
PRINT  889 

889  FORMAT ( 16X28H  LOCATION  OF  NEW  STREAMLINES//) 

PRINT  882 

XOEL!  1  )  =  C . 

XDEL! MO) =0.  .... 

PRINT  890,(XDEL(M),M*1,5) 

890  FORMAT ( 1 CH  CELTA  X  =,5F10.4) 

PRINT  891 , ( ZBETA! M) ,M=1 .5  ) 

891  FORMAT ( 1CH  BETA(X)  =,5Fl0.4) 

PRINT  888, (RX(M),M=1 ,5) 

PRINT  892 

092  Format (2i xi7H  velocity  prcfile//) 

PRINT  882 

PRINT  893, (WM(M),M=1 ,5) 

893  FORMAT ( 1 CH  WM(P)  =,5F1C.4) 

PRINT  894 , ( WX( M), M= 1 ,5) 

894  FORMAT!  1CH  WM(X)  =,5F1C.4) 

PRINT  895 , ( W( M) ,M= 1 ,5 ) 

895  FORMAT!  1CH  REL  VE L  *,5F]0.4) 

PRINT  87  6,  (DWCO(M  )  ,M=1,5) 

875  FORMAT! 1CH  CW  COEF  =,5F1C.4) 

PRINT  874 , (CWFUNC (M) ,M=1 , 5) 
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87U  FORMAT I 10H  CW  FUNC  =,5FIC.O) 

896  FORMAT  ( 6R  M  =  8X2H  69X1H7  -9X1H8  9X1H9//) 
PRINT  88C 
PRINT  883 
PRINT  884, NS 
PRINT  896 

PRINT  886, (CNSTAR (M) ,M=6,9> 

PRINT  897, (GML(M) ,M=6,9) 

PRINT  885, (GAMMA! M) ,M=6, 9 ) 

PRINT  887, ( APRETA (M) ,M  =  6, 9) 

PRINT  878, ( AR ( M  ) ,  N=  6 , 9  ) 

87e  FORMAT ( ICH  RADIUS  =,4F10.4///) 

PRINT  889 
PRINT  896 

PRINT  89C,(X0EL!M ) ,  M  =  6, 9 ) 

PRINT  891,(ZRETA(M),M=6,9) 

PRINT  P78 , ( RX ( M ) , M=6 , 9 ) 

PRINT  892 
PRINT  896 

PRINT  893t ( WM(M),  M  =  6, 9) 

PRINT  P94 , (WXIM),  M=6,9) 

PRINT  895, (W(M) ,M=6,9) 

PRINT  875, (DWCCIM ) ,M=6,9) 

PRINT  874 , ( CWFUNC  (M  )  , M=6 f9) 

PRINT  86C 

860  FORMAT (1H1 ,4H  END  J 
55  CONTINUE 
ENC 


134- 


Proa  ram  ROTOR  2 


Subroutine  PS  TAR 
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SUBROUTINE  PST AR( NN , JO , B ) 

CIKENSICN  CLAM {9,  1C),R(9,1C),Z(9,10),A8ETA(9,10), 

1  AL  AM!  (  9  ,  1 C  )  ,  8  C  T  A  {  9  ,  1C)  ,  DELTA (9, 10) ,ADELTA(9,1C), 

2AR( 9) ,PLAM(9) ,  GAMMA (9) , PG AMM  A ( 9 ),PBETA{9) »  DELN ( 9 ) » 
3PDELTA ( 9) ,0NSTAR(9),CM{9),01(9),02(9) ,CURV( 9,  10 ) , 

4  PCUR  V  (  9  )  t  VI  (9  )  tY2  (9  ) ,  Y3(  9  )  ,  D  H  F  D  N  (  9  ) ,  CFL  (  9  )  ,WM  (  9  ) , 

5DWSC(  9  )  ,DX  (9)  ,  AXI  9)  ,FX(9)  ,APP.ETA(  9)  ,0CELTA(9)  ,XLAM  (9)  , 
6R  X ( 9 )  » W  X ( 9 ) »  X  G  A  MM  A ( 9 )  ,XBETA(9) , GAMX (9 ) , M 9 ) , XCEL (9 ) , 

7 ZB ET A ( 9 ) »  GML ( 9 ) ,Y DELTA! 9) ,CWCG(9) ,OkFUNC(9) ,  TKEGV(  9) 
COMMON  M,DELN  ,  ABE  TAt  NC,  PB  ETA  ,  POE  LTA ,  PGAMM.A  ,PLAM, 

1  AR,R,CM,CNSTAR, ADELTA ,NOGC, ALAM,NS,Bl  ,  C2  , CURV , PCUR V , MO , 
2Y1 ,Y2,Y3,UM,GCEL, CWSG , OMEGA, CX , AX , G A, CC , 

3WX,RX  ,XGAMMA, XBET A,YDELTA 
DO  30C  1=1,2 
DEL  =  CELN  {  M. ) 

DO  301  K  =  1  , N N 
NA=NS* JO*  C  K— 1 ) 


NB  =  NS+JC*K 

X  =  ( CNSTAR(M)-B*DEL) /DELN(M) 

I F { ABSF(ONSTAR(K)  ) -CEL) 302,303,304 

304  CEL  =  DELN(M)»FL0ATF(K+1 ) 

GO  TO  3C1 

302  AM  =  FLGATFl M-l )/ FLOATF( MC- 1 ) 

CURV(M,NA)  =  CURV ( 1 , NA )+(CURV (MO,NA)-CURV ( 1 ,NA) ) «AM 
CUR V ( M, NC )  =  CURV( 1 , NB ) + ( CURV ( MO, NB l-CURV ( 1 , N8 ) ) «AM  * 
PCURV(M)  =  CURV(M ,NB )-B« (CURV( M, NA ) -CURV ( M ,NA ) )*X 
PBETA(M)  =  ABETA(M»NH)-B*(ABETA(M,NA)-ABETA(M,NB))*X 
AR ( M  )  =  R(K,NB)-B«(R(M,NA)-R(M,NB) )*X 

PDELTA(K)  =  AOELT  A(M,NB)-B*(ACELTA(M,NA)-AOELTA(M,  NB) )*X 
GO  TO  3C5 

303  PBETA(M)  =  ARETA(M,NB) 

AR ( M  )  =  R(M,NH) 

PDELTA(M)  =  ADELTA(M.NB) 

AM  =  FLCATF (M-l )/ FLOATF( MC-1 ) 

PCUR V ( M)  =  CURV ( 1  ,  N6>)  +  (  CURV  ( MO , NB  )-CURV  (  1  »N8  )  ) * AM 

305  CONTINUE 

TGAMMA= (S INF ( PBET  A ( M ) )  **2  )  *  T AN F ( PDEL T A ( M )  ) 

PGAMMA(M)  =  ATANF ( TGAMMA ) 

PLAM(M)  =  ALAM (M, NP  )-B*( ALAMIM,  NA )-ALAM(M,NB ) ) «X 
I  F (  I-  1  )  306,306,300 

306  TOAMD  .=  T  AN  F  (  PG  AM  M  A  (  M  )  )  -  T  ANF  (  PG  AMMA  ( M  - 1  )  ) 

IF(TGAMC)  3  C  7 , 3  CO , 3  C  7 

3C7  EN  =  (  D  M  ( »M  )  *•  T  G  A  MO  )  /  2  . 

CNSTAR(M)  =  D  N  S  T  AR ( M )  +  E  N 


301 


GO  TO  3C0 
CONTINUE 
NOGO  =  I 


300  CONTINUE 
END 


13  6. 


Subroutine  DDL 
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subroutine  ccl 

DIMENSION  DLAM.(9,  1C)  ,R(9,  10)  ,Z(  9,  10)  ,  ARETA(9,  1C)  , 

1 ALAM ( 9 , 1C)  , eETA(9 , 1C) ,DELTA(9, 10) , ADELTA (9, 10), 

2AR  ( 9  )  »  P  L  A  M  (  9  )  ,GAMMA(9)»PGAMMA(9),PBETA(9)  ,CELM9), 
3PDELT  A ( 9 )  , CNST AR( 9 ) , CM ( 9 ) , D 1 ( 9 ) , D2 ( 9 ) , CURV ( 9 ,  1 C ) , 

4PCURV (9 ) , Y1 (9 ) , Y2 (9) , Y3 ( 9 ) , OH EON ( 9 ) ,DFL(9),NM(9), 

5QWSGI 9), CXI  9) , AX( 9 ) , F X ( 9 ) , APC ET A ( 9 ) , CCELTA (9 ) , XL AM ( 9 ) , 

6  R  X  (  9  ) »  to  X  (  9  ) ,  X  G  A  MM  A  (  9  ) ,XBETA(9) ,GAMX{9) ,W(9),XUEL(9), 
7ZBET  A ( 9 ) » GVL ( 9 ) , YDELT  A (9 ) ,  DWCC ( 9 ) , DtoFUNC (9),TKECV(9) 
COMMON  M,DELN, ABE  T A, NO, P PET A, PC ELT A , ?G AMM A , PL  AM , 

1  AR,R, CM.CNSTAR, ADELTA,NCGC,ALAM,NS, 01 ,C2,CLRV,PCURV,MC  , 
2Y1 , Y2,Y3,WM,QDEL,  CtoSC , OMEGA , UX , AX, CA , CC , 

3WX , RX ,XGAMMA, XB ETA, YD ELT A 
CO  40C  M= 1 , MO 
NX  =  CNSTAR(K)/DELN(M) 

IF(M-MC)  401,402,400 

4  C 1  MM  =  M+  1 
MX  =  MM 
EP=1 . 

GO  TO  403 

4C2  MM=M- 1 
MX  =  M 
8B=-1 . 

403  I F ( DNST  AR { M ) )  404,408,405 

404  N A=NS  +NX 
N8=NA- 1 

X  =  -DNSTAR(M)+FLCATF(NX)*DELN(M) 

DEL  =  1  . -X/DELN  (M ) 

GO  TO  4  C6 

405  NB»NS+NX 
NA=NB ♦ 1 

X  =  ONSTAR(M)-FLOATF(NX)*CELN(M) 

DEL  =  X /C  ELN ( M ) 

406  IF ( X )  407 ,408,407 

407  ADN  »  (TANF(ABETA(MM,NA)  ) -TANF ( ABETA ( M , N A  )  )  )  /  CM  (  MX  )  *8  P. 
BUN  =  (TANF (ABETA  (MM, NR)  )-TANF(ABETA(M»NB) ) )  / CM ( MX ) «P6 
DTCN  =  BDN+ ( ADN-B  DN  )  «CEL 

DTBM  =  (TANF( A6ETA(M,NA) J-TANF (ABETA( M,NE) ) )/CELN(M) 

GO  TO  4C9 

408  NC  =  NS-»NX 

DTRN  =  (7ANF(ABET A( MM,NC ) )-TANF( ABETA (M,NC) ) ) /DM(MX) *BB 

NA  =  NS-»NX+  1 

NB=NS+NX-1 

DTBM  =  (TANF(AEETA(M,NA) J-TANF ( ABETA( M,NB) ) )/ 

1 ( DELN (M) «2  .  ) 

4 09  CONTINUE 

D 1  ( M)  =  CCSF(PDELTA(M)-PLAM(M) ) 

D2(M)  =  TANF( PDEL TA(M ) ) *0 TPM+ CTBN+ ( D 1 (M) *TANF ( PBET A 
KM)  )  )  /  (AR(M)*CCSF  (PCELTA(M)  )  ) 

4 OC  CONTINUE 
END 
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Subroutine  RELVEL 


S.R.  FLON 
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SUBROUTINE  RELVEL 

DIMENSION  C  L  A  N  (  9 , 1 C  ) , R( 9, 1C ) , Z ( 9. 1 0 > , ABETA (9, 10) , 

1  ALAN ( 9  t 1C ) , BETA (9  t 1C) ,DELTA( 9, 10) * ADELT  A  ( 9  ♦ 1C ) t 
2AR ( 9 )  ,PLAN(9) ,GAMNA(9  ),PGAMMA(9) ,PBETA( 9),CELN(9), 
3P0ELTA19) ,DNSTAR( 9),DM(9) ,C1 (9) ,02(9) ,CURV(9, 1C) , 

4  PCURV (  9  )  ,  Y  1 (9) ,Y2 (9),Y3(9),DHEDN(9),0FL(9),WM(9>, 
5CWSC(9),CX(9),AX(9> ,FX(9),APBETA(9) ,0DELTA(9>  ,  XL  AM  (9) 
6RX ( 9 ) ,WX(9) »XGANMA(9),XBETA(9J , GAKX ( 9 ) , W ( 9 ) ,XDEL(9 ), 
7ZB6TA (9)  ,GML(9) ,YCELTA(9)  ,OWCO (9) ,OWFUNC (9 ) , T KECVl 9 > 
COMNCN  M,OELN, ABE TA , NO, PBET A , PCELTA , PGANNA, PL AM, 

1  AR,R,CM,CNSTAR,ADELTA ,NQGO, ALAK,NS,D1 , D2 » CLR V , PCUR V , 
2Y1 ,Y2  ,Y3,UVtQCEL, CWSC , OMEGA, CX, AX ,CA , CC , 
3UX,RX,XGAKMA,XBETA, YDELTA 
A  X  t 1 )  =C . 

CO  500  M=2,MC 

DX ( M )  =OM(M)/CCSF ( PGAMMA ( M- 1 ) ) 

AX(M)  =  AX(N-1 )+OX(N) 

DWSC(M-l)  =  Y3  (  M-  1  )  -  (  toiX  (  N  - 1  )  ** 2 ) * Y 1 (  M-  1  ) 

1  -  W  N  (  M-  1  )  *  CN.EQA »Y2  (  M-  1  ) 

WSC  =  VM(M-1)«*2  +C'aSC(M-1  )«DX(M)/12. 

WM(M)  =  SCRTF(WSQ) 

CWSQ(M)  =  Y3(M)-WSC*Yl 1 M)-WM( M) *0MEGA«Y2 ( M) 

WSQ  =  hM(N-1)**2  +( (DWS0(M-1 J4DWSCIN)  )/2.  )«DX(N)/1 2. 
WM(M)  =  SQRTF(KSQ) 

5 OC  CONTINUE 

CALL  FLOMMC,C.) 

QDEL  *  CC-QA 
END 


14-0 


Subrout jnm  FLCM 


141 


600 

6C6 

607 


01 


IF IDELTAX  ) 

JD  =  G 
GO  TO  640 
JU=  1 

CONTINUE 

RX(MM)  =  AR(^M)+CELTAX*COSF(PGAMKA(MM)-PLAN(MM) ) 
WX(MM)  =  to.Y(MM)+l  CELTAX«CWSQIN’M)  )/I24.*WfMMM)  ) 
XGANMAIKM)  =  PGAMN'A(I^M)  +  DELTAX*(PGAMMA(MM)-PGANPA 

1  inm-i )  j/cximw+jo) 

XBETAlf'M)  =  PEETA (MM) +DELT AX* l PBET A { NM )- PB ET A 
1 (MM- 1 ) ) /CX( MM+JO) 

YOELT  A  ( )  =  PDELTA(MM)  +  nELTAX*{PDELTA(NM)-Pl)ELTA 
1 ( MM- 1 ) )/CX(MM+JD) 

FX(MM)  =  RX(mm)*WX(MM)«CCSF(XGAMMA(KM  )) 

DO  6CC  M=1,NM-1 

FXIM)  =  A R {  M)  *KM(  M )  *COSF  (  PGAPNA  (  M  )  J 

CONTINUE 

CC  =  0. 

I  =  I 
IQ  =  1+2 

I F  C IG-MM)  601  *607,602 
AA  =  F  X ( I  ) 

xi  =  oxi  i  +  n 

X2  =  Xl+DXI 1+2) 

RB  =  IX2*X2*FXII+1)-  X1*X1  *FX( I +2 ) -A  A * l X2* * 2-X 1 **2) ) 
1/{Xl*X2*X2-X2*Xl*Xl ) 

CC  =1X2*1 F  X { 1  +  1  ) —  F  X { I ) )  —  X  1 »  l  F  X I I  +  2 ) -F  X  I I ) ) )/ 

1  { X  1  *X  1*X2-X2*X2*X 1 ) 

XAX  =  AX ( I +2)  +  CEL  TAX 

DQ  =  ( AA* ( XAX-AXI  I  )  )  +B5* IXAX**2-AX( I)  *  *  2  )  /2 . 

1+CC*  I  XAX**3-AXI I) *«  3 ) /3. ) /144 . 

GO  TO  6C3 
AA  =  FX(I) 

XI  *  CXI  1  +  1) 

X2  =  X 1+DXI 1+2  ) 

RB  =  ( X2*X2*FX  (  1+ 1  )-  Xl*Xl  *FX  (  I  +2  )-A  A*  I  X2«  *2 -X  1  «*  <; )  ) 
1 / ( XI* XT* X2-X?*X  1* X  1  ) 

CC  =  I  X  2  *  I  FX(  1  +  1  )-FX(  I  )  ).-X1*IFX(  I+21-FXI  I  )  )  )./ 
i (xi»xi*x2-xa«x2«x.r)  ‘  ■ .  .  .  •  „ 

DO  =  { A  A  * ( AX (  I +2 )- AX  I  I)  )  +  P6*  (AXI  I  ♦2.')  *  *  c-  A  X  ( I  )  *  *  2  )  /  2  . 
1+CC*  (  AX  (1+2)  *«3r-A-X(  I  )  **3 ) /3. ) /  1  44« 

GO  TO  603 
A  =  F  X  I  I  ) 

xy  =  cxi  i  +  n 

XX  =  CX  I  I  +  n+UELT  AX 

00=^ I  A  *XX  +  PM  I  AXI  I  +  1)+CELTAX)**2-IAX(  I))**2)/2.)/144. 

or  s  cc  +  i:c*6.  ?C3l  p 


cc  =  cc  +  i;c*6.? 

IFIIC-f-T)  604, 
1=1+2 
GO  TC  6C6 
CONTINUE 
END 
END 


P3i 
6  C  5 , 6  C  5 
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PROGRAM  BLADE 


VARIABLE  NAMES 

Name  Equivalent  to  Units 


DELN 

AN 

in 

DWCO 

DWcoef 

DWFUNC 

DWfunc 

DDM 

d(DWfunc)/dm 

W 

w 

ft. /sec 

WDSL 

AW 

ft. /sec 

WSUC 

ws 

ft. /sec 

WPR3SS 

WP 

ft. /sec 

INDEX  NAMES 


MO  s  Number  of  streamlines 

NMAX  3  Maximum  number  of  data  points 

NN  s  Number  of  data  point  on  a  streamline 
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Program  bL*0£ 
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PROGRAM  BLADE 

C I  MENS  ION  DWCCI  10  ,3C ) .DWFUNC ( 1C . 30 1 . W C 1C, 30 > .NN ( 10 ) , 

1 WDEL(  1C,30),WSUC(  1 C , 3C 1 , WPRESSf 10,30) ,OELN( 10  I  ,00M( 30 ) 
READ  10.HO.NMAX 
10  FORMAT ( 6 I  10) 

READ  10, (NN(M) ,N= 1 ,MC ) 

READ  20, (OELN(M), H=1,M0) 

20  FORMAT ( 6F 10.0) 

RE  AO  2C,( (DwCO(M. N),N=1,NMAX) ,M«1,M0) 

RE AC  20, ( ( DWFUNC ( M ,N  )  ,N= 1 , NMAX ) . M» 1,M0) 

READ  20, ( (W(M,N),N«1,NMAX),M*1,M0I 

CO  30  M=1,M0 
PRINT  40 

40  FORMAT! 1H 1 /////// /25X9H  TABLE  ///21X 
118H  VELOCITY  PROFILE  ///» 

PRINT  50, M 

50  FORMAT ( 16X.26H  MERICICNAL  STREAMLINE  NO. ,12///) 

PRINT  60 

60  FORMAT ( 4H  M ICXl HW 1 4X2HDW8X9HW  SUCTIGN4X1CHW  PRESSURE// ) 

DELTA  =  12.*DELN(  M) 

NM  =  NN ( M ) 

DO  30  N=  1  • NM 
I  F  (N-2  )4.1 , 4  1  ,42 

41  DDM(N)  =  (-25.«OWFUNC(M,N)+46.*DWFUNC(M,N*l 1-36. • 

1  DWFUNC  ( M, N*2  )♦  16.  •  DWFUNC  ( M, NO  >-3.«0WFUNC  (  M  ,N>4  )  l/DELTA 
GO  TO  45 

42  NRFM  =  NM-N 

IF (NR EM- 1 >43,43,44 

43  COM ( N )  =  (3.*DWFUNC(M,N-4)-16.»DWFUNC(M,N-3>+36.« 

1 DWFUNC (M.N-2) -48.  •DWFUNC ( H, N- 1 ) *2 5. »DWFUNC ( M, N ) ) /DELTA 

GO  TO  45 

44  DOM ( N  )  =  (DWFUNC(  M, N-2 > -DWFUNC ( M, N+ 2 ) -9 .•( DWFUNC ( M ,N- 1 ) 

1 -DWFUNC ( M, N* 1 )))/ DELTA 

45  CONTINUE 

X  =  DELN(M) «FLOATF(N- 1 ) 

WDEL(M,N)  *  DWCO( M»N)#DDM(N)/24. 

WSUC(M.N)  *  W(M,N ) ♦WDEL( M,N) 

WPRESS(M.N)  =  W(M .N)-WDEL(M.N) 

PRINT  100,X,W(M,N ).WOEL(M,N),WSUC(M,N),WPRESS(M,N) 

100  F0RMAT(F6.2,4F14.4) 

3C  CONTINUE 
PRINT  70 

70  FORMAT ( 1H1.5H  END) 

END 

END 
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